Ikhtisar teknik
Bendungan multi-lengkungan Foziling, bendungan multi-lengkungan beton bertulang pertama setelah berdirinya Republik Rakyat Tiongkok, dibangun pada tahun 1952, terletak di Sungai Bi di Provinsi Anhui. Ini adalah waduk pengatur tahunan dengan manfaat pemanfaatan komprehensif untuk pengendalian banjir, irigasi, pembangkit listrik, dan pasokan air. Total kapasitas tampungan adalah 1,94 × 10 8 m 3 , dan tinggi tampungan normal waduk, tinggi banjir rencana, dan tinggi banjir periksa masing-masing adalah 125,56 m, 125,97 m, dan 129,80 m. Gambar 1 menunjukkan pemandangan bendungan bagian atas, hulu, dan hilir. Puncak bendungan memiliki panjang 510 m dan lebar 1,8 m. Seperti yang ditunjukkan pada Gambar 2, bendungan multi-lengkungan Foziling terdiri dari 20 penopang, 21 lengkungan, dan bagian bendungan gravitasi. Ketinggian puncak bendungan adalah 131,06 m dan tinggi maksimum bendungan adalah 75,9 m, termasuk tinggi dinding gelombang 1,5 m.
Lengkungan dan penopang bendungan multi-lengkungan merupakan struktur berdinding tipis dan ringan. Penopang dibentuk dengan menghubungkan dua dinding vertikal berbentuk segitiga dan dinding partisi, seperti ditunjukkan pada Gambar 3. Dinding penopang tipis di bagian atas tetapi tebal di bagian bawah, dengan ketebalan 0,6 m pada puncak dan 1,97 m pada ketinggian 70,0 m dari puncak. Bagian lengkung berbentuk setengah lingkaran dengan jari-jari bagian dalam 6,75 m, sedangkan jari-jari bagian luar bervariasi menurut ketebalan lengkungan. Dari puncak bendungan ke bawah, bagian lengkungan berangsur-angsur menebal, dengan ketebalan 0,6 m di puncak bendungan dan 2,0 m pada ketinggian 70,0 m dari puncak. Karena kapasitas penyimpanan yang kecil, variasi ketinggian air yang besar, medan yang kompleks dan kondisi geologi pondasi bendungan, serta pengendalian suhu yang tidak memadai selama konstruksi, bendungan telah menjalani tiga kali pemeliharaan dan perkuatan skala besar sejak dibangun. layanan 35 .
Penguatan pertama dilakukan pada bulan Juni 1965 sampai dengan Oktober 1968. Untuk mengatasi permasalahan kavitasi dan getaran pipa baja pada penopang dan kebocoran retakan, maka reservoir dikosongkan, dan lengkungan yang diperkuat dituangkan di bagian bawah No. –16 lengkungan, pipa baja diperkuat pada penopang, dan sambungan susut pada setiap penopang dipasang. Selain itu, dilakukan grouting konsolidasi dan perkuatan tirai pada pondasi bendungan, serta penambahan fasilitas pemantauan bendungan. Oleh karena itu, standar pengendalian banjir waduk ditingkatkan.
Penguatan kedua dilakukan pada bulan Juni 1982 sampai dengan Juni 1986. Akibat kecelakaan luapan banjir pada bulan Juli 1969, dimana muka air waduk tertinggi melebihi puncak dinding gelombang sebesar 1,08 m, maka bendungan ditambah 1,5 m dan satu saluran pelimpah tambahan dibangun. Sementara itu, pondasi bendungan dan tirai pondasi tepi kanan bendungan diperkuat untuk meningkatkan keamanan bendungan.
Penguatan ketiga dilakukan pada bulan Oktober 2002 hingga April 2005. Setelah musim banjir tahun 1993, waduk beroperasi pada permukaan air yang lebih tinggi dalam waktu yang lama. Arus dingin yang tiba-tiba terjadi pada akhir November, menyebabkan perpindahan setiap penopang melebihi nilai ekstrim historisnya. Meskipun bendungan telah diperkuat dua kali, retakan yang parah masih terjadi. Berdasarkan hasil penilaian keselamatan tahun 2002, keamanan waduk dinilai Kelas III, termasuk dalam waduk berbahaya dan memerlukan penguatan kembali. Isi perkuatan utama dapat dibagi menjadi empat bagian: (1) Proyek pengolahan pondasi, termasuk pembersihan dan penggalian terak, penimbunan beton ekspansi mikro pada penopang, grouting konsolidasi tekanan tinggi untuk pondasi bendungan, grouting tirai untuk pondasi bendungan hulu, dilakukan dan lubang drainase pondasi serta lubang pengamatan rembesan ditambahkan pada kedua tepian. (2) Proyek perkuatan, termasuk perbaikan retakan badan bendungan, perkuatan ketebalan untuk daerah bertekanan tinggi di bagian hulu penopang bendungan dengan beton bertulang serat baja (SFRC) yang dicor di tempat, partisi beton bertulang horizontal baru di dalam penopang, dan penyemprotan SFRC permukaan hulu dan hilir lengkungan No.2 dan 22. (3) Lapisan insulasi termal poliuretan disemprotkan di zona fluktuasi ketinggian air hulu (ketinggian hulu di atas 110 m dari dinding) dan permukaan bendungan hilir (lihat Gambar. perkuatan ketebalan untuk daerah tegangan tinggi di bagian hulu penopang bendungan dengan beton bertulang serat baja (SFRC) yang dicor di tempat, partisi beton bertulang horizontal baru di dalam penopang, dan penyemprotan SFRC pada permukaan hulu dan hilir No.2 dan 22 lengkungan. (3) Lapisan insulasi termal poliuretan disemprotkan di zona fluktuasi ketinggian air hulu (ketinggian hulu di atas 110 m dari dinding) dan permukaan bendungan hilir (lihat Gambar. perkuatan ketebalan untuk daerah tegangan tinggi di bagian hulu penopang bendungan dengan beton bertulang serat baja (SFRC) yang dicor di tempat, partisi beton bertulang horizontal baru di dalam penopang, dan penyemprotan SFRC pada permukaan hulu dan hilir No.2 dan 22 lengkungan. (3) Lapisan insulasi termal poliuretan disemprotkan di zona fluktuasi ketinggian air hulu (ketinggian hulu di atas 110 m dari dinding) dan permukaan bendungan hilir (lihat Gambar. 3 ). (4) Sistem pemantauan otomasi keamanan bendungan dan fasilitas pengelolaan lainnya ditambahkan.
Instrumen untuk pemantauan bendungan
Fluktuasi ketinggian air, perpindahan bendungan, dan variasi suhu termasuk beton bendungan, air di depan bendungan, dan udara sekitar, terus dipantau selama pengoperasian waduk.
Dua termometer (TQ13_1, TQ13_2) yang memantau suhu udara di dekatnya masing-masing dipasang di bagian depan dan ekor rongga yang dibentuk oleh lengkungan No.13 dan penopang yang terhubung dengannya. Termometer lain (T) dipasang di kotak kisi-kisi di puncak dinding gelombang bendungan gravitasi di tepi kanan untuk memantau suhu udara di dekatnya. Di sepanjang sisi timur panel hulu penopang No.11 dipasang 14 termometer (TS11085-TS11124) untuk memantau suhu air pada ketinggian berbeda. Dua termometer (T_wall1, T_wall2) dipasang untuk memonitor suhu di dalam dan di luar lapisan insulasi masing-masing pada dinding barat penopang No.21, begitu pula untuk lapisan insulasi di sisi timur lengkungan No.21 (T_arch1, T_arch2). Untuk memantau suhu beton dan udara di penopang No.13, 3 . Termometer bernomor T13-1–T13-16 baru ditambahkan pada penguat ini, dan termometer bernomor T13011–T13102 adalah asli.
Seperti ditunjukkan pada Gambar 4 , 20 garis tegak lurus normal (PL2–PL5, PL6–PL12, PL13-1, PL13-2, PL14–PL21) yang digunakan untuk memantau perpindahan puncak bendungan dipasang di No. 2–5, 7– 21 penopang, di antaranya terdapat dua garis tegak lurus normal pada penopang No. 13 (lih. Gambar 3 ). Tiga garis tegak lurus terbalik (IP3, IP13, IP21) yang digunakan untuk memantau deformasi pondasi dipasang pada penopang No.3, 13, dan 21.
Hasil pemantauan
Ketinggian & suhu air
Gambar 5 menunjukkan fluktuasi tinggi muka air setelah dilakukan perkuatan. Untuk irigasi dan pembangkit listrik yang optimal, ketinggian air dioperasikan pada ketinggian antara 105,0 dan 125,0 m.
Data pemantauan suhu udara ditunjukkan pada Gambar 6 . Dapat ditemukan bahwa suhu udara berkorelasi signifikan dengan musim. Di musim dingin, suhu di dekat bagian hilir lengkungan lebih rendah dibandingkan di sisi hulu, namun di musim panas justru sebaliknya. Dari data termometer di dalam dan di luar lapisan isolasi yang ditunjukkan pada Gambar 7 , terlihat bahwa variasi tahunan suhu di dalam lapisan isolasi lebih kecil daripada di luar, dan respon waktu suhu di lapisan isolasi adalah sedikit histeris dibandingkan di luar. Fenomena ini menggambarkan bahwa lapisan insulasi mempunyai insulasi termal yang baik.
Gambar 8 menunjukkan suhu air pada ketinggian reservoir yang berbeda selama pengoperasian. Dapat ditemukan bahwa semakin besar variasi suhu tahunan yang terjadi pada ketinggian yang lebih tinggi. Suhu udara dan suhu beton dari titik pengamatan tipikal yang dipilih ditunjukkan pada Gambar 9 . Dapat dilihat bahwa variasi suhu tahunan di dekat permukaan penopang cenderung lebih besar daripada variasi suhu di penopang, dan respons waktu suhu di bagian hulu sedikit lebih histeretik dibandingkan dengan respons waktu di dekat permukaan penopang. Yang perlu diperhatikan adalah beberapa termometer tidak berfungsi, sehingga data dari 57 termometer dipilih sebagai masukan data suhu untuk pembentukan model HTT dan MLSR.
Perpindahan
Dalam makalah ini, perpindahan bendungan mengacu pada perpindahan horizontal penopang bersamaan dengan perpindahan hulu dan hilir (yaitu perpindahan radial) dan perpindahan kiri-kanan penopang sepanjang tepi kiri dan kanan (yaitu perpindahan tangensial). Berdasarkan orientasi spasial penopang, titik pengamatan perpindahan dapat dikelompokkan menjadi tiga bagian: titik pengamatan dasar sungai (PL8–PL16), titik pengamatan tepi kiri (PL2–PL5,) dan titik pengamatan tepi kanan (PL17–PL21). Karena terdapat banyak titik pengukuran perpindahan, titik pengamatan khusus dipilih untuk ditampilkan dan dianalisis.
Gambar 10 menunjukkan perpindahan radial dari penopang No. 3, 13, dan 21 dari waktu ke waktu. Dapat dilihat bahwa: (1) Terdapat hubungan yang signifikan antara perpindahan radial dengan musim. (2) Perpindahan radial penopang di dasar sungai (misalnya PL13_1 dan PL13_2) lebih besar dibandingkan dengan perpindahan di tepi sungai (misalnya PL3 dan PL21). (3) Perpindahan radial penopang di tepi kanan (misalnya PL21) lebih kecil dibandingkan dengan perpindahan di tepi kiri (misalnya PL3). (4) Varians perpindahan radial pada PL3, PL13_1, PL13_2, dan PL21 masing-masing kira-kira 5,0 mm, 3,5 mm, 4,5 mm, dan 2,5 mm.
Gambar 11 menunjukkan perpindahan tangensial penopang terhadap waktu. Hal ini dapat dilihat sebagai berikut: (1) Dengan mengambil dasar sungai penopang (penopang No. 8–16) sebagai titik tengah, maka perpindahan tangensial penopang di tepi kiri (penopang No. 2–6) berkorelasi erat dan positif dengan suhu, sedangkan perpindahannya berlawanan dengan perpindahan tangensial penopang di tepi kanan (penopang No. 17–21). (2) Perpindahan tangensial bersifat siklis, terutama pada penopang di tepian, dan semakin dekat ke tepian, semakin besar perpindahan yang terjadi. (3) Varians perpindahan tangensial pada PL3, PL5, PL12, PL13_1, PL13_2, PL14, PL17, dan PL21 kira-kira 5,0 mm, 2,0 mm, 1,0 mm, 1,5 mm, 1,0 mm, 1,0 mm, 2,0 mm, dan masing-masing 6,0 mm.
Banyaknya instrumen pemantauan memberikan jaminan untuk menguasai perubahan ketinggian air, suhu, dan kondisi kesehatan bendungan setiap saat. Perpindahan adalah efek komprehensif dari tekanan air, suhu dan waktu. Tingginya kemiripan antara perpindahan dan suhu menunjukkan bahwa perpindahan sangat dipengaruhi oleh suhu. Lalu bagaimana faktor individu mempengaruhi perpindahan dan seberapa besar kontribusinya terhadap perpindahan dipelajari pada bagian berikut.
Model statistik untuk pemantauan perilaku bendungan
Model pemantauan hidrostatik, suhu, waktu-perpindahan
Model pemantauan kesehatan bendungan berdasarkan perpindahan merupakan metode yang efektif untuk menganalisis dan memperoleh kondisi kerja bendungan. Perpindahan bendungan terdiri dari dua bagian, yaitu perpindahan reversibel yang disebabkan oleh fluktuasi tinggi muka air waduk dan suhu. Dan bagian lain yang disebut perpindahan ireversibel disebabkan oleh mulur, reaksi agregat alkali dan efek inelastis lainnya 6 , 29 . Oleh karena itu, perpindahan bendungan dapat dituliskan seperti yang dinyatakan di bawah ini:
dimana δ adalah perpindahan terukur; δ H adalah komponen tekanan air; δ T adalah komponen suhu dan δ θ adalah komponen penuaan.
Komponen tekanan air δ H
Komponen tekanan air δ H mewakili perpindahan badan bendungan dan pondasi yang disebabkan oleh tekanan hidrostatik reservoir, dan dihitung sebagai
dimana α i adalah koefisien statistik; n diambil sebagai 3 untuk bendungan gravitasi dan 4 untuk bendungan lengkung; H adalah kedalaman air hulu.
Komponen suhu δ T
Komponen suhu δ T menggambarkan perpindahan panas yang disebabkan oleh variasi suhu badan bendungan dan pondasi. Jadi pilihan terbaik adalah mengambil suhu tubuh bendungan dan pondasi bendungan sebagai faktor perhitungan. Dalam teknik praktis, ada dua situasi berikut mengenai tata letak termometer:
(1) Ketika termometer yang cukup dipasang pada badan dan pondasi bendungan, komponen suhu δ T juga dapat dihitung dengan suhu terukur dalam model HTT
dimana b i , b 1i , b 2i adalah koefisien statistik; m dan m 1 masing-masing adalah jumlah termometer efektif dan lapisan yang dilengkapi termometer; T adalah suhu yang diukur; dan β adalah suhu rata-rata dan gradien masing-masing lapisan.
(2) Bila tidak ada termometer atau sedikit termometer pada badan dan pondasi bendungan, perpindahan panas umumnya diasumsikan mengikuti perubahan musim. Oleh karena itu, model hidrostatis-musim-waktu (HST) banyak digunakan, di mana komponen suhu δ T biasanya disimulasikan dengan fungsi harmonik periode satu tahun atau setengah tahun.
dimana m 2 diambil sebagai 1 atau2; b 1i dan b 2i adalah koefisien statistik; , dan t adalah jumlah hari akumulasi sejak tanggal pemantauan awal.
Bidang suhu bendungan pada dasarnya disebabkan oleh pergantian suhu dengan musim selama periode pengoperasian, sehingga perpindahan termal juga dapat dihitung dengan suhu udara terukur 1 , 4 , 21 . Karena efek konduksi panas, terjadi histeresis temporal pada respon suhu beton terhadap perubahan suhu udara sekitar. Ketika suhu udara yang diukur digunakan untuk mensimulasikan perpindahan termal, komponen suhu δ T dapat dinyatakan sebagai
dimana b dan d adalah koefisien statistik; T adalah suhu udara terukur atau rata-rata; d adalah jumlah faktor suhu. Karena bentuk bendungan multi lengkung yang tipis, pengaruh suhu udara sebulan yang lalu diperhitungkan dalam penelitian ini. Untuk menyederhanakan perhitungan, T d ditetapkan sebagai T 0 , T 1 , T 2 , T 3–5 , T 6–15, dan T 16–30 , di mana T 0 , T 1 , T 2adalah suhu udara terukur pada hari pemantauan, masing-masing 1 dan 2 hari sebelumnya; T 3–5 , T 6–15 dan T 16–30 masing-masing adalah suhu udara rata-rata pada 3 hingga 5, 6 hingga 15, dan 16 hingga 30 hari sebelumnya. Dalam hal ini, d diambil sebagai 6.
Komponen penuaan δ θ
Komponen penuaan δ θ mewakili efek waktu, biasanya dinyatakan sebagai
dimana θ = t/100 , t adalah jumlah hari sejak awal analisis, dan c 1 dan c 2 adalah koefisien.
Regresi bertahap linier berganda
Perlu diperhatikan bahwa fungsi harmonik dalam model HST tidak dapat secara efektif menginterpretasikan perpindahan termal yang disebabkan oleh variasi suhu dinamis jangka pendek. Oleh karena itu, berdasarkan Persamaan. ( 1 ), ( 2 ), ( 3 ), ( 5 ), dan ( 6 ), model pemantauan kesehatan perpindahan HTT multi-arch dapat dinyatakan dalam dua rumusan berikut:
M 1 dan M 2 adalah model pemantauan kesehatan perpindahan HTT berdasarkan pengukuran suhu udara dan suhu tubuh bendungan dan pondasi bendungan . Dimana adalah suhu air pada ketinggian reservoir yang berbeda; a 0 adalah konstanta; δ 0 , H 0 masing-masing adalah perpindahan terukur dan kedalaman air hulu pada tanggal awal dalam rangkaian pemodelan; , , , , ,
T_{w0} T_{0j} T_{0}^{0} T_{1}^{0} T_{2}^{0} T_{3 - 5}^{0} T_{6 - 15}^{0} T_{16 - 30}^{0} \theta_{0} \ln \theta_{0}, , , dan masing-masing adalah suhu dan faktor penuaan tanggal awal dalam rangkaian pemodelan.
Jika kesalahan acak dipertimbangkan, Persamaan. ( 7 ) dan ( 8 ) meluas ke bentuk umum model MLSR, yang dapat ditulis dalam notasi matriks sebagai
dimana adalah vektor variabel terikat ( m × 1) dari perpindahan bendungan dan m adalah jumlah titik pengukuran perpindahan. adalah matriks ( m × ( i + j + 3)) dari variabel bebas yang dinyatakan dalam , , dan . adalah vektor (( n + i + 3) × 1) dari koefisien regresi. adalah ( m
{\mathbf{X}} \delta_{H} \delta_{T} \delta_{\theta } {{\varvec{\upalpha}}} {{\varvec{\upvarepsilon}}}× 1) kesalahan sisa antara perpindahan yang diukur dan perkiraan perpindahan model.
Langkah-langkah penyelesaian model MLSR adalah sebagai berikut: (1) Pemilihan faktor model yang optimal dimulai dari tidak adanya variabel dalam model, pengujian penambahan setiap variabel dengan uji F. (2) Variabel ditambahkan jika penyertaannya memberikan peningkatan kesesuaian yang paling signifikan secara statistik. (3) Proses ini diulangi sampai tidak ada yang memperbaiki model secara signifikan secara statistik. (4) Koefisien regresi Persamaan. ( 9 ) dapat diperoleh dengan metode kuadrat terkecil.
Solusi di atas dapat diwujudkan dengan program MATLAB. Selain itu, komponen tekanan air , komponen suhu , dan komponen penuaan juga dapat diperoleh. Dengan demikian, perilaku bendungan multi-lengkungan Foziling dapat dianalisis lebih lanjut.
\delta_{\theta }
Kriteria kinerja untuk evaluasi kinerja model mencakup koefisien determinasi ( ), rata-rata kesalahan absolut ( ), dan akar rata-rata kesalahan kuadrat ( ).
dimana n adalah jumlah sampel; dan masing-masing adalah perpindahan yang diukur dan dihitung; dan adalah nilai rata-rata perpindahan terukur.
Perilaku analisis perpindahan bendungan multi-lengkungan oleh MLSR
Kumpulan data perpindahan radial, periode 13 tahun (Mei 2006–Juli 2018, 624 observasi) dipilih untuk pemodelan. Observasi bulan Mei 2006 hingga November 2016 dipilih untuk pelatihan model, sedangkan observasi bulan Desember 2016 hingga Juli 2018 dipilih untuk pengujian model. Data perpindahan radial titik pengamatan (PL3, PL13-1, PL21), yang diambil sebagai titik tipikal, dipilih untuk pemodelan dan analisis. Mengingat termometer yang mengukur suhu beton hanya dipasang pada penopang 13, berdasarkan Persamaan. ( 7 ), model pemantauan kesehatan perpindahan HTT ( M 1 ) berdasarkan suhu beton yang diukur ditetapkan. Sedangkan untuk penopang No. 3, 21, model pemantauan kesehatan perpindahan HTT lainnya ( M 2) ditetapkan berdasarkan suhu udara dan suhu air yang diukur. Menurut langkah penyelesaian model MLSR, faktor dan koefisien optimal tercantum pada Tabel 1 . Dan parameter kinerja model pemantauan kesehatan bendungan berdasarkan MLSR tercantum pada Tabel 2 . Perpindahan yang diukur dan dimodelkan ditunjukkan pada Gambar 12 . Semua ini menunjukkan bahwa metode MLSR mempunyai kinerja yang baik dalam penyesuaian dan prediksi.
Komponen tekanan air, komponen temperatur dan komponen penuaan titik pengamatan perpindahan (PL3, PL13-1) dihitung berdasarkan Persamaan. ( 2 )–( 6 ). Selanjutnya, kontribusi ketinggian air, suhu dan waktu terhadap amplitudo perpindahan tahunan juga diperoleh. Gambar 13menunjukkan kontribusi masing-masing komponen terhadap amplitudo tahunan perpindahan radial dari tahun 2007 hingga 2017. Terlihat bahwa: (1) Pada amplitudo tahunan perpindahan radial penopang No.3, komponen tekanan air menyumbang 4,13 –6,66%, komponen suhu menyumbang 91,89–94,37%, dan komponen penuaan menyumbang 1,3–1,5%, dengan nilai rata-rata masing-masing sebesar 93,11%, 5,45% dan 1,44%. (2) Pada amplitudo tahunan perpindahan radial penopang No.13, komponen tekanan air menyumbang 18,50–27,90%, komponen suhu menyumbang 69,58–79,82%, dan komponen penuaan menyumbang 1,21–3,43%, dengan nilai rata-rata masing-masing sebesar 22,58%, 74,65% dan 2,77%. (3) Perpindahan radial penopang di dasar sungai lebih dipengaruhi oleh tekanan air tetapi tidak terlalu dipengaruhi oleh suhu dibandingkan penopang di tepian sungai, dan pengaruh waktu terhadapnya hanya sedikit. (4) Perpindahan radial paling dipengaruhi oleh suhu, diikuti oleh ketinggian air dan waktu.
Gambar 14 menunjukkan hubungan suhu udara dengan komponen suhu yang dipisahkan berdasarkan model HHT di PL13-1. Ini menggambarkan bagaimana suhu mempengaruhi perpindahan radial, yaitu puncak bergeser ke hilir ketika suhu naik, dan puncak bergeser ke hulu ketika suhu turun. Perubahan perpindahan radial tertinggal dari perubahan suhu udara. Di sini, pengaruh suhu merupakan kombinasi perubahan suhu pada lengkungan dan penopang. Namun jika dilihat secara terpisah, pengaruhnya terhadap perpindahan bendungan sangatlah berbeda.
Sedangkan untuk lengkungan, ketika suhu naik, penopang terdorong ke hilir oleh gaya dorong yang disebabkan oleh ekspansi termal dari lengkungan yang terhubung (lih., Gambar 15 a ) , sedangkan suhu turun, penopang ditarik ke atas karena penyusutan. lengkungan (lih. Gambar 15 b). Sedangkan untuk penopang, puncak penopang dipindahkan ke hulu ketika suhu naik (lihat Gambar 15 c), dan puncak penopang dipindahkan ke hilir ketika suhu turun (lihat Gambar 15 d), yang merupakan hal yang sama. untuk bendungan lengkung dan bendungan gravitasi.
Karena proporsi komponen tekanan air relatif kecil, komponen tekanan air yang dipisahkan dari model yang dibuat dengan menggunakan seluruh data menjadi bias. Oleh karena itu, untuk menyoroti pengaruh tekanan air pada perpindahan radial sebanyak mungkin, nilai terukur yang sesuai dengan suhu udara antara 20 dan 23 °C dipilih, dan komponen tekanan air dipisahkan dengan metode MLSR. Gambar 16 mencerminkan hubungan antara komponen tekanan air dan perubahan ketinggian air. Terlihat ketika level reservoir diatas H 0(113,92 m), perpindahan radial berkorelasi positif dengan tinggi muka air, yaitu puncaknya bergeser ke arah hilir seiring dengan naiknya muka air, dan sebaliknya. Ketika tinggi muka air berada di bawah H 0 , perpindahan radial berkorelasi negatif dengan tinggi muka air, yaitu puncaknya bergeser ke arah hulu seiring dengan menurunnya muka air.
Alasannya adalah ketika ketinggian air lebih rendah ( H < H 0 ) (lih. Gambar 17 a), garis resultan aksi, yang diteruskan dari tekanan air hulu melalui lengkungan ke penopang, melewati hulu air. titik tengah permukaan pondasi penopang, timbul momen lentur berlawanan arah jarum jam, dan puncaknya bergeser ke arah hilir. Ketika ketinggian air lebih tinggi ( H > H 0 ) (lih. Gambar 17 b), garis aksi resultan melewati titik tengah permukaan pondasi penopang, momen lentur searah jarum jam dihasilkan, dan puncak dipindahkan ke hilir. .
Mirip dengan analisis perpindahan radial, kontribusi ketinggian air, suhu dan waktu terhadap amplitudo perpindahan tangensial tahunan (PL13-1, PL21) dihitung berdasarkan model HTT. Rata-rata kontribusi masing-masing komponen dalam 11 tahun terakhir (2007 hingga 2017) ditunjukkan pada Gambar 18 . Terlihat bahwa: (1) Pada amplitudo tangensial perpindahan tangensial penopang No.13 dan 21, komponen suhu masing-masing menyumbang rata-rata 96,52% dan 98,50%, dan komponen penuaan menyumbang 3,48% dan 1,50 % rata-rata, masing-masing. (2) Perpindahan tangensial hampir seluruhnya disebabkan oleh suhu.
Arah tekanan air tegak lurus terhadap perpindahan tangensial, dan bendungan telah beroperasi selama bertahun-tahun, sehingga tekanan air dan waktu tidak banyak berpengaruh terhadap perpindahan tangensial. Pengaruh ekspansi dan kontraksi termal merupakan penyebab mendasar perpindahan tangensial bendungan. Ketika suhu turun, lengkungan tersebut mengecil sehingga menarik penopang pada kedua tepian ke tengah, sehingga penopang pada tepi kiri bergerak ke arah tepi kanan, sedangkan penopang pada tepi kanan bergerak ke arah tepi kiri. Apabila suhu naik maka badan bendungan memanas dan mengembang sehingga mendorong penopang ke tepian, sehingga penopang di tepi kiri bergerak ke arah tepi kiri, sedangkan penopang di tepi kanan bergerak ke arah tepi kanan.
Hasil analisis di atas memberikan gambaran bahwa suhu mempunyai pengaruh paling besar terhadap perpindahan. Oleh karena itu, pengamatan harus ditingkatkan pada saat suhu cenderung meningkat atau menurun secara tiba-tiba. Karena efek superimposisi perpindahan tangensial, perhatian lebih harus diberikan pada lengkungan di tepi kanan dan kiri jika terjadi keretakan. Bendungan multi-lengkungan Foziling adalah jenis struktur yang ringan dan tipis, perubahan suhu yang drastis pasti akan menyebabkan tekanan suhu yang besar. Kemudian, simulasi keadaan tegangan bendungan pada suhu yang bervariasi akan menjadi sangat penting untuk mengevaluasi lebih lanjut keamanan bendungan.
Analisis numerik perilaku tegangan bendungan multi-lengkungan
Model numerik
Model numerik metode elemen hingga (FEM) merupakan teknik yang efektif untuk mempelajari perilaku bendungan dan menyelidiki keamanan bendungan. Seperti yang ditunjukkan pada Gambar 19 , jaring elemen hingga tiga dimensi dari bendungan multi-lengkung Foziling dibuat, yang mencakup 508.094 elemen dan 568.417 node. Dalam simulasi numerik, model elastis digunakan untuk sistem pondasi bendungan, dan kondisi keamanan tegangan elemen dinilai dibandingkan dengan kuat tarik ijin dan kuat tekan yang diperoleh berdasarkan kode relatif.
Model MLSR menunjukkan bahwa, dibandingkan dengan komponen tekanan dan suhu air, komponen penuaan sangat kecil dan cenderung stabil. Oleh karena itu, hanya variasi ketinggian air dan suhu yang dianggap sebagai beban utama setelah keadaan geostatis model.
Parameter beton badan bendungan dan batuan pondasi, yang menentukan kekuatan tegangan dan karakteristik deformasi, merupakan landasan penting untuk memperoleh hasil yang dapat diandalkan melalui pemodelan numerik. Parameter mekanik dan termodinamika yang sesuai 35 dalam simulasi FEM dirangkum dalam Tabel 3 .
Ketinggian air normal, tinggi muka banjir rencana, dan tinggi muka banjir periksa bendungan masing-masing adalah 125,56 m (di masa depan), 125,97 m, dan 129,97 m. Badan bendungan multi-lengkungan tipis dan bidang suhu terutama bergantung pada perubahan periodik suhu udara sekitar, sehingga tidak ada bidang suhu yang stabil di badan bendungan. Oleh karena itu, selisih antara suhu rata-rata bulanan maksimum dan suhu rata-rata tahunan selama periode operasi diambil sebagai beban kenaikan suhu. Selisih antara suhu rata-rata bulanan minimum dan suhu rata-rata tahunan diambil sebagai beban penurunan suhu. Prosedur simulasi numerik dapat dibagi menjadi tiga bagian: (1) Model perpindahan panas dibuat, di mana metode kesetaraan parameter 36digunakan untuk mensimulasikan lapisan isolasi termal. (2) Berdasarkan data yang diukur termasuk suhu air dan udara, bidang suhu tertinggi, bidang suhu terendah, dan bidang suhu kuasi-stabil diperoleh dengan menerapkan suhu batas yang berbeda. (3) Model ekspansi linier dibuat dengan mesh yang sama dengan model perpindahan panas. Setelah model mencapai keseimbangan di bawah gravitasi, tekanan air diterapkan pada permukaan hulu melalui pembebanan hidrostatik. Bidang suhu ditetapkan sebagai bidang yang telah ditentukan sebelumnya di setiap langkah penghitungan. Kemudian keamanan bendungan dievaluasi berdasarkan tegangan hasil simulasi numerik.
Analisis hasil numerik
Ketinggian air yang tinggi, suhu air yang tinggi, dan kenaikan dan penurunan suhu yang tiba-tiba merupakan kasus pembebanan yang merugikan bagi keamanan bendungan karena kerentanan suhu bendungan dan kapasitas penyimpanan yang kecil. Dalam simulasi, ketinggian penyimpanan normal pada 125,56 m, ketinggian banjir terkendali pada 129,96 m, serta suhu tinggi dan rendah dipertimbangkan. Di sini, suhu tinggi dan rendah mengacu pada rata-rata suhu maksimum dan minimum bulanan. Jadi, ada empat kombinasi pembebanan, yang tercantum pada Tabel 4 .
Diasumsikan tegangan tekan bernilai negatif dan tegangan tarik bernilai positif. Distribusi tegangan utama dalam kombinasi pembebanan (C1 dan C2) ditunjukkan pada Gambar. 20 dan 21 masing-masing, dimana penopang No.12 dan 13 diambil sebagai contoh. Hal tersebut dapat dilihat sebagai berikut:
-
Tegangan tarik di sebagian besar area lengkungan dan penopang lebih kecil dari tegangan tarik ijin (0,9 Mpa) beton bendungan di C1, kecuali untuk lengkungan No. 2 dan 22. Pada daerah tengah dan bawah penopang No. 12 dan 13 terjadi tegangan tarik. Karena keterbatasan pondasi, tegangan tarik maksimum pada penopang No. 12 dan 13 masing-masing mencapai 1,54 Mpa dan 1,55 Mpa (lihat Gambar 20 a), melebihi tegangan tarik ijin, namun lebih kecil dari tegangan tarik ultimit ( 3,53 MPa). Hal ini juga terjadi di bagian bawah lengkungan.
-
Tegangan tekan maksimum lengkungan No. 3–21 terutama terjadi pada ketinggian mendekati 1/3 kedalaman air di C1, dan nilai maksimumnya mencapai 3,54 Mpa, yang sangat mirip dengan distribusi tegangan tekan maksimum penopang. Tegangan tekan maksimum penopang No. 12 dan 13 masing-masing adalah 3,41 Mpa dan 3,46 Mpa (lihat Gambar 20 b), yang jauh lebih rendah daripada tegangan tekan izin (12,45 Mpa).
-
Tegangan tarik lengkungan dan penopang umumnya menurun pada C2 dibandingkan dengan C1. Tegangan tarik yang lebih besar terjadi pada bagian tengah panel penopang bagian hulu, dan tidak ada satupun yang melebihi tegangan tarik ijin. Terlebih lagi, konsentrasi tegangan pada bagian kontak antara lengkungan dan pondasi menghilang, namun masih ada pada kaki penopang (lihat Gambar 21 a). Tegangan tekan maksimum penopang No. 12 dan 13 masing-masing sebesar 4,26 Mpa dan 4,14 Mpa.
-
Tegangan tarik yang lebih besar terjadi pada bagian tengah puncak lengkungan No. 2 dan 22 pada C1 dan C2. Dan nilai maksimumnya masing-masing mencapai 1,46 Mpa dan 1,27 Mpa. Di dekat bagian bawah sisi lengkungan No. 2 dan 22 yang dihubungkan dengan bendungan gravitasi, tegangan tarik yang lebih besar (1,52 Mpa) terjadi pada C1, sedangkan tegangan tekan yang lebih besar (4,14 Mpa) terjadi pada C2.
Kasus beban khusus juga diterapkan pada bendungan di C3 dan C4. Distribusi tegangan utama sangat mirip dengan C1 dan C2, namun nilai tegangan umumnya meningkat. Hal ini tidak terulang di sini. Tegangan tarik maksimum dan tegangan tekan maksimum tercantum pada Tabel 5 .
Akibat besarnya tegangan tarik yang terjadi di dekat puncak lengkungan No. 2 dan 22, serta pondasi lengkungan dan penopang di C1 dan C3, melebihi kuat tarik yang diijinkan, maka kasus beban tinggi muka air + suhu Drop dianggap sebagai kombinasi beban operasi yang paling tidak menguntungkan bagi keamanan bendungan.
Dibandingkan dengan tegangan yang disimulasikan sebelum perkuatan (35) , tegangan umumnya berkurang, yang menunjukkan bahwa keamanan bendungan meningkat. Area yang melebihi tegangan tarik yang diijinkan dianggap sebagai area rawan kerusakan, yang ditunjukkan dalam zona persegi panjang (lihat Gambar 20 dan 21). Karena sensitivitas bendungan multi-lengkungan terhadap suhu dan efek kumulatif deformasi, perpindahan tangensial penopang No. 2 dan No. 22 lebih besar daripada yang lain. Selain itu, salah satu ujung lengkungan No. 2 dan 22 dihubungkan dengan bendungan gravitasi, dan ujung lainnya dihubungkan dengan penopang. Perbedaan kekakuan di antara keduanya begitu besar sehingga ujung yang terhubung ke bendungan gravitasi lebih terkekang dan terkena gaya yang lebih besar. Oleh karena itu, jika ketinggian air dan suhu terus berubah, retakan lelah cenderung terjadi pada lengkungan No. 2 dan 22. Akibatnya, lengkungan No. 2 dan 22 serta bagian bawah lengkungan dan penopang dianggap sebagai bagian lemah dari bendungan, dan perhatian lebih harus diberikan dalam patroli keselamatan.
Oleh karena itu, area lengkungan No. 2 dan 22 dari pondasi hingga ketinggian 127,46 m diperkuat dengan beton serat baja sebagai tulangannya. Area permukaan hulu dan hilir lengkungan serta perpanjangan kiri dan kanannya selebar 0,5 m, masing-masing dipertebal sebesar 0,15 m dan 0,1 m. Untuk daerah tegangan tarik yang besar pada muka beton bagian hulu dan dinding samping tiap penopang, dilakukan penyemprotan beton serat baja untuk penebalan dari bagian dalam penopang. Muka beton bagian hulu dipertebal 0,5 m, muka beton bagian bawah (di bawah elevasi 79,56 m) dipertebal 0,7 m, dan dinding samping kiri dan kanan dipertebal 0,4 m. Apalagi kuat tekan beton serat baja 48,0Mpa, kuat tarik 3,5Mpa, dan kuat rekat dengan beton lama 2,5–3,0Mpa 35. Proyek perkuatan ini dimaksudkan untuk mencegah keretakan.
Kesimpulan
Dalam studi ini, deformasi bendungan multi-lengkung Foziling dikarakterisasi berdasarkan data pemantauan. Kontribusi ketinggian air, suhu, dan waktu terhadap perpindahan radial bendungan diukur dengan model HTT dan metode MLSR. Kemudian dianalisis pengaruh karakteristik tekanan dan suhu air terhadap perpindahan. Akhirnya, model numerik dibuat untuk menyelidiki keadaan tegangan bendungan dalam berbagai kasus beban. Kesimpulan utama yang dapat diperoleh dari penelitian ini adalah sebagai berikut:
-
Bendungan multi-lengkung Foziling berada dalam keadaan elastis setelah perkuatan terakhir berdasarkan kinerja MLSR yang baik pada titik-titik pemantauan perpindahan tipikal. Pengaruh karakteristik ketinggian air dan suhu terhadap perpindahan bendungan diinterpretasikan dengan memisahkan variabel berdasarkan model HTT. Penopang dipindahkan ke hilir (hulu) ketika suhu naik (turun). Ketika tinggi muka air naik dari rendah ke tinggi, hubungan antara perpindahan dan tinggi muka air mula-mula bersifat negatif dan kemudian positif.
-
Suhu memberikan kontribusi paling besar terhadap perpindahan radial penopang, diikuti oleh ketinggian air, dan waktu. Perpindahan tangensial hampir seluruhnya disebabkan oleh suhu. Dapat disimpulkan bahwa bendungan multi-lengkungan cukup sensitif terhadap fluktuasi suhu, dan fluktuasi suhu yang drastis merupakan faktor yang merugikan keamanan bendungan.
-
Tegangan tarik yang lebih besar terutama terjadi di dekat puncak lengkungan No.2 dan 22 serta bagian bawah lengkungan dan penopang, terutama pada kasus beban suhu rendah dan ketinggian air tinggi. Selain itu, kerusakan akibat kelelahan dapat terjadi pada area dengan variasi tegangan yang besar, yang mungkin menyebabkan timbulnya retakan. Meskipun penerapan beton bertulang serat baja di wilayah tegangan tinggi dapat secara efektif mengurangi timbulnya retakan lelah, perhatian lebih harus diberikan pada bagian bendungan yang lemah.
-
Hasil analisis perpindahan dan tegangan menggambarkan bahwa bendungan masih beroperasi dengan baik sejak perkuatan terakhir, dan efektivitas proyek perkuatan, memberikan landasan yang baik untuk pengoperasian di masa depan pada permukaan air yang lebih tinggi. Di sisi lain, hal ini juga menunjukkan bahwa metode yang diusulkan dalam makalah ini memiliki nilai untuk mempopulerkan teknik dan penerapan pada bendungan beton lainnya.
Referensi
Wu, Teori Pemantauan Keamanan ZR dan Penerapannya pada Struktur Hidraulik (Higer Education Press, 2003).
Jeon, J., Lee, J., Shin, D. & Park, H. Pengembangan sistem manajemen keamanan bendungan. Adv. bahasa Inggris perangkat lunak. 40 , 554–563. https://doi.org/10.1016/j.advengsoft.2008.10.009 (2009).
Loh, C.-H., Chen, C.-H. & Hsu, T.-Y. Penerapan metode statistik tingkat lanjut untuk mengekstraksi tren jangka panjang dalam data pemantauan statis dari bendungan lengkung. Struktur. Pemantauan Kesehatan. 10 , 587–601 (2011).
Kang, F., Li, J., Zhao, S. & Wang, Y. Pemantauan kesehatan struktural bendungan beton menggunakan suhu udara jangka panjang untuk simulasi efek termal. bahasa Inggris Struktur. 180 , 642–653. https://doi.org/10.1016/j.engstruct.2018.11.065 (2019).
Mata, J. Interpretasi perilaku bendungan beton dengan jaringan saraf tiruan dan model regresi linier berganda. bahasa Inggris Struktur. 33 , 903–910. https://doi.org/10.1016/j.engstruct.2010.12.011 (2011).
Salazar, F., Morán, R., Toledo, M. Á. & Oñate, E. Model berbasis data untuk prediksi perilaku bendungan: Tinjauan dan beberapa pertimbangan metodologis. Lengkungan. Hitung. Metode E 24 , 1–21. https://doi.org/10.1007/s11831-015-9157-9 (2015).
Tatin, M., Briffaut, M., Dufour, F., Simon, A. & Fabre, JP Perpindahan termal bendungan beton: Menghitung suhu air dalam model statistik. bahasa Inggris Struktur. 91 , 26–39. https://doi.org/10.1016/j.engstruct.2015.01.047 (2015).
Su, H., Li, X., Yang, B. & Wen, Z. Wavelet mendukung model prediksi deformasi bendungan berbasis mesin vektor. Mekanisme. sistem. Proses Sinyal. 110 , 412–427. https://doi.org/10.1016/j.ymssp.2018.03.022 (2018).
Stojanovic, B., Milivojevic, M., Ivanovic, M., Milivojevic, N. & Divac, D. Sistem adaptif untuk pemodelan perilaku bendungan berdasarkan regresi linier dan algoritma genetika. Adv. bahasa Inggris perangkat lunak. 65 , 182–190. https://doi.org/10.1016/j.advengsoft.2013.06.019 (2013).
Salazar, F., Toledo, MA, Oñate, E. & Morán, R. Perbandingan empiris teknik pembelajaran mesin untuk pemodelan perilaku bendungan. Struktur. Saf. 56 , 9–17. https://doi.org/10.1016/j.strusafe.2015.05.001 (2015).
Yu, Y., Liu, X., Wang, E., Fang, K. & Huang, L. Evaluasi keamanan bendungan berdasarkan regresi linier berganda dan simulasi numerik. Mekanisme Batu. Batu Eng. 51 , 2451–2467. https://doi.org/10.1007/s00603-018-1435-z (2018).
Xi, G.-Y., Yue, J.-P., Zhou, B.-X. & Tang, P. Penerapan algoritma imun buatan pada model statistik perpindahan bendungan. Hitung. Matematika. Aplikasi. 62 , 3980–3986. https://doi.org/10.1016/j.camwa.2011.09.057 (2011).
Stojanovic, B., Milivojevic, M., Milivojevic, N. & Antonijevic, D. Sistem penyesuaian mandiri untuk pemodelan perilaku bendungan berdasarkan jaringan saraf tiruan yang berkembang. Adv. bahasa Inggris perangkat lunak. 97 , 85–95. https://doi.org/10.1016/j.advengsoft.2016.02.010 (2016).
Kao, C.-Y. & Loh, C.-H. Pemantauan data deformasi statis jangka panjang bendungan lengkung Fei-Tsui menggunakan pendekatan berbasis jaringan syaraf tiruan. Struktur. Kontrol Kesehatan 20 , 282–303. https://doi.org/10.1002/stc.492 (2013).
Wei, B., Yuan, D., Li, H. & Xu, Z. Model prakiraan kombinasi untuk perpindahan bendungan beton dengan mempertimbangkan koreksi sisa. Struktur. Pemantauan Kesehatan. 18 , 232–244. https://doi.org/10.1177/1475921717748608 (2017).
Nourani, V. & Babakhani, A. Integrasi jaringan saraf tiruan dengan interpolasi fungsi basis radial dalam pemodelan rembesan bendungan timbunan tanah. J.Komputasi. Teknik Sipil. 27 , 183–195. https://doi.org/10.1061/(asce)cp.1943-5487.0000200 (2013).
Cheng, L. & Zheng, D. Dua model pemantauan keamanan bendungan online berdasarkan proses penggalian dampak lingkungan. Adv. bahasa Inggris perangkat lunak. 57 , 48–56. https://doi.org/10.1016/j.advengsoft.2012.11.015 (2013).
Ranković, V., Grujović, N., Divac, D. & Milivojević, N. Pengembangan model identifikasi regresi vektor pendukung untuk prediksi perilaku struktural bendungan. Struktur. Saf. 48 , 33–39. https://doi.org/10.1016/j.strusafe.2014.02.004 (2014).
Bui, K.-TT, Tien Bui, D., Zou, J., Van Doan, C. & Revhaug, I. Pendekatan kecerdasan buatan hibrid baru berdasarkan model inferensi fuzzy saraf dan optimasi kawanan partikel untuk pemodelan perpindahan horizontal pembangkit listrik tenaga air bendungan. Komputasi Neural. Aplikasi. 29 , 1495–1506. https://doi.org/10.1007/s00521-016-2666-0 (2016).
Salazar, F., Toledo, M. Á., González, JM & Oñate, E. Deteksi dini anomali dalam kinerja bendungan: Metodologi yang didasarkan pada pohon regresi yang ditingkatkan. Struktur. Kontrol Kesehatan. https://doi.org/10.1002/stc.2012 (2017).
Kang, F., Liu, J., Li, J. & Li, S. Model prediksi deformasi bendungan beton untuk pemantauan kesehatan berdasarkan mesin pembelajaran ekstrim. Struktur. Kontrol Kesehatan. https://doi.org/10.1002/stc.1997 (2017).
Dai, B., Gu, C., Zhao, E. & Qin, X. Model statistik mengoptimalkan model regresi hutan acak untuk pemantauan deformasi bendungan beton. Struktur. Kontrol Kesehatan. https://doi.org/10.1002/stc.2170 (2018).
Yuan, D., Gu, C., Qin, X., Shao, C. & He, J. Model DHM bendungan lengkung super tinggi berbasis TSVR yang ditingkatkan kinerjanya menggunakan suhu udara terukur. bahasa Inggris Struktur. https://doi.org/10.1016/j.engstruct.2021.113400 (2022).
Ribeiro, LS, Wilhelm, VE, Faria, É. F., Correa, JM & dos Santos, ACP Analisis komparatif model deformasi beton jangka panjang bendungan penopang. bahasa Inggris Struktur. 193 , 301–307. https://doi.org/10.1016/j.engstruct.2019.05.043 (2019).
Qu, X., Yang, J. & Chang, M. Model pembelajaran mendalam untuk prediksi deformasi bendungan beton berdasarkan RS-LSTM. J. Sens.1–14 , 2019. https://doi.org/10.1155/2019/4581672 (2019 ) .
Liu, W., Pan, J., Ren, Y., Wu, Z. & Wang, J. Model prediksi kopling untuk perpindahan bendungan lengkung jangka panjang berdasarkan jaringan memori jangka pendek. Struktur. Kontrol Kesehatan. https://doi.org/10.1002/stc.2548 (2020).
Yang, D. dkk. Metode prediksi deformasi bendungan beton berbasis LSTM dengan mekanisme perhatian. Akses IEEE 8 , 185177–185186. https://doi.org/10.1109/access.2020.3029562 (2020).
Léger, P. & Leclerc, M. Hidrostatik, suhu, model perpindahan waktu untuk bendungan beton. J.Eng. Mekanisme. 133 , 267–277 (2007).
Mata, J., Tavares de Castro, A. & Sá da Costa, J. Membangun model statistik untuk deformasi bendungan lengkung. Struktur. Kontrol Kesehatan 21 , 423–437. https://doi.org/10.1002/stc.1575 (2014).
Valente, S., He, Q. & Capriulo, C. Analisis Seismik Bendungan Beton Datar Pinus (Analisis Numerik Bendungan, 2021).
Valente, S., Capriulo, C. & He, Q. Retakan geser kompresi yang tumbuh pada bendungan gravitasi lengkung. teori. Aplikasi. pecahan. Mekanisme. 101 , 178–184. https://doi.org/10.1016/j.tafmec.2019.02.011 (2019).
Gan, L., Shen, Z.-Z. & Xu, L.-Q. Analisis deformasi jangka panjang pada bendungan timbunan batu berwajah beton Jiudianxia. Arab. J.Ilmu. bahasa Inggris 39 , 1589–1598. https://doi.org/10.1007/s13369-013-0788-6 (2013).
Hariri-Ardebili, MA & SEED-Kolbadi, SM Retakan seismik dan ketidakstabilan bendungan beton: Pendekatan retakan oles. bahasa Inggris Gagal. Dubur. 52 , 45–60. https://doi.org/10.1016/j.engfailanal.2015.02.020 (2015).
Wang, Z., Liu, S., Vallejo, L. & Wang, L. Analisis numerik penyebab retakan pelat muka di bendungan timbunan batu Gongboxia. bahasa Inggris Geol. 181 , 224–232. https://doi.org/10.1016/j.enggeo.2014.07.019 (2014).
Shi, Q. Penguatan Waduk Foziling (China Water & Power Press, 2008).
Wang, Q., Wang, J.-T. & Jin, F. Simulasi setara isolator panas bendungan dan analisis efeknya. Adv. Sains. Teknologi. Sumber Daya Air. 27 , 4 (2007) ( (dalam bahasa Cina) ).
Penutup
Sekian Penjelasan Singkat Mengenai Bendungan multi lengkungan (Foziling). Semoga Bisa Menambah Pengetahuan Kita Semua.