Đề tài Ứng Dụng lý thuyết độ tin cậy và phương pháp thiết kế ngẫu nhiên trong đánh giá an toàn ổn định đê kè biển

kiệt ở phía cửa Nam Triệu hầu như không biến động mạnh. Bảng 2. Sai số quân phương tương ứng với sự khác nhau của thuỷ triều hình thái N 0 Con triều Hình dạng triều RMSE (Mùa lũ) N0 Con triều Hình dạng triều RMSE (Mùa kiệt) 1 08 11 : 09 12 0.1561 1 05 18 : 06 19 0.0687 2 09 12 : 10 13 0.0858 2 16 18 : 17 19 0.0491 3 23 12 : 00 14 0.0977 3 06 19 : 07 20 0.0507 4 10 13 : 11 14 0.0766 4 17 19 : 18 20 0.0415 5 11 14 : 12 15 0.0892 5 06 20 : 07 21 0.0428 6 23 13 : 00 15 0.1096 6 18 20 : 19 21 0.0283 7 00 15 : 01 16 0.1124 7 07 21 : 08 22 0.0377 8 19 21 : 20 22 0.0505 RMSE = 0.0767 xây dựng công trình thủy lợi 18 Hình 5a: Kết quả mô phỏng thuỷ triều hình thái mùa lũ Hình 5b: Kết quả mô phỏng thuỷ triều hình thái mùa kiệt IV. Nhận xét và kết luận - Từ kết quả của mô hình cho thấy sự phản ứng của hình thái và kiểu chuyển tải bùn cát có liên hệ mật thiết với thuỷ triều, đặc biệt là biên độ và hình dạng triều, trong khi góc pha triều không ảnh hưởng đến chúng lắm. - Trong trường hợp này thuỷ triều hình thái ứng với triều cường và triều kém đều cho sai số lớn hơn so với triều trung bình. Con triều cho kết quả gần giống với mô phỏng cho cả mùa có đỉnh triều xuất hiện trước chân chiều. - ở đây, mô hình không ió và sự biến đổi độ mặn trong nước, mà chỉ qu Kết quả mô phỏng của các ở đây. Tuy nhiên đòi hỏi p cũng như chạy với nhiều k quả tốt hơn. - Từ những kết quả trê sông ven biển Hải Phòng, (Xem tiếp trang 36) MorphologicalMùa kiệt RMSE = 0.0284 Hình 5a: Kết quả mô phỏng thuỷ triều hình thái mùa lũ được xem xét với sự ảnh hưởng hỗn hợp của sóng, g an tâm đến ảnh hưởng của thuỷ triều, dòng triều và dòng chảy trong sông. mô hình đối với vùng Hải Phòng cho thấy nó là công cụ tốt cho nghiên cứu hải có những nghiên cứu đo đạc đầy đủ địa hình và số liệu quan trắc khác ịch bản khác nhau để kiểm chứng thuỷ triều hình thái chắc chắn sẽ cho kết n có thể thấy được khả năng ứng dụng mô hình Delft-3D đối với vùng cửa cũng như cho các vùng khác thuộc dải ven biển Việt Nam. xây dựng công trình thủy lợi 19 Phân tích động lực đập vật liệu địa phương T.S Nguyễn Quang Hùng - Trường Đại học Thủy Lợi 1. Mở đầu: Quá trình vận hành công trình luôn chịu tác động của tải trọng động, đặc biệt là tác động của các vụ nổ mìn hoặc là động đất. Khi chịu tác động của tải trọng động, công trình xây dựng nói chung và công trìnhthủy lợi nói riêng phản ứng như thế nào? Tải trọng động có ảnh hưởng như thế nào đến vấn đề ổn định công trình? Để có lời giải đáp, bài báo này dựa trên lý thuyết cơ bản cố kết Biot động lực, giới thiệu phương pháp tính toán động lực. Thông qua các chuẩn tắc phán đoán hóa lỏng công trình để tính toán mức độ ổn định công trình dưới tác dụng của tải trọng động, Thông qua đó có thể hiểu rõ hơn cơ chế tác động của tải trọng động cũng như phản ứng của công trình trong toàn bộ quá trình xảy ra dao động. Kết quả nghiên cứu góp phần hoàn thiện thiết kế công trình thủy lợi cũng như có ý nghĩa quan trọng trong việc thiết kế kháng chấn. Việc giải các ẩn số: chuyển vị, tốc độ, gia tốc cũng như ứng suất, áp lực lỗ rỗng … tại 1 điểm nào đó trong công trình (hoặc đất nền) được gọi là quá trình phân tích động lực. Thông qua phân tích động lực, một mặt có tác dụng giúp cho quá trình phân tích chấn động trong môi trường đất và hiểu rõ được quy luật truyền tải trọng động trong môi trường. Mặt khác chính xác hóa quá trình thiết kế kháng chấn, đối với thiết kế kháng chấn, tính toán ổn định công trình đất có nhiều lợi ích. Trong quá trình phân tích động lực công trình đất cần thiết phải hiểu rõ được đặc tính động lực của đất, loại hình tải trọng, phân bố chất điểm, điều kiện biên phản ứng như thế nào dưới tác động của tải trọng động. Dựa vào đó để có thể đơn giản hóa bài toán một cách hợp lý, đưa ra được mô hình tính thích hợp, lựa chọn được phương pháp phân tích ứng suất hiệu quả dưới tác dụng của tải trọng động một cách chính xác. 2. Phương trình cố kết động lực cơ bản Biot [1,2,6] Phương trình cân bằng động lực: 0, ==+ ugiiij &ρρσ (1) Phương trình liên tục thấm 0/, =Γ++− piiij &&& ωε (2) Nguyên lý ứng suất hiệu quả Teraghi: ijijij pδσσ += ' (3) Quan hệ Biến dạng ~ Chuyển vị: )( 2 1 ,, ijjiij uu +−=ε (4) Phương trình cân bằng động lực áp lực lỗ rỗng: )( , * ififiijj ugpk && ρρω +−−= (5) Biến đổi ba phương trình (1), (2), (4) dưới dạng số gia: Phương trình cân bằng: jjijiiij ugp &∆=∆+∆+∆ ρρδσ ,,' (6) Phương trình liên tục thấm: 0/, * , * ,, * , =Γ∆+∆−∆+∆−∆ pukgkpku jifijjifijijiijii &&& ρρδ (7) Quan hệ ứng suất- Biến dạng : klijklij D εσ ∆=∆ ' (8) xây dựng công trình thủy lợi 20 Trong đó: σij: ứng suất tổng cộng. σ’ij: ứng suất hiệu quả εij: Biến dạng. ωi: chuyển vị tương đối của nước trong lỗ rỗng so với cốt đất. p: áp lực lỗ rỗng. ρ: mật độ của đất ρs: mật độ của hạt đất ρf: mật độ của nước trong lỗ rỗng của đất 3. Phương pháp giải phương trình động lực Thông thường khi giải phương trình động lực thường dùng phương pháp phần tử hữu hạn (FEM). Quá trình giải thường có 2 giai đoạn: giai đoạn 1 là tiến hành rời rạc hóa miền tính toán và phương trình động lực, giai đoạn 2 là dùng phương pháp phù hợp tiến hành giải bài toán trong miền biến thời gian hoặc biến dao động. Khi phân tích bài toán động (tải trọng tác dụng là tải trọng dao động), ngoài việc xét tác dụng của tải trọng như: trọng lượng bản thân, áp lực nước cũng như các tải trọng ngoài khác còn phải xét tới lực quán tính cũng như lực cản (dao động tắt dần) của đất. Dưới tác dụng của lực cản của đất, năng lượng sóng bị phát tán mà tiêu hao. Đối với công trình nói chung và công trình thủy lợi nói riêng, thường phân lực cản này thành lực dính và lực ma sát cản. Các học giả trước đây thường cho rằng tốc độ biến dạng thay đổi là do tần số chuyển động của miền tính toán. Hiện nay, cũng với miền tính toán như vậy đã chứng minh được rằng tốc độ biến dạng thay đổi có liên quan đến mức độ tiêu hao năng lượng nội bộ trong miền tính toán mà không có liên quan gì đến tốc độ chuyển động. Dựa vào mức độ lớn nhỏ của chuyển vị và biến dạng, dùng mô hình đàn hồi dẻo để nghiên cứu tính chất dính của đất đá và phân tích ứng suất hiệu quả của nó[1,2]. Do sự lan truyền của sóng trong đất đá quyết định sự hao tán năng lượng mà sự lan truyền sóng này có liên quan đến tần suất của sóng cho nên để giải quyết và tính toán phân tích tải trọng động người ta thường chủ yếu nghiên cứu lực cản dính của đất đá[1,5]. 3.1. Phương trình không gian rời rạc miền tính toán. Khi lựa chọn phương pháp phần tử hữu hạn để giải quyết bài toán tất cả lực, điều kiện cân bằng lực phương trình cân bằng động lực đều viết cho điểm nút. Dùng chuyển vị chất lỏng tổng cộng Ui (có cùng tham số với ui) để thay thế chuyển dịch Darcy i: n uU iii ω+= (mô hình cơ học đất bão hòa động lực u-U do Zienkiewicz & Shinomic đưa ra năm 1984). Lựa chọn phương pháp toán học Galerkin đưa ra được phương trình khống chế viết dưới dạng ma trận của phương pháp phần tử hữu hạn: }{}{}]{[}]{[}]{[ FPdKvCaM p =+++ (9) Trong đó: {a}: vector gia tốc      = pj pj U u a & & }{ {v}:vector vận tốc      = pj pj U u v & & }{ xây dựng công trình thủy lợi 21 {d}:Vector chuyển vị   = pj pj U u d}{ [M]: ma trận khối lượng tổng thể; [C]: ma trận cản; [Kp]: ma trận cứng tổng thể; {P}: Vector nội lực; {F}: Vector ngoại lực; 3.2. Giải phương trình động lực. Khi giải phương trình động lực (9), thường dùng phương pháp Newmark-β . Tuy nhiên khi dùng phương pháp này vẫn còn hạn chế: không thể đồng thời đạt được cả độ chính xác về thành phần cao tần của tổn hao năng lượng và độ chính xác của nghiệm giải phương trình. Hiber và đồng nghiệp (1977) đã đưa ra phương pháp Hilber-Hughes-Taylor- để giải quyết vấn đề này. Phương pháp này dùng 3 tham số α, β, γ, với điều kiện: 0 3 1 ≤≤− α , ( ) ( )21 4 1 ,21 αβαγ −=−= 2 1 , đối với bài toán tuyến tính, không những ổn định vô điều kiện mà đối với thành phần cao tần của lực cản còn đạt được độ chính xác bậc 2, đối với thành phần thấp tần thì những hiện tượng dị thường tương đối ít hơn so với các phương pháp khác. Trong miền thời gian, phương pháp Newmark- β xác định vector chuyển vị {d}và vector vận tốc {v}như sau: Với biến thiên thời gian t của bước thời gian thứ n: {d}n+1={d}n+{v}nt+({a}n+1+(1/2-){a}n)t2 (10) {v}n+1={v}n+({a}n+1+(1-){a}n) t (11) Dựa vào phương pháp Hilber-Hughes-Taylor-α phương trình (9) được biến đổi thành: [M]{a}n+1+(1-α )[C]{v}n+1-α[C]{v}n + (1+α)[Kp]{d}n+1–α[Kp]{d}n+(1+α){P}n+1-α{P}n =(1+α){F}n+1-α{F}n (12) Đối với mỗi bước thời gian đều tiến hành tính toán theo các bước dưới đây để tìm ra được số gia gia tốc {a}i+1. Từ đó có thể tiến hành tính toán được số gia chuyển vị , biến dạng, ứng suất cũng như áp lực lỗ rỗng tại bước thời gian tính toán: [M]ieff{a}i+1=Ψ(i) (13) Trong đó: i: bước tính lặp thứ i trong bước tính toán n. n i nnp i np i n i nnn i PPdKdKvCaMFF }{}){1(}]{[}]{)[1(}]{)[1(}]{[}{}){1( )( 1 )( 1 )( 1 )( 11 )( αααααααψ ++−++−−−−−+= +++++ (14) )(22)( ][)1(][)1(][)1(][][ iTp i eff KtKtCtMM ∆++∆++∆++= βαβαγα (15) Trong đó: [KT] Ma trận tiếp tuyến cứng tổng thể. Từ phương trình (13) giải ra được số gia gia tốc {a}i+1. Tiến hành điều chỉnh lại chuyển vị, vận tốc và gia tốc tại bước lặp tiếp theo: xây dựng công trình thủy lợi 22 2)1( 1 1 1 )1( 1 1 1 )1()( 1 1 1 }{}{}{ }{}{}{ }{}{}{ tadd tavv aaa i n i n i n i n ii n i n ∆+= ∆+= ∆+= + + + + + + + + + + + + β γ (16) Dùng kết quả 11 1 1 1 1 }{;}{;}{ + + + + + + i n i n i n dva tính toán được từ phương trình (16) thay vào phương trình (14) tính được giá trị Ψ (i+1) tại bước tiếp theo: n i nnp i np i n i nnn i PPdKdKvCaMFF }{}){1(}]{[}]{)[1(}]{)[1(}]{[}{}){1( )1( 1 )1( 1 )1( 1 )1( 11 )1( αααααααψ ++−++−−−−−+= ++++++++++ (17) Điều kiện hội tụ: 2)1( 1 )1( 1)0( )1( }{ }{ ε εψ ψ ≤∆ ≤ + + + + i n i i a a (17) Tiến hành lặp từ phương trình(13)-(16) đến khi nào thỏa mãn điều kiện hội tụ (17) đồng nghĩa với việc giải xong bài toán tại bước thời gian n. 4. Hóa lỏng công trình. Đất bão hòa chịu tác động của tải trọng cắt lặp đi lặp lại động đất theo phương ngang hoặc là tải trọng chấn động lặp đi lặp lại theo phương đứng dẫn tới có sự thay đổi sắp xếp của các hạt đất, lỗ rỗng giảm nhỏ, đất bị ép chặt lại, một phần ứng suất của các hạt đất truyền sang nước trong lỗ rỗng, áp lực nước trong lỗ rỗng vượt qua áp lực nước tĩnh dẫn tới ứng suất hiệu quả bị giảm nhỏ. Khi áp lực lỗ rỗng đạt đến giá trị trên ứng suất hiệu quả của đất cát, cường độ kháng cắt động lực hoàn toàn bị triệt tiêu sẽ xảy ra hiện tượng hóa lỏng. Dưới tác dụng của tải trọng chấn động kéo dài, quá trình phát triển của ứng suất là hai quá trình hoàn toàn khác nhau nhưng giai đoạn phát triển liên quan đến nhau. Một là do ứng suất từ cốt đất dần dần truyền qua cho nước trong lỗ rỗng và hai là do nước trong lỗ rỗng truyền qua cho cốt đất. Giai đoạn phát triển trước là quá trình chấn động hóa lỏng, giai đoạn sau là chấn động nén ép. Nghiên cứu đặc trưng động lực của đẩt có nghĩa là ngăn chặn, phòng tránh nguy hại của hóa lỏng do chấn động đối với công trình. Lợi dụng chấn động nén ép (giai đoạn sau của chấn động) để giải quyết vấn đề phòng tránh hóa lỏng và gia cố (xử lý) nền công trình. Đất bão hòa phát sinh hóa lỏng luôn thỏa mãn hai điều kiện. Một là tác động của chấn động đạt được đến mức đủ phá vỡ kết cấu (tác động của tải trọng động đủ lớn hoặc là cường độ kết cấu của đất cát tương đối yếu). Hai là sau khi kết cấu phát sinh phá hoại, hạt đất phát sinh chuyển dịch mà xu thế phát sinh là nén ép chứ không phải là trương nở. Ví dụ: cát chặt chịu tác động của chấn động, một mặt, kết cấu của nó rất dễ bị phá hoại nhưng mặt khác kết cấu của nó không dễ dàng bị chấn động phá hoại hoàn toàn. Dưới tác động của chấn động, các hạt cát có xu thế là trương nở chứ không phải là nén ép, không thỏa mãn điều kiện hóa lỏng cơ bản. Do vậy nên cát chặt không dễ phát sinh hóa lỏng. Chỉ có tải trọng động trong tình trạng vượt qua gia tốc giới hạn (kéo dài trong một thời gian nhất định hoặc là vô hạn) hoặc là vượt quá giới hạn số lần chấn động mới có thể xuất hiện tình trạng áp lực lỗ rỗng và biến dạng gia tăng đột biến. Khi gia tốc của tải trọng động đạt đến gia tốc phá hoại hoặc xây dựng công trình thủy lợi 23 là số lần của tải trọng tuần hoàn đạt đến số lần chấn động phá hoại, áp lực lỗ rỗng tăng đến giá trị lớn nhất, ứng suất tăng đến giá trị lớn nhất sẽ xảy ra hiện tượng hóa lỏng. Để phán đoán đất nền hoặc là đập vật liệu địa phương đã phát sinh hóa lỏng hay chưa, một mặt cần phải tính toán ứng suất động ( dzxdyzdxydzdydx ,,,,,, ,,,,, τττσσσ ) của đất nền hoặc là thân đập dưới tác dụng của tải trọng động. Mặt khác phải xác định khả năng kháng hóa lỏng (cường độ kháng cắt tuần hoàn : Nτ ) của đất nền hoặc thân đập dưới tác dụng của tải trọng động. Cường độ kháng cắt vòng Nτ cho biết mẫu cát chịu chấn động đạt đến sự phá hoại hóa lỏng hoặc là ứng suất cắt trên mặt phá hoại của giá trị biến dạng nào đó. Công thức tính cường độ chống cắt: ctgud ′+′−= φστ )( (18) Trong đó: du : áp lực lỗ rỗng động. φ ′ : Góc nội ma sát trong động. c′ : Lực dính đơn vị động. Đối với chuẩn tắc phá hoại hóa lỏng do tải trọng động, năm 1961 Hoàng văn Hi đưa ra “tỉ lệ áp lực lỗ rỗng chấn động 3σ du ” làm chỉ tiêu xác định độ hóa lỏng trong thí nghiệm chấn động ba chiều[3]. Năm 1966, Seed khi làm thí nghiệm cắt tuần hoàn, lấy biến dạng dọc trục đạt đến 10% để làm tiêu chuẩn phá hoại hóa lỏng[4]. Năm 1978, Phân hội động lực học hiệp hội công trình đất đưa ra định nghĩa hóa lỏng như sau : “Quá trình chuyển đổi từ trạng thái ổn định sang trạng thái hóa lỏng của vật thể bất kể là do tĩnh lực, động lực, sóng biển, mạch xung (sự thay đổi lên xuống của dòng nước) hoặc là sự thay đổi của áp lực nước dưới đất đều gọi là hóa lỏng ”. Sự hóa lỏng của đất cát là do sự gia tăng của áp lực lỗ rỗng và sự suy giảm của ứng suất hiệu quả tạo thành. Ba tiêu chuẩn phá hoại hóa lỏng: (1) Tiêu chuẩn 1 Đối với thí nghiệm chấn động ba trục : 1 3 =σ du :xảy ra hiện tượng hóa lỏng Đối với thí nghiệm cắt chấn động đơn : 1 0 =σ du : xảy ra hiện tượng hóa lỏng (19) Trong đó: ud : áp lực lỗ rỗng động σ3 : áp lực cố kết xung quanh của thí nghiệm. σ0 : áp lực cố kết trục đứng của thí nghiệm. (2) Tiêu chuẩn 2 Trong thí nghiệm cắt tuần hoàn, khi biến dạng dọc trục (thí nghiệm chấn động 3 trục) hoặc biến dạng cắt (thí nghiệm cắt 1 trục) đạt đến 5% hoặc 10%, có thể cho là hóa lỏng. (20) Dựa vào những định nghĩa trên đây về hóa lỏng, ở đây không thể gọi là hóa lỏng mà đúng hơn phải gọi là phá hoại. Thí nghiệm đã chứng minh, khi cK =1/3 ≤1.5, hai tiêu chuẩn trên đây về cơ bản là tương đồng. Trước khi xảy ra chấn động, khi hiệu số ứng suất chính trong đất tương đối lớn (các ứng suất chính khác biệt tương đối lớn) thì áp lực lỗ rỗng không đạt được 100% độ hóa lỏng. Tại đường biên của đập vật liệu địa phương, mặt nghiêng của đập (mái đập) cũng thuộc vào tình trạng này. hình thức phá hoại hóa lỏng của đất cát bão hòa dưới mặt đất là đất cát sủi nước (nước bị ép ra xây dựng công trình thủy lợi 24 ngoài) nhưng mái dốc hoặc nền đập không phải là hình thức này. Tuy nhiên mặt nền hoặc là than đập chỉ cần phát sinh 30~40% độ hóa lỏng, mái đập hoặc thân đập có thể phát sinh trượt ở quy mô lớn. Sở dĩ như vậy là vì khi đào mái đất nghiêng hoặc là thiết kế mái đập, nếu không xét đến áp lực lỗ rỗng chấn động, hệ số ổn định trong khoảng 1.3~1.4 thì khi phát sinh 30%-40% độ hóa lỏng, mái đập sẽ mất ổn định. (3) Tiêu chuẩn 3 Lấy áp lực lỗ rỗng giới hạn của cân bằng cực hạn làm tiêu chuẩn, khi du đạt đến giá trị cru coi là hóa lỏng. Giá trị cru được tính như sau: Trong đó: dϕ′ là góc ma sát động, các kí hiệu khác như đã giải thích ở trên. Trong các tiêu chuẩn để trình bày ở trên, cường độ chống cắt Nτ tính được từ tiêu chuẩn 3 nhỏ hơn rất nhiều so với tiêu chuẩn 1 và tiêu chuẩn 2. Có thể nói rằng: cùng với 1 giá trị Nτ , dùng tiêu chuẩn 3 thu được số lần chấn động đạt đến mức độ phá hoại là nhỏ nhất. Do vậy, vì tiêu chẩn hóa lỏng không thồng nhất, việc các định chính xác Nτ là hết sức không hợp lý. 5. Tính toán ứng dụng 5.1. Sơ đồ tính 5.2. Chỉ tiêu cơ lý Bảng 1. Chỉ tiêu cơ lý dùng trong tính toán Thứ tự Loại đất Ex=Ey Mpa µx=µy γ KN/m3 kx=ky10-7 m/s 1 I 30 0.25 19.1 3 2 II 30 0.25 19.1 4.5 3 III 30 0.25 21.0 1 4 IV 30 0.25 21.0 1 5.3. Lưới phần tử tính toán Hình 1. Sơ đồ tính toán hóa lỏng công trình dưới tác dụng của tải trọng động xây dựng công trình thủy lợi 25 Hình 2. Sơ đồ lưới phần tử tính toán Lướ i phần tử t ính toán bao gồm: 1546 nút đỉnh, 723 phần tử . Phổ động đấ t được dùng trong tính toán là trận động đấ t Loma Prieta vớ i gia tốc động đấ t lớn nhấ t là: 1.2g trong thờ i gian 6s (xảy ra tạ i San Francisco năm 1989)[9 ] . 5.4. Kết quả tính toán Hình 3. Tỉ lệ áp lực lỗ rỗng tại một số phần tử Hình 4. Vùng hóa lỏng công trình tại thời điểm t= 6s 5. Kết luận Bài báo này cơ bản dựa trên lý thuyết cố kết động Biot, sau đó dùng phương pháp tích phân từng bước (số gia từng bước) để giải phương trình động lực. So sánh với với phương pháp tăng chồng chấn động, phương pháp phân tích phổ, phương pháp này có những ưu điểm sau: Không chỉ giải quuyết được bài toán tuyến tính mà còn có thể giải quyết được bài toán phi tuyến. Có thể không nhất thiết phải giải tần suất chấn động bản thân của công trình mà trực tiếp giải ra được chuyển vị và ứng suất tại các thời điểm bất kì. Tuy nhiên tính toán chịu hạn chế bởi tính chính xác và mức độ ổn định, bước thời gian tính toán không được phép quá lớn. Tuy nhiên hạn chế này có thể được khắc phục bởi việc chọn lựa 1 cách thích đáng các hình thức tham số sai phân đưa thêm vào để cái tiến. Ngoài ra phương pháp giải này còn có thể dùng để giải bài toán phân tích ứng suất hiệu quả. Do vậy, trước mắt phương pháp tích phân từng bước là một phương pháp chủ yếu để giải bài toán động lực và đặc biệt là bài toán tải trọng động trong thời gian hết sức ngắn (ví dụ như tải trọng động do nổ mìn, động đất…) . Hơn nữa, khi dùng phương pháp Hilber-Hughes-Taylor-α có thể khắc phục được hạn chế đã nêu ở trên: đồng thời có thể đạt được độ chính xác về tổn hao năng lượng cao tần và độ chính xác của nghiệm giải phương trình, đối với thành phần cao tần của lực cản còn đạt được độ chính xác bậc 2, đối với thành phần thấp tần thì những hiên tượng dị thường đương đối ít hơn so với các phương xây dựng công trình thủy lợi 26 pháp khác. Không những thế, thông qua phân tích công trình đất dưới tác dụng của tải trọng động đã góp phần làm sang tỏ được những phản ứng của công trình dưới tác dụng của động đất – là một hiện tượng tự nhiên thường xảy ra các nguy hại cho công trình. Dưới tác dụng của tải trọng động, thông qua 3 tiêu chuần phá hoại hóa lỏng, có thể tính đoán phán đoán sự phá hoại của công trình hay đất nền. Từ đó có thể hoàn thiện hơn quá trình thiết kế ổn định, thiết kế kháng chấn công trình nói chung và công trình đất nói riêng. Tài liệu tham khảo [1 Steven L. Kramer. Geotechnical Earthquake Engineering. University of Washington, 1996 [2] Cung Hiểu Nam. Tính toán phân tích công trình đất. Nhà xuất bản kiến trúc1999. (bản tiếng Trung). [3] Hoàng Văn Hi, Nghiên cứu hóa lỏng nền cát và mái dốc cát. Kỹ thuật công trình thủy lợi thủy điện.1959 – kỳ 15. [4] Seed H B and Lee K L. Liquefaction of Saturated Sands During Cyclic Loading. Proc. ASCE. J. SMFD. 1966. 92(SM6): 105~134. [5] R.B.J.Brinkgereve Plaxis B.V., Netherlands and P.A. Vermeer. University of Stuttgart .Finite Element code for Soil and Rock Analyses. 1998 [6] Biot M.A.General theory of Three Dimensional Consolidation. J.App.Physis, Vol. 12 ,155. 1941 [7] Geudehus, G., Finite Elements in Geomechanics, Jonhn Wiley&Sons, 1977 [8] Nguyễn Quang Hùng, Phúc Thiếu Quân, Trần Sinh Hồng, dùng phương pháp FEM nghiên cứu bước thời gian tự thích ứng trong quá trình cố kết. Báo Cơ học đất Trung Quốc, 2005-3. (bản tiếng Trung) [9] Seed.R.B. Dickenson.S.E. Reimer.M.F. Bray.J.D. Sitar.N. Mitchell.J.K. Idriss.I.M. Kayen.R.E. Kropp.A. Harder. L.F and Power.M.S. Preliminary report on the principal geotechnical aspects of the October 17,1989 Loma Prieta earthquake. Report UCB/EERC-90/05, Earthquake Enginerring Research Centre, University of California, Becrkeley. 1990. Chống bùn cát bồi lắng hồ Ba Bể PGS.TS. Lưu Như Phú I. Đặc điểm hồ Ba Bể: Hồ Ba Bể là hồ tự nhiên lớn nhất nước ta, là một trong 20 hồ tự nhiên lớn nhất trên thế giới. Hồ nằm giữa Vườn Quốc gia Ba Bể, một vùng rừng nguyên sinh, có độ đa dạng sinh học cao. Vườn có diện tích 7611 ha, nếu kể cả khu đệm có diện tích 32.200ha. Hồ Ba Bể có diện tích gần 600ha, dài gần 9000m, rộng 200 đến 1200m, chứa hơn 90 triệu m3 nước (Hình 1). Mặt hồ cao hơn mực nước biển 150m. Có 4 sông đổ vào hồ Ba Bể là sông Chợ Lèng, Bó Lù và Tà Han có diện tích lưu vực F = 420km2, chiều dài dòng chính từ 19÷27km, tổng lưu lượng lũ là Q2 = 840m3/s. 3 sông này đổ trực tiếp vào hồ về phía Nam. Sông Năng đổ vào hồ ở phía Bắc. Tính đến cửa hồ, sông Năng có F = 1420km2; Qlũ= 980m3/s, Qkiệt = 4,2m3/s. Do có dung tích lớn, hồ Ba Bể có khả năng điều tiết hơn 40 triệu m3 nước cho sông Gâm, điều hoà lưu lượng và mực nước cho một vùng rộng lớn. Giảm bớt tác hại của lũ dồn về của 4 con sông với lưu vực gần 2000km2. Mặt khác, hồ còn bổ sung nước mùa kiệt vào sông Năng phục vụ cho sản xuất và đời sống nhân dân với lưu lượng hàng ngày Q=4÷10m3/s. Hồ Ba Bể còn là lá phổi của Vườn Quốc gia, điều hoà khí hậu cho một vùng rộng lớn, giữ nhiệt độ ôn hoà trung bình năm là 22oC. II. Bồi lắng bùn cát và hiểm hoạ bồi lấp hồ Ba Bể: xây dựng công trình thủy lợi 27 Đã từ nhiều năm nay, từ 4 hướng, phù sa bồi lắng thành nhiều bãi lớn đang lấn lấp lòng hồ. Năm 2002 lượng bồi từng cửa sông như sau: Cửa sông Chợ Lèng: 184.000 m3 Cửa sông Bó Lù: 111.000 m3 Cửa sông Tà Han: 97.000 m3 Cửa Hồ: 30.000 m3 Lượng bồi 3 sông phía Nam chiếm 93% tổng lượng, với diện tích bãi nổi 144 ha bãi lấn vào hồ 10÷60m. III. Các yếu tố tác động đến bồi lắng: Sự bồi lấp bùn cát ở các cửa sông là hậu quả của dòng chảy bùn cát trên mặt lưu vực do mưa rào rơi xuống và dòng chảy bùn cát trong sông. Hai vấn đề này có liên quan mật thiết với nhau, nhưng lại có quy luật vận động khác nhau: Quy luật xói mòn mặt lưu vực và quy luật vận chuyển bùn cát trong sông. Có 2 loại nhân tố tác động đến quá trình xói, bồi gồm: + Yếu tố thiên nhiên: Mưa, địa hình, thổ nhưỡng, thảm thực vật, hình dạng và độ dốc sông suối .v.v. + Yếu tố con người: Trồng rừng, phá rừng, khai hoang, canh tác, công trình các loại. Yếu tố về rừng, lớp phủ thực vật đóng vai trò quan trọng. Nước mưa rơi xuống, một phần ngấm vào đất thành dòng chảy ngầm, chỉ còn một phần chảy trên mặt đất. Hệ số dòng chảy khu có rừng nhỏ hơn nhiều so với khu trống, trọc. Một thực nghiệm ở độ dốc địa hình cấp 5 (i≥35°) cho thấy: - Rừng tự nhiên: Hệ số dòng chảy α = 4,5% - Cây bụi: Hệ số dòng chảy α = 16.0% - Lúa nương: Hệ số dòng chảy α = 31.0% - Đất trống, trọc: Hệ số dòng chảy α =50÷80% Các vùng đất trống, trọc, mưa rơi trực diện làm tan rã, kéo các hạt đất, dần tạo thành những rãnh sâu chằng chịt, đó là những dòng bùn cát đầu nguồn của sông suối. Địa hình lưu vực Ba Bể phức tạp, chủ yếu là đồi núi bị chia cắt (Bảng 1); 90% diện tích này thuộc vùng núi ở độ cao 300÷1500m. Phần lớn địa hình rất dốc. Đa phần ở độ dốc i ≥16°. Hầu hết đất sườn núi thuộc loại không chặt, độ rỗng n = 31÷37%. Vì thế khi mưa xuống, đất dễ bị trương nở, dễ bị bào mòn, xói lở. xây dựng công trình thủy lợi 12.09 217.9F (ha) 0 14.8 112.7125 120 0 0 60 10 180 0 120 20 6 3w(x10 m ) 0 300 50 240 40 3.070.33 420 70 360 60 155 150 145 140 1350 130 Z (m) §iÓm lò xuÊt hiÖn n¨m 1986 120Z (m) 135130125 W ~ z 273.7 343 481.5 614 640.88 w (x 10 m ) 61.3640.0824.69 540 90 480 80 96.189.69 F (ha)600 100 36 150145140 156155 F ~ Z Hình 1. Biểu đồ diện tích mặt nước và dung tích hồ Ba Bể Bảng 1. Đặc điểm địa hình lưu vực T T Kiểu địa hình Độ cao (m) Diện tích (km2) % Cấp độ dốc (0) % Khả năng xói lở 1 Núi cao >1000 92,3 16,8 V >35 48.1 Dễ xói lở 2 Núi đá 700÷800 115,2 21,0 IV 26÷35 28,7 Dễ xói lở 3 Núi thấp 300÷700 283,2 51,6 III 16÷25 15 Dễ bào xói 4 Đồi <300 20,3 3,7 II 8÷15 1,7 Dễ rửa trôi 5 Hồ và thung lũng 150÷300 38,0 6,9 I <8 6,5 Hạn chế xói mòn Tổng cộng 549,0 100 100 Trước năm 1970 toàn bộ khu vực quanh hồ Ba Bể là rừng nguyên sinh. Gần đây, rừng bị chặt phá mạnh. Phân tích ảnh vệ tinh, viễn thám loại LANDSAT, SPOT đa phổ màu tỷ lệ 1/50.000 và 1/100.000 ở 2 thời điểm 1990 và 2000 thấy được sự biến đổi lớp phủ thực vật lưu vực này (Bảng 2). Bảng 2. Sự biến đổi đất, rừng trên các lưu vực sông đổ vào hồ Ba Bể (km2) Lưu vực sông Chợ Lèn Bó Lù Tà Han Thời điểm Đất trống + lúa màu Cây bụi, cỏ Rừng Đất trống + lúa màu Cây bụi, cỏ Rừng Đất trống + lúa màu Cây bụi, cỏ Rừng F 13,829 42,85 9 137,41 4 8,931 25,048 102,16 2 2,938 1,944 83,18 5 199 0 %/ 1970 7,1 22,1 70,8 6,6 18,4 75,0 3,3 3,3 93,4 xây dựng công trình thủy lợi 29 F 35,821 73,78 6 83,784 17,998 58,992 59,609 15,600 10,63 1 63,33 6 200 0 % /1970 18,5 38,2 43,3 13,2 43,2 43,6 17,4 11,9 70,7 Vùng dự án gồm 14 xã, 98% là đồng bào dân tộc sống chủ yếu về nông lâm nghiệp. Các mảnh ruộng nương phân tán nhỏ lẻ, canh tác tự nhiên một vụ nhờ nước trời, nước chảy tràn bờ, đất bị rửa trôi. Hầu hết các công trình thủy lợi là công trình tạm, mùa lũ bị nước cuốn trôi ở khu vực này. Đời sống nhân dân không ổn định nên tăng cường chặt phá rừng, đốt nương làm rẫy, tạo nhiều vùng đất trống, trọc tạo điều kiện cho mưa lũ bào xói lưu vực mạnh mẽ hơn. Tốc độ phá rừng đã thống kê được như sau: + Trong 20 năm (1970 ÷1990), rừng đã bị chặt phá mất 23% diện tích lưu vực, tốc độ trung bình 1,15% năm. + Trong 10 năm tiếp theo (1990 ÷2000), rừng tiếp tục bị chặt phá mất 27,7%, tốc độ trung bình 2,77%/năm. Tức là tốc độ phá rừng gấp 2,4 lần so với giai đoạn trước. IV. Diễn biến bồi lấp hồ Ba Bể: Dựa vào những tài liệu điều tra, khảo sát đo đạc các yếu tố liên quan, đã phân tích và thiết lập 2 mô hình toán để tính xói mòn trên toàn lưu vực các sông đổ vào hồ Ba Bể. Trong đó đã sử dụng những tài liệu thủy văn và bùn cát thực đo trong 2 năm để chỉnh sửa mô hình. + Tính xói mòn lưu vực hồ Ba Bể, sử dụng mô hình SEAGIS dựa trên nguyên lý phương trình mất đất phổ dụng USLE với sự trợ giúp của hệ thống thông tin địa lý. Chia lưu vực ra 30 ô, tính lượng xói mòn từng ô rồi tính lượng bùn cát chuyển tải ở các ô về các mặt cắt nhập lưu của các suối nhánh với sông. Tổng lượng bùn cát xói mòn trên lưu vực vào các năm 1970, 1990, 2001, 2002 tương ứng là 11,5; 29,07; 67,42; 86,88 vạn tấn/năm. Từ 1970÷2001, rừng bị chặt phá mất 50,7% diện tích lưu vực, lượng đất bị xói mòn tăng 58,5 lần. + Tính quá trình dòng lũ của mạng các sông với hồ và dòng bùn cát chuyển vận từ các cửa nhập lưu về bồi lấp các cửa sông bằng mô hình NAM+ MIKE11 + GIS. Kết quả tính lượng bùn cát xói mòn vận chuyển về lấn lấp hồ năm 2002 ghi trên Bảng 3. Bảng 3. Lượng bùn cát xói mòn lưu vực và bồi lắng các cửa sông năm 2002 Lượng bùn cát Lượng xói mòn Lượng bồi cửa sông tính toán Lượng bồi thực đo Sai số Tỷ lệ bồi tính toán Cửa sông (vạn tấn) (vạn tấn) (vạn tấn) % % Chợ Lèng 42,464 33.,97 31,229 8,8 80,0 Bó Lù 24,039 20,.2 18,802 7,4 84,0 Tà Han 20,373 17,.9 16,490 8,6 87,9 Tổng cộng 86,87 72,.07 66,520 8,3 83,0 Ta thấy 4/5 lượng đất cát bị bào xói trên lưu vực vào mùa mưa đã bị dòng lũ cuốn về các cửa sông, lấn lấp hồ. Vì vậy phải có những giải pháp trên sông tác động trực tiếp vào dòng lũ bùn cát mới triệt giảm hiểm hoạ cho hồ. Kết quả tính tổng lượng bồi hàng năm, trong 28 năm (1975÷2002) ghi trên Hình 2. Sự bồi lấp bùn cát ở 3 cửa sông đổ vào hồ Ba Bể có thể chia ra 2 thời kỳ: + 1975÷1989: Lượng bồi tăng từ từ, trung bình 1.280 tấn/năm. + 1990÷2002: Lượng bồi tăng nhanh; 3.4000 tấn/năm (gấp 2,7 lần giai đoạn trước). Nguyên nhân xây dựng công trình thủy lợi 30 chính là do tốc độ phá rừng giai đoạn này gấp 2,4 lần giai đoạn trước, dẫn đến xói mòn bồi lấp mạnh mẽ. Dựa vào tốc độ lấn lấp hồ hiện nay (đo và tính toán) có thể tính ra khoảng 90 năm nữa hồ Ba Bể sẽ bị lấp đầy. V. Giải pháp chống bồi lấp hồ Ba Bể: Có 2 nhóm giải pháp chủ yếu để chống bồi lấp hồ Ba Bể: 1. Giải pháp lâm nghiệp, xã hội: + Bảo vệ và trồng rừng: Cho đến năm 2020 phải đạt độ che phủ trở lại như năm 1990 (chiếm 77%), phải trồng được 11.000 ha (600 ha/năm), dự định đầu tư 45 tỷ đồng. Đây là một kế hoạch phải đầy khó khăn và lâu dài. + Đầu tư cơ sở vật chất, ổn định và nâng cao đời sống nhân dân, điều kiện tiên quyết để gìn giữ và chống chặt phá rừng. 2. Giải pháp thủy lợi: Lưu vực hồ Ba Bể rất rộng, gần 500km2, thuộc miền núi cao, sông suối rất dốc, lũ lớn. Trên sông suối không có các hồ đập kiên cố để triệt giảm tốc độ dòng lũ, chặn giữ bùn cát; mùa mưa lũ cuốn 4/5 lượng đất bị bào xói trên lưu vực về bồi lấp hồ. Giải pháp chủ động chống bồi lấp hồ có hiệu quả trong một số năm tới là đầu tư xây dựng các công trình thủy lợi trên các sông đổ vào hồ Ba Bể. Dự án đã tiến hành khảo sát, tính toán, thiết kế hệ thống các công trình thủy lợi trên 3 sông này với mục tiêu đạt được lợi ích tổng hợp: - Chống úng ngập cho 2200 ha ruộng, đưa lên cấy 2 vụ ăn chắc. - Cấp nước, sinh hoạt cho 8.500 dân. - Triệt giảm lượng bùn cát lấn lấp hồ, đảm bảo hồ tồn tại lâu dài. - Cải tạo môi trường sinh thái và cảnh quan toàn vùng. 3. Giải pháp công trình: + Trên sông Chợ Lèng: xây các hồ đập kiên cố: Đập Bản Chán; đập Hin Khao; hồ đập Bản Phẹc; đê bãi Pắc Ngòi. + Trên sông Bó Lù và Tà Han: xây hồ Đồng Lạc (14 triệu m3) trên sông Bó Lù xây hồ Xuân Lạc (25 triệu m3) trên sông Tà Han, nối thông 2 hồ này bằng 1 đường hầm dài gần 3 km. Nước 2 hồ này được xả qua một kênh vào sông Năng, sau Thác Đầu Đẳng chảy về sông Gâm. Lũ, bùn cát chủ yếu của 2 sông này xả ra ngoài lưu vực Ba Bể, không cho chảy vào hồ. xây dựng công trình thủy lợi 31 0 10 20 30 40 50 60 70 80 1975 1980 1985 1990 1995 2000 2005N¨m T æn g l− în g (v ¹n tÊ n/ n¨ m ) Tæng l−îng S«ng Chî LÌng S«ng Bã Lï S«ng Tµ Han Hình 2. Tổng lượng vận tải bùn cát các sông đổ vào hồ Ba Bể + Đầu tư cơ sở vật chất, ổn định và nâng cao đời sống nhân dân, điều kiện tiên quyết để gìn giữ và chống chặt phá rừng. 2. Giải pháp thủy lợi: Lưu vực hồ Ba Bể rất rộng, gần 500km2, thuộc miền núi cao, sông suối rất dốc, lũ lớn. Trên sông suối không có các hồ đập kiên cố để triệt giảm tốc độ dòng lũ, chặn giữ bùn cát; mùa mưa lũ cuốn 4/5 lượng đất bị bào xói trên lưu vực về bồi lấp hồ. Giải pháp chủ động chống bồi lấp hồ có hiệu quả trong một số năm tới là đầu tư xây dựng các công trình thủy lợi trên các sông đổ vào hồ Ba Bể. Dự án đã tiến hành khảo sát, tính toán, thiết kế hệ thống các công trình thủy lợi trên 3 sông này với mục tiêu đạt được lợi ích tổng hợp: - Chống úng ngập cho 2200 ha ruộng, đưa lên cấy 2 vụ ăn chắc. - Cấp nước, sinh hoạt cho 8.500 dân. - Triệt giảm lượng bùn cát lấn lấp hồ, đảm bảo hồ tồn tại lâu dài. - Cải tạo môi trường sinh thái và cảnh quan toàn vùng. 3. Giải pháp công trình: + Trên sông Chợ Lèng: xây các hồ đập kiên cố: Đập Bản Chán; đập Hin Khao; hồ đập Bản Phẹc; đê bãi Pắc Ngòi. + Trên sông Bó Lù và Tà Han: xây hồ Đồng Lạc (14 triệu m3) trên sông Bó Lù xây hồ Xuân Lạc (25 triệu m3) trên sông Tà Han, nối thông 2 hồ này bằng 1 đường hầm dài gần 3 km. Nước 2 hồ này được xả qua một kênh vào sông Năng, sau Thác Đầu Đẳng chảy về sông Gâm. Lũ, bùn cát chủ yếu của 2 sông này xả ra ngoài lưu vực Ba Bể, không cho chảy vào hồ. Giới thiệu một số kết quả khảo sát nền và mang cống dưới đê bằng phương pháp Rađa KS. Ngô Trí Côi - KS. Đỗ Anh Chung - KS. Nguyễn Văn Lợi I. Đặt vấn đề : Việt Nam có hàng ngàn cống dưới đê phục vụ cho mục đích tưới, tiêu và lấy sa. Chỉ tính riêng xây dựng công trình thủy lợi 32 cho hệ thống đê sông Hồng và sông Thái Bình trong phạm vi đề tài nghiên cứu đã có hơn 1.000 cống. Công việc xây dựng và vận hành đã trải qua nhiều giai đoạn lịch sử khác nhau nên hiện tượng hư hỏng cũng xảy ra khác nhau. Các hư hỏng có thể phân loại theo các nhóm: hỏng cửa van; hỏng khớp nối; lún, nứt, chuyển vị tường cánh thượng, hạ lưu; lún, nứt, chuyển tường thân, trên cống; bồi, xói lở, bể tiêu năng phía đồng, phía sông; thấm qua nền, mang cống; sạt mái đá thượng, hạ lưu… Có nhiều phương pháp khảo sát đánh giá chất lượng nền và mang cống như: phương pháp khoan lấy mẫu, phương pháp Rađa, điện, đánh dấu đồng vị phóng xạ… cho kết quả khả quan. Trong khuôn khổ bài giới thiệu này, các tác giả giới hạn ở phương pháp Rađa địa và thuỷ thám. II. Địa điểm và kết quả nghiên cứu: 1. Địa điểm nghiên cứu: Địa điểm được chọn để thử nghiệm khảo sát nền và mang cống là cống Nhâm Lang thuộc huyện Hưng Hà-Thái Bình. Cống được thiết kế xây dựng từ năm 1936 và 1982, sau khi xây dựng bổ sung phía đồng, cứ sau 5-7 năm lấy sa lại xẩy ra hiện tượng sụt đất ở 2 mang cống. Gần đây nhất là vào năm 2000 đã xảy ra sụt đất ở mang cống vai trái phía đồng tạo thành hố sụt hình phễu có đường kính 2,0m và sâu gần 4,0m. Đơn vị quản lý phải xử lý lấp hết hơn 4,0m3 đất. Điều này chứng tỏ đất phía dưới tiếp giáp mang cống đã bị rửa trôi. Nguyên nhân nào dẫn đến hiện tượng này, hiện tại vẫn chưa lý giải được. 2. Kết quả nghiên cứu bằng Rađa địa thám: Nghiên cứu sóng điện từ nửa không gian trong lòng đất cho thấy khả năng thu nhận tín hiệu phản xạ trên bề mặt phụ thuộc vào tính chất điện của đối tượng gây phản xạ và môi trường truyền sóng. Tín hiệu thu nhận được tuân thủ theo biểu thức: 21 21 εε εε + −=R Trong đó: R: hệ số phản xạ ε1: hằng số điện môi của môi trường. ε2: hằng số điện môi của đối tượng gây phản xạ. Sự sai biệt hằng số điện môi càng lớn thì hệ số phản xạ càng cao, theo đó sự phân dị trên giản đồ sóng càng rõ. Ngoài ra, kích thước của đối tượng nghiên cứu được phân dị trên giản đồ phụ thuộc vào tần số trung tâm của ăngten ứng dụng. Tần số ăngten càng cao thì độ phân dị đối tượng có kích thước bé và độ sâu nghiên cứu giảm trong cùng môi trường và ngược lại. Độ phân dị đối tượng phụ thuộc vào biểu thức: d = 1/4 λ Trong biểu thức: d: kích thước của đối tượng nghiên cứu λ: bước sóng Từ công thức trích dẫn cho thấy kích thước của đối tượng nghiên cứu tối thiểu phải bằng 1/4 bước sóng. Kích thước thiết bị càng lớn hơn 1/4 bước sóng thì độ phân dị càng rõ trên giản đồ sóng Rađa. Theo nguyên lý đã nêu trên, nhóm tác giả đã ứng dụng các loại hình ăngten có tần số trung tâm là 400MHz, 200Mhz, 100MHz và 80MHz khảo sát trên cống Bùng-Nghệ An, cống Vạn Phúc- Bắc Ninh và cống Nhâm Lang-Thái Bình. Hệ thiết bị sử dụng là SIR-10B do hãng GSSI-Mỹ chế tạo. Xử lý số liệu trên chương trình phần mềm RADAN for WINDOW bằng module không gian 3 xây dựng công trình thủy lợi 33 chiều (3D) do hãng cung cấp thiết bị lập chương trình cho kết quả chính xác vị trí, kích thước của đối tượng nghiên cứu (Hình 1). Hình 1- Lát cắt 3D Rađa địa thám trên mang cống trái Nhâm Lang-Thái Bình Trên giản đồ sóng Rađa xác định đối tượng tại vị trí cách mang cống 4,0m có kích thước 3,5÷4m ở độ sâu 1,5m. Đây là vị trí hố sụt xảy ra năm 2000 và đã được xử lý lấp bịt. Những thông tin này giúp cho cơ quan quản lý là Công ty khai thác công trình thuỷ lợi Bắc Thái Bình đề ra các biện pháp xử lý thích hợp. 3. Kết quả nghiên cứu bằng Rađa thuỷ thám: Việc ứng dụng Rađa thuỷ thám là dựa vào nguyên lý cơ bản của Rađa địa thám. Sự khác biệt ở đây là cải tiến ăngten địa thám, mở rộng mạch xung để vượt qua lớp nước đồng nhất nhằm khảo sát các đối tượng nằm sâu dưới lớp nước. Phương pháp được ứng dụng tại cống Nhâm Lang-Thái Bình. Đây là cống hở nên việc thao tác tại hiện trường có nhiều thuận lợi. Trên giản đồ sóng Rađa (Hình 2) quan sát từ phía đồng ra phía sông cho thấy: - Nền cống có xu hướng lún về phía đồng với biên độ dao động từ 0,3÷0,5m. Đặc biệt tại vị trí tiếp giáp nền cống. - Dưới lớp bê tông tại tiếp giáp giữa tường ngực và tường cánh có dị thường sâu 1,0÷1,5m, chiều rộng 2,5÷3,0m. - Bể tiêu năng phía đồng bị lắng bồi với chiều dài 8÷10m. Các hiện tượng bất thường ở nền cống cần được làm sáng tỏ bằng biện pháp khoan lấy mẫu thí nghiệm các chỉ tiêu cơ lý lực học để thiết kế xử lý. Tại tiếp giáp tường ngực và tường cánh bố trí tuyến cắt ngang. Trên giản đồ sóng Rađa (Hình 3) cho thấy dưới lớp bê tông có một số dị thường hình hypecbol và dị thường phân lớp. Các dị thường tiếp giáp với tường trái có dạng hang rỗng (từ 1,5÷4,5m). Khoảng giữa của tuyến đoạn từ 7,5÷15m dị thưòng phân lớp không đồng đều. Đoạn cuối tuyến từ 17,5÷25m có dạng phân lớp. xây dựng công trình thủy lợi 34 MÆt Hình 2. Giản đồ Rađa thuỷ thám cắt dọc nền cống từ phía đồng ra phía sông Hình 3. Giản đồ Rađa thuỷ thám cắt ngang nền cống MÆt n-íc III. Kết luận: Kết quả nghiên cứu nền và mang cống bằng Rađa địa thám và Rađa thuỷ thám giải quyết các vấn đề: - Xác định vị trí và kích thước dị vật phần đất tiếp giáp mang cống. - Xác định dị vật dưới lớp bê tông nền cống. - Xác định mức độ bồi lắng bể tiêu năng. Cá c phương pháp Tính toán cọc chịu tải trọng ngang và Momen - Kiến nghị một phương pháp thực hành thiết kế KS. Trần Văn Thái - KS. Phan Đình Tuấn Trong thiết kế công trình thuỷ lợi thường gặp bài toán tính toán cọc chịu tải trọng ngang và moment. Các phương pháp tính toán hiện nay có rất nhiều, bài báo này nêu tổng quan một số phương pháp tính và đề xuất một phương pháp thực hành đơn giản, thuận tiện. I. Mở đầu Một đặc điểm quan trọng đối với móng cọc trong công trình thuỷ lợi là cọc vừa chịu tải trọng đứng vừa chịu tải trọng ngang và mô men do áp lực nước gây ra. Dưới tác dụng của lực ngang cọc có 2 khả năng bị phá hoại: - Đối với cọc ngắn chế tạo bằng vật liệu cứng thì khả năng chịu lực thường bị hạn chế bởi biến xây dựng công trình thủy lợi dạng ngang quá lớn, khi đó đất ở mặt bên bị phá hoại. - Đối với cọc dài, cọc thường bị phá huỷ do ứng suất quá lớn trước khi đất ở mặt bên bị phá hoại. Người thiết kế phải kiểm tra an toàn của cọc theo 3 điều kiện: 1) N < [ N]; 2) M < [M]; 3) [ ]yy ∆<∆ . Trong đó: N, M - lực thẳng đứng và mô men tác dụng lên đầu cọc. [N], [M] - khả năng chịu tải trọng đứng và mô men của cọc. [∆y] - chuyển vị đầu cọc cho phép. Quá trình kiểm tra theo điều kiện thứ nhất thường là đơn giản. Vấn đề khó khăn nhất đối với người thiết kế là kiểm tra theo điều kiện 2, 3. Dưới đây là tổng quan các phương pháp tính toán M, ∆y đối với cọc chịu tải ngang và mô men, trên cơ sở đó đề xuất một phương pháp tính toán nhanh và hợp lý. II. Tổng quan các phương pháp tính toán sức kháng bên của cọc thẳng đứng [1] Hình 1. Sức kháng bên với cọc đứng có đầu tự do chịu tải trọng ngang. 1. Phương pháp Phương pháp n ngắn và cứng. Phương pháp n đặt tải trọng và cho qvx KPux .= δ Trong đó: −vxδ áp lực qu C - Lực dính củ Kc, Kq - hệ số; l e Q M = Q. e §iÓm uèn xrdx Q35 Brinch Hansen (1961) [1] ày dựa trên lý thuyết về áp lực đất đối với đất có C và ϕ, chỉ áp dụng đối với cọc ày xác định tâm quay bằng cách lấy mô men của tất cả các lực đối với một điểm = 0. Sức kháng cực hạn của đất ở độ sâu nào đó: cKC.+ á tải hiệu quả thẳng đứng; a đất; Kc, Kq =    B xf ;ϕ . B p ux xây dựng công trình thủy lợi 36 Đối với các tải trọng ngắn hạn, lấy độ bền thấm không thoát nước Cu và ϕ (= 0). Đối với tải trọng quan trọng và dài hạn, phải dùng giá trị độ bền hiệu quả thoát nước (c’, ϕ’). 2. Phương pháp Brom (1964) [1] Đã đề xuất phương pháp xác định sức kháng bền đối với cọc thẳng đứng dựa trên lý thuyết áp lực đất nhưng đã thừa nhận một số giả thiết đơn giản: - Đất sử dụng xem đơn thuần là đất không dính (C = 0) hoặc đơn thuần là đất dính (ϕ = 0). - Chỉ xem xét riêng từng loại cọc: cọc ngắn và cứng, cọc dài và mềm. - Cọc ngắn đầu tự do bị quay quanh một tâm quay với những cọc có đầu cố định lại bị dịch chuyển ngang (hình 2a,b). l T©m quay Qu §Çu tù do l Qu a, Đầu tự do b, Đầu cố định Hình 2: Chuyển vị quay và dịch chuyển ngang với cọc ngắn chịu tác dụng ngang. Dạng biến dạng của cọc dài khác với cọc ngắn vì cọc dài sẽ không bị xoay và dịch chuyển do sức kháng bị động của đất rất lớn ở phần dưới cọc (hình 3 a, b). Do vậy, khả năng chịu tải trọng bên của 2 loại cọc này được xác định bằng các phương pháp khác nhau. Hình 3: a, Đầu tự do b, Đầu cố định 3. Cọc dài trong đất dính phương pháp Brichhasen Trong đất dính, phản lực nền của đất giảm ngoài khoảng (1,5 B ∼ xo). Với cọc dài thì sự phân bố của phản lực đất theo chiều sâu (hình 4a, b). xây dựng công trình thủy lợi 37 γ. Β. xo.K γ. Β. xo.K Hình 4: a, Đầu tự do b, Đầu cố định III. Phương pháp chuyển vị cho phép khi chịu tải trọng bên làm việc [1] Khi thiết kế cọc chịu tải trọng bên thường dựa trên chuyển vị ngang cho phép hơn là dựa trên khả năng chịu tải trọng ngang cực hạn. Chính vì vậy trong bài báo này, tác giả chỉ đi sâu vào nghiên cứu phương pháp này, đó là phương pháp xấp xỉ phản lực nền (Reese và Matlock 1956 và 1960). 1. Tính toán cọc chịu tải trong ngang theo TCXD 205 - 1998 [3] Đất quanh cọc được xem như môi trường đàn hồi biến dạng tuyến tính đặc trưng bằng hệ số nền CZ (T/m3) (hình 5a, b). Khi không có những số liệu thí nghiệm, CZ xác định theo công thức. CZ = K. Z Trong đó: K - Hệ số tỉ lệ (T/m4) (lấy từ bảng G.1 - TCXD 205-1998). Z - Độ sâu của vị trí tiết diện cọc, kể từ mặt đất với đài cao, kể từ đáy đài với cọc đài thấp. Tất cả tính toán được thực hiện theo chiều sâu tính đổi của vị trí tiết diện cọc trong đất Le, xác định theo công thức. Ze = αbd. Z; Le = αbd . L ; JE bK b cS bd . .=α Trong đó: Z và L - Chiều sâu thực tế vị trí tiết diện cọc trong đất và chiều sâu hạ cọc thực tế (mũi cọc) (m);J - Mô đun quán tính tiết diện ngang dọc (m4); K - Hệ số tỉ lệ; Eb - Mô đun đàn hồi ban đầu của bê tông khi nén và kéo (T/m2), lấy theo tiêu chuẩn thiết kế cấu bê tông cốt thép; bC - Chiều rộng qui ước của cọc (m) * Tính toán chuyển vị ngang của cọc ở mức đáy đài và góc xoay theo công thức sau: JE Ml JE Hly b O b O OOn 23 23 +++=∆ ϕ ; JE Ml JE Hl b O b O O ++= 2 2 ϕϕ lO - Chiều dài đoạn cọc (m). (bằng khoảng cách từ đáy đài cọc đến mặt đất). xây dựng công trình thủy lợi 38 yO và ϕo - chuyển vị ngang (m), và góc xoay của tiết diện ngang của cọc (radian), xác định theo công thức: yO = HO . δHH + MO . δHM. ϕO = HO . δMH + MO .δHM. ∆n ϕ ϕo H M N δMH δHH H0 =1 δMMM0 =1 Hình 5: a- Sơ đồ tải trọng b- Sơ đồ chuyển vị của cọc trong đất trên cọc Trong đó: HO - Giá trị tính toán của lực cắt (T), lấy HO = H. MO - Mô men uốn (Tm), lấy MO = M + HlO. δHH - Chuyển vị ngang của tiết diện (m / T), bởi HO = 1. δHM - Chuyển vị ngang của tiết diện (1 / T), bởi MO = 1. δmM - Góc xoay của tiết diện (1 / T), bởi HO = 1. δMM - Góc xoay của tiết diện (1 / T), bởi MO = 1. Chuyển vị δHH , δMH δHM và δMM xác định theo công thức : δHH = JEbbd . 1 3α . AO. ; δMH = δHM = JEbbd . 1 2α . BO; δMM = JEbbd . 1 α . CO. (1 - 16) Trong đó: AO, BO, CO - Những hệ số không thứ nguyên. Từ đây có thể xác định được: áp lực ngang tính toán lên đất σZ (T/m2), lực cắt QZ (T), MZ (Tm):     +−= 13121 .DJE H C JE AyZK bbd O bbd O Oc bd Z αα ϕ ασ ( )3 bd O 3O3Obbd3Ob 2 bdZ D H CMB.J.EA.y.JEM α++ϕα−α= xây dựng công trình thủy lợi 39 QZ = αbd3 Eb .J.yO .A4 - αbd2. Eb J.ϕO .B4 + + αbd MO. C4 + HO.D4 ; NZ = N. Trong đó: K = Hệ số tỉ lệ xác định bảng G.1; A1, B1, C1, D1; A3, B3, C3, D3; A4, B4, C4, D4: là các hệ số Nhận xét. Trong tiêu chuẩn TCXD 205-1998 cho khối lượng tính toán quá cồng kềnh và phức tạp nhưng chỉ xác định được momen và chuyển vị đầu cọc, còn để tính biểu đồ phân bố momen trong thân cọc thì hầu như khó thực hiện được . Như vậy, từ những phân tích trên trong báo cáo này tác giả kiến nghị dùng phương pháp xấp xỉ phản lực nền của tác giả (Reese và Matlock - 1956) để tính toán. 2. Phương pháp xấp xỉ phản lực nền (Reese và Matlock - 1956): kiến nghị làm phương pháp thực hành. 2.1. Nội dung của phương pháp. Trong phương pháp này thừa nhận rằng đất tác dụng giống như một dãy các lò xo đàn hồi độc lập tuyến tính. Q M P y x x y Q M P Hình 6: a. Cọc chịu tải trọng ngang trong đất b. Cọc chịu tải trọng ngang tựa trên các lò xo Theo mô hình này đất đàn hồi được thay thế bằng một loạt các lò xo đàn hồi độc lập nhau, có độ cứng Kh thể hiện hình 6a, b và công thức: y pK y = Palmer và Thompson (1948) đã sử dụng công thức sau để tính mô đun phản lực nền. Kx = Kh    L x 4 Với:Kh - Giá trị của Kx tại x = L hoặc ở đỉnh cọc. x - Độ sâu một điểm nào đó dọc theo cọc. n ≥ 0. Phản lực của đất Cọc trước khi chịu tải trọng P = K. yx x xây dựng công trình thủy lợi 40 2.2 Nhận xét: * Trong phương pháp này, việc xác định nội lực trong cọc khi chịu tác dụng của sức chịu tải trọng ngang chủ yếu phụ thuộc vào giá trị của mô đun phản lực nền Kh của lò xo. Việc xác định giá trị Kh cho hợp lý là một yêu cầu quan trọng và nhiều nhà cơ học đất đã kiến nghị, đề ra những công thức để tính Kh như : Tezaghi (1955), Davisson (1960). Palmer và Thompson (1948)… - Với đất dạng hạt, Tezaghi kiến nghị Kh tỉ lệ với độ sâu.Kh = nh. x. - Với đất sét chất tải bình thường, Davisson chỉ ra rằng Kh tỉ lệ với độ sâu. - Với đất bồi lắng cũng được lý luận tương tự như vậy với Peck và Davisson (1962). * Những điểm còn tồn lại của (Reese và Matlock): đó là độ chính xác của Kh và nh không có giới hạn, nh luôn cần thiết biến đổi để tạo ra sự biến đổi mômen trong. Do vậy việc xem xét hợp lí các hiện tượng để dẫn tới kết luận Kh chủ yếu là trị số kinh nghiệm trong tự nhiên và thay đổi theo một số yếu tố như: Chuyển vị, độ sâu, đường kính cọc, loại gia tải, tốc độ chất tải, số tải trọng tác dụng. 2.3. Kiến nghị của tác giả. Qua phân tích và tìm hiểu, tác giả đã đi sâu vào kiến nghị dùng công thức [2] sau để xác định mođun phản lực nền dùng để tính cọc: KS = S1. AS + S2. BS. Zn Trong đó: S1 = S2 - 1 - với cọc vuông; S1 = 1,3 ~ 1,7 - với cọc tròn; S2 = 2 ~ 4,4 - với cọc tròn; n - hệ số (n = 0,4); AS. BS - là các hệ số được biểu diễn; AS = 40 ( C. NC. SC + 0,5. γ. B. Nγ. S γ); BS = 40 (γ. Nq. Sq) Với: SC , Sγ , Sq - các hệ số; Nq , NC , Nγ = f (ϕ , C, γ); γ - Trọng lượng riêng của đất. - Kh = KS. B ; B - Bề rộng của cọc. - Sau khi xác định giá trị KS bằng cách lập bảng tính. Sử dụng phần mềm Sap 2000 đưa thông số KS vào để tính các giá trị M , N , Q. Ví dụ tính toán: Vẽ biểu đồ momen dọc thân cọc 35x35x1600, biết rằng cọc có đầu không xoay được và chịu tải trọng ngang H=5T.m. Sau đây là tính chất của lớp đất dọc thân cọc: Lớp đất STT Cao độ đầu lớp Cao độ đáy lớp Chiều dày lớp Lực dính c (kN/m2) (độ) (kN/m 3) Nc Nq Ng 1 0 4 4 10 4 14,3 6,98 1,48 0,40 2 4 8 4 12 6 14,1 7,76 1,82 0,64 3 8 12 4 14 8 14,4 8,68 2,26 0,92 4 12 16 4 22 8 14,2 8,68 2,26 0,92 Modun phản lực nền. c gSố TT Độ sâu Nc Nq Ng (kN/m2) (kN/m3) As Bs ks(Kn/m3) KS (KN/m) 0 0,00 1 6,98 1,48 0,4 10 14,3 2832 847 2832 991 xây dựng công trình thủy lợi 41 1 1,00 1 6,98 1,48 0,4 10 14,3 2832 847 3679 1288 2 2,00 1 6,98 1,48 0,4 10 14,3 2832 847 4525 1584 3 3,00 1 6,98 1,48 0,4 10 14,3 2832 847 5372 1880 4 4,00 1 6,98 1,48 0,4 10 14,3 2832 847 6218 2176 5 5,00 2 7,76 1,82 0,64 12 14,1 3788 1026 8920 3122 6 6,00 2 7,76 1,82 0,64 12 14,1 3788 1026 9947 3481 7 7,00 2 7,76 1,82 0,64 12 14,1 3788 1026 10973 3841 8 8,00 2 7,76 1,82 0,64 12 14,1 3788 1026 12000 4200 9 9,00 3 8,68 2,26 0,92 14 14,4 4954 1302 16669 5834 10 10,00 3 8,68 2,26 0,92 14 14,4 4954 1302 17971 6290 11 11,00 3 8,68 2,26 0,92 14 14,4 4954 1302 19273 6746 12 12,00 3 8,68 2,26 0,92 14 14,4 4954 1302 20575 7201 13 13,00 4 8,68 2,26 0,92 22 14,2 7730 1284 24418 8546 14 14,00 4 8,68 2,26 0,92 22 14,2 7730 1284 25701 8995 15 15,00 4 8,68 2,26 0,92 22 14,2 7730 1284 26985 9445 16 16,00 4 8,68 2,26 0,92 22 14,2 7730 1284 28269 9894 Biểu đồ momen trong thân cọc 35x35 x1600cm, chịu tải trọng ngang H=5T(thể hiện trong hình 7). xây dựng công trình thủy lợi 42 Hình 7. Kết quả: Momen max đầu cọc là 7.8T.m; Chuyển vị ngang đầu cọc là 1cm. IV. Kết luận. Qua nghiên cứu một số phương pháp xác định ″sức chịu tải ngang của cọc”. Tác giả chỉ đưa ra một số nhận xét trên cơ sở so sánh phương pháp tính toán giá trị lực ngang tiêu chuẩn tính theo tiêu chuẩn Việt Nam 205-1998 ứng với chuyển vị ngang cho phép và phương pháp xấp xỉ phản lực nền là: - Tính toán theo quy phạm thường thiên về an toàn, xét về mặt kinh tế là không có lợi (lãng phí cọc). - Trị số lực ngang tiêu chuẩn tính theo chuyển vị ngang cho phép còn phụ thuộc vào yêu cầu của thiết kế: Nếu ∆n ≤ 1cm thì tính toán theo tiêu chuẩn. Nếu ∆n > 1cm thì phải xác định theo thí ngiệm tải trọng tĩnh (khối lượng công việc là tương đối lớn và mất nhiều thời gian ). - Tính theo tiêu chuẩn 205-1998 sẽ gặp khó khăn khi xét đến sự phức tạp của địa chất (như nền nhiều lớp..v.v.) - Phương pháp xấp xỉ phản lực nền cho khả năng: + Tính toán nhanh, tiết kiệm thời gian cho người thiết kế nhờ sử dụng phần mềm SAP2000. + Cho ra được hình dạng của biểu đồ Momen để từ đó có thể bố trí cốt thép cho hợp lí.❐ xây dựng công trình thủy lợi 43 ứng dụng Mô Hình Delft-3D... (Tiếp theo trang 16) Tài liệu tham khảo 1. Nguyễn Thị Phương Thảo (2004) Msc Thesis- Estuarine Morphology, case study in Hai Phong, Viet Nam, UNESCO – IHE Institute for water education – Delft- The Netherlands. 2. Nguyễn Thị Thảo Hương (2000) Luận án tiến sỹ khoa học, Nghiên cứu diễn biến cửa sông thuỷ triều phục vụ cho khai thác thuỷ lợi, giao thông, Trường ĐH Thuỷ Lợi. 3. De Vriend et al (1993), approach to long-term modeling of coastal morphology. Elsevier Science publishers B.V, Amsterdam. 4. WL | Delft Hydraulics (2003), Delft3D User Manual, Version 3.10, March 2003.