余勁松 安浩然 郭來棟

以某3層框架結(jié)構(gòu)為研究對象，基于流固耦合分析理論，采用ANSYS構(gòu)建了山洪與結(jié)構(gòu)的耦合分析模型，分析結(jié)構(gòu)在不同深度洪水沖擊時的響應(yīng)情況，同時比選分析了“過流堰”防護(hù)措施、短柱擾流防護(hù)措施下結(jié)構(gòu)的響應(yīng)情況。研究結(jié)果表明：山洪沖擊結(jié)構(gòu)時，水深未達(dá)到底層樓面時，結(jié)構(gòu)位移及應(yīng)力均較小，滿足結(jié)構(gòu)安全要求，到達(dá)樓面板后，結(jié)構(gòu)最大應(yīng)力突增，迎水側(cè)底層柱柱底混凝土開裂；“過流堰”防護(hù)措施可以顯著降低結(jié)構(gòu)穩(wěn)定繞流階段的穩(wěn)定位移，但應(yīng)力降低不明顯，迎水側(cè)混凝土仍有開裂風(fēng)險；短柱擾流防護(hù)措施能夠降低結(jié)構(gòu)應(yīng)力，滿足結(jié)構(gòu)安全要求，可以考慮在山區(qū)新建房屋時采取此種方式。相關(guān)經(jīng)驗可供類似山區(qū)低層框架結(jié)構(gòu)抗洪設(shè)計借鑒。

洪水沖擊； 山區(qū)房屋； 流固耦合； 動力響應(yīng)； 防護(hù)措施

TV131.2 A

［定稿日期］2022-03-31

［作者簡介］余勁松（1996—），男，碩士，研究方向為橋梁防災(zāi)減災(zāi)。

洪澇是因降雨、融雪、冰凌、潰壩（堤）、風(fēng)暴潮、熱帶氣旋等造成的江河洪水、漬澇、山洪、滑坡和泥石流等，以及由其引發(fā)的次生災(zāi)害。山洪災(zāi)害是我國洪澇災(zāi)害引起的主要次生災(zāi)害，危害極大，2020年，全國（港澳臺除外）因洪澇受災(zāi)7 861.5萬人次，倒塌房屋9.0萬間，造成直接經(jīng)濟(jì)損失2 669.8億元［1］。山洪具有突發(fā)性強、破壞力大、難以預(yù)警、發(fā)生在偏遠(yuǎn)山區(qū)等特點。國內(nèi)外學(xué)者對山洪及其引起的次生災(zāi)害進(jìn)行了大量研究，除了對暴雨山洪的形成原因和演進(jìn)規(guī)律進(jìn)行研究外［2-5］，也有學(xué)者對洪水、泥石流等對結(jié)構(gòu)的影響，及其防治措施進(jìn)行了研究。肖詩云等［6-10］通過水槽試驗，研究了到潰口不同距離、不同墻體開洞率下的山區(qū)鄉(xiāng)村房屋抗洪能力。楊哲豪等［11-12］建立了數(shù)值模型研究潰壩洪水對結(jié)構(gòu)物的沖擊荷載。王友彪等［13-15］通過水槽試驗和數(shù)值模擬的方法，以流深和流速作為主要研究參數(shù)，研究了洪水、泥石流對橋墩的影響。橋梁工程中，常采用設(shè)置擾流短柱的方式減小洪水對橋墩基礎(chǔ)的沖刷和沖擊［16］。而目前針對山區(qū)房屋防山洪的工程措施的研究還很少［17-18］。

目前，對于橋梁、大壩等結(jié)構(gòu)在洪水或泥石流作用下的動力響應(yīng)研究水平不斷提高，研究內(nèi)容也不斷深入，但對于山區(qū)房屋，研究相對滯后。在我國西部山區(qū)，山洪影響更為突出，山區(qū)房屋多為低層，常采用整體框架的結(jié)構(gòu)形式。在地質(zhì)災(zāi)害多發(fā)地，房屋常面臨洪水和泥石流的威脅，因此研究低層框架結(jié)構(gòu)在山洪沖擊作用下的動力響應(yīng)很有必要。

本文依托某洪水易發(fā)地3層框架結(jié)構(gòu)民房開展研究，結(jié)構(gòu)首層層高4 m，2～3層層高3 m。為確保該房屋在洪水沖擊作用后的正常使用，保障人民生命財產(chǎn)安全，本文通過數(shù)值模擬計算，研究在一定洪水流速沖擊作用下，不同水深工況下洪水對結(jié)構(gòu)的影響，通過分析結(jié)構(gòu)的位移及應(yīng)力，提出合理的防護(hù)措施，以期對相似民建工程提供參考。

1 流固耦合模型與方法

1.1 原理

山洪沖擊房屋屬于流固耦合作用問題，即建筑結(jié)構(gòu)在流場作用下的振動和變形會反過來改變流場，改變后的流場又反作用于結(jié)構(gòu)，導(dǎo)致結(jié)構(gòu)的荷載幅值及分布發(fā)生改變。流固耦合過程中結(jié)構(gòu)的動力學(xué)控制方程基于達(dá)朗貝爾原理（DAlemberts principle）見式（1）。

［M］{}+［C］{a·}+［K］{a}={P（t）}（1）

式中，［M］為結(jié)構(gòu)質(zhì)量矩陣，［C］為阻尼矩陣，［K］為剛度矩陣；{}、{a·}、{a}分別為結(jié)構(gòu)的加速度、速度及位移；{P（t）}為結(jié)構(gòu)所受合力，t為結(jié)構(gòu)與流體耦合時間參數(shù)。

流體介質(zhì)控制方程為質(zhì)量守恒方程、動量守恒方程和能量守恒方程。其通用形式表示為式（2）。

（ρφ）t+（ρuφ）x+（ρvφ）y+（ρwφ）z=

xΓφx+yΓφy+zΓφz+S（2）

式中：t是時間；φ是通用變量，例如，表示為質(zhì)量守恒方程時取1，表示動量守恒方程時取速度量；u、v、w分別表示速度矢量在x、y、z方向的分量；ρ為計算域流體密度；Г為廣義擴散系數(shù)；S為廣義源項。

以水作為流體介質(zhì)，考慮其為黏性不可壓縮流體，不考慮熱量交換。采用有限體積法對控制方程和流體區(qū)域進(jìn)行離散，結(jié)構(gòu)域和流體域的計算結(jié)果通過流固耦合面互相傳遞。待當(dāng)前時刻的計算結(jié)果迭代收斂后，進(jìn)行下一個時間步的計算。具體計算過程通過調(diào)用ANSYS- FLUENT的瞬態(tài)計算模式實現(xiàn)。

1.2 數(shù)值模型

分析對象為3層框架結(jié)構(gòu)，如圖1所示。首層層高4 m，2、3層層高3 m，樓板厚度0.1 m，立柱0.4 m×0.4 m，梁高0.5 m、寬0.3 m。梁板柱均采用C30混凝土，抗拉強度ft=1.43 MPa，彈性模量E=3×104 MPa，泊松比μs=0.18，容重G=23 kN/m3?；炷帘緲?gòu)關(guān)系采用FLUENT默認(rèn)本構(gòu)。山區(qū)洪水含沙量較大［12］，會少量增加水的密度及黏度。本模型忽略含沙量對洪水的影響，采用黏性不可壓縮的清水作為流體材料。密度ρ=1000 kg/m3，動力粘度μ=1.003×10-3 Pa·s。

設(shè)置洪水流域為矩形斷面，為減弱計算域壁面對模型流場的影響，設(shè)置計算域?qū)?0 m、長50 m，高為h （h=1～4 m），同時設(shè)置邊界層10層，每層5 mm；河流洪水流速范圍一般為4～15 m/s（這是對于常規(guī)河流流域洪水而言，河流斷面常呈“V”型，河流寬度通常在20～50 m），對于山區(qū)聚居地，山洪演進(jìn)過程中，洪水?dāng)嗝娓?，能量被地表植被、低矮建筑耗散較大，且山區(qū)山谷交錯。因此，山洪洪峰到達(dá)民建建筑時其能量已大幅減少，在沖擊結(jié)構(gòu)時流速相對于河流流速較小，故偏危險的將其流速考慮為4 m/s。設(shè)置進(jìn)口（inlet）邊界為速度條件4 m/s；地面、側(cè)面及流固耦合面設(shè)置為壁面（Wall）；頂面及出口設(shè)置為開放邊界（Opening）且相對壓強為0。流固耦合模型如圖2所示。

數(shù)值計算采用標(biāo)準(zhǔn)k-epsilon兩方程模型，近壁面函數(shù)采用標(biāo)準(zhǔn)壁面函數(shù)。

2 計算結(jié)果分析

2.1 結(jié)構(gòu)響應(yīng)分析

在一定流速下，洪水沖擊結(jié)構(gòu)的位移、應(yīng)力響應(yīng)與水位有關(guān)。入口流速恒定為4 m/s，計算域水深為1～4 m，得到結(jié)構(gòu)承受的沖擊力時程曲線、上游側(cè)柱底迎水面正應(yīng)力時程曲線及底層柱柱頂位移時程曲線，如圖3～圖5所示。

流體沖擊結(jié)構(gòu)一般可分為瞬時沖擊、減速衰減、穩(wěn)定繞流3個階段。如圖3所示，水深為4 m時，結(jié)構(gòu)所受沖擊力在0.2 s內(nèi)從A點的0迅速增加到B點的398.7 kN，結(jié)合圖4、圖5亦可發(fā)現(xiàn)該階段的結(jié)構(gòu)應(yīng)力及位移響應(yīng)值迅速增大，為瞬時沖擊階段。此外，水流與結(jié)構(gòu)物接觸的瞬間，速度變化率達(dá)到最大值，故其沖擊力達(dá)到最大值。在B點到C點間，不難看出其沖擊力逐漸衰減，圖像均呈下凹狀，沖擊力衰減速率降低，沖擊力平均衰減速率為64.13 kN/s，在B點時結(jié)構(gòu)的自身能量（勢能）達(dá)到最大值，因此結(jié)構(gòu)出現(xiàn)回彈，位移減小，做往復(fù)運動，出現(xiàn)振動現(xiàn)象。由于結(jié)構(gòu)與流場的流固耦合作用，水流會耗散結(jié)構(gòu)的彈性勢能及動能，對其產(chǎn)生阻尼力，隨著結(jié)構(gòu)總能量的減小，這種阻尼力也隨之減小，因此結(jié)構(gòu)響應(yīng)在峰值點后呈現(xiàn)減速衰減的現(xiàn)象。C點過后，當(dāng)結(jié)構(gòu)最終動能及勢能被水流耗散，結(jié)構(gòu)自身產(chǎn)生一個基本穩(wěn)定的變形，其回彈力與水流沖擊力平衡，該階段為穩(wěn)定繞流階段。在本模型中，最大變形值僅為0.21 mm，因此不采用位移作為結(jié)構(gòu)安全評估因素。

由圖4～圖6可知，水深為1～3 m時，結(jié)構(gòu)最大應(yīng)力出現(xiàn)在上游側(cè)柱底迎水面，分別為0.05 MPa、0.27 MPa、0.58 MPa，最大位移分別為0.005 mm、0.072 mm、0.210 mm。當(dāng)水深為4 m時，結(jié)構(gòu)最大位移、應(yīng)力分別為0.157 mm、2.12 MPa，結(jié)構(gòu)最大拉應(yīng)力出現(xiàn)在上游側(cè)柱底迎水面，此工況下上游側(cè)柱底穩(wěn)定繞流下的拉應(yīng)力（1.68 MPa）大于混凝土抗拉強度設(shè)計值，柱底出現(xiàn)拉裂縫。不難看出，水深為4 m時結(jié)構(gòu)應(yīng)力響應(yīng)較3 m水深大1.54 MPa，但其最大位移響應(yīng)小于3 m水深時結(jié)構(gòu)最大位移，2個工況在穩(wěn)定繞流階段的位移響應(yīng)值接近。通過分析認(rèn)為，當(dāng)洪水洪峰沖擊到柱迎水面及樓面板的瞬間，接觸面水流流速突變?yōu)?，其加速度達(dá)到最大值，結(jié)構(gòu)產(chǎn)生最大響應(yīng)，該階段為瞬時沖擊階段。洪水沖擊樓面板，樓板會增大結(jié)構(gòu)的抗側(cè)力剛度，直接導(dǎo)致結(jié)構(gòu)的順流向位移降低，樓板具有較大的阻水面積且豎向剛度較低，導(dǎo)致其產(chǎn)生向上的位移，洪水對面板有一個向上的舉托力，樓面板產(chǎn)生向上翹曲，最大位移出現(xiàn)在樓面。同時，該舉托力導(dǎo)致柱底拉應(yīng)力加大，限制了結(jié)構(gòu)沿水流方向的位移，因此水深4 m時結(jié)構(gòu)應(yīng)力較大，位移卻較小。

綜上所述，在一定的流速下，當(dāng)水深低于首層樓面高度時，結(jié)構(gòu)所受的沖擊力和位移響應(yīng)都隨著洪水深度的增大而增大。當(dāng)水深達(dá)到首層樓面高度，洪水直接作用于梁和樓面板時，結(jié)構(gòu)受到的沖擊力急劇增加，上游側(cè)柱底處的應(yīng)力也大幅增加。但是由于樓面板對水流的影響，結(jié)構(gòu)的位移響應(yīng)較小。可以看出，在保證結(jié)構(gòu)安全性的基礎(chǔ)上，適當(dāng)增加結(jié)構(gòu)首層層高，可以避免結(jié)構(gòu)薄弱位置出現(xiàn)較大拉應(yīng)力，有助于防止該處混凝土出現(xiàn)拉裂縫。

2.2 結(jié)構(gòu)表面沖擊荷載分析

研究結(jié)構(gòu)受到最大沖擊力時，提取上游側(cè)底層中柱迎流面壓應(yīng)力分布，如圖7所示。

從圖7容易看出，h=1～3 m時，洪水對柱體的最小沖擊壓力均出現(xiàn)在h處，最小沖擊壓力分別為 1 621 Pa、1 643 Pa、1 690 Pa，隨著洪水深度的增加，其對結(jié)構(gòu)的最小沖擊壓力雖有增長，但增加幅度不大，基本保持不變。h=1～3 m時，水對柱體的最大沖擊壓力出現(xiàn)在約0.5h的位置，分別為5 400 Pa、5 950 Pa、6 100 Pa。同樣的，隨著洪水深度的增加，其對結(jié)構(gòu)的最大沖擊壓力的增加幅度不大。

h=4 m時，3 m以下的柱表面沖擊壓力在1 000 Pa左右，大于3 m時，隨著高度的增加，柱表面所受沖擊壓力持續(xù)增加，最大沖擊壓力出現(xiàn)在柱子最高處約為3 448.1 Pa，較其余3個工況小。分析其原因，水深h=4 m時，洪水已經(jīng)淹沒底層樓面，樓面和梁分擔(dān)了部分洪水沖擊力。同時由于樓板的阻水效應(yīng)，導(dǎo)致部分水流回流，與來流的動能相互抵消，因此其柱子的最大表面沖擊壓應(yīng)力和最小表面沖擊壓應(yīng)力均較小。

綜上所述，當(dāng)洪水水深低于樓面板板高度時，柱體表面所受的沖擊壓力分布呈倒“C”型。最小沖擊壓力出現(xiàn)在水深高度處，最大沖擊壓力出現(xiàn)在0.5倍水深處。而且洪水水深對沖擊壓力的最大值、最小值影響較小。當(dāng)洪水水深剛剛漫過樓面板時，靠近水面部分的柱體表面所受的沖擊壓力較大，其余部分沖擊壓力較小，分布也較均勻。

3 山洪防護(hù)措施分析

由前文可知，當(dāng)水深h=4 m時，結(jié)構(gòu)最大拉應(yīng)力為2.12 MPa，穩(wěn)定繞流階段的拉應(yīng)力為1.67 MPa，都出現(xiàn)在迎流側(cè)柱底，均超過混凝土抗拉強度設(shè)計值1.43 MPa。因此，將針對水深h=4 m的情形，對該3層民用建筑進(jìn)行了山洪防護(hù)設(shè)計。

3.1 防護(hù)措施

山洪防護(hù)措施主要包括非工程措施和工程措施［10］。非工程措施包括調(diào)查評價、監(jiān)測預(yù)警系統(tǒng)、群測群防體系及加強山丘區(qū)人員搬遷管理，為主要防護(hù)措施。工程措施為被動防護(hù)方式，采取護(hù)岸、堤防，河道護(hù)底、跌水等消能措施，但此類防御措施范圍廣、投資大，且在水流長期作用下防護(hù)結(jié)構(gòu)容易出現(xiàn)基底脫空，混凝土斷裂，結(jié)構(gòu)沖毀等問題。因此，一種簡單、高效、投資小的防護(hù)措施對改善山區(qū)山洪沖擊影響具有積極意義。

本文提出2種適用于普通民宅的洪水防護(hù)方式——增設(shè)小型“過流堰”或短柱，如圖8所示。因山村地區(qū)部分民宅門外會修建一個小斜坡，以保證室內(nèi)外高差來防潮。當(dāng)洪水沖擊時，若漫過斜坡，該斜坡即為過流堰，第一種洪水防護(hù)措施即來源于此。此模型設(shè)置的“過流堰”截面尺寸為3 m×0.5 m見圖8（a）。第二種防護(hù)措施來源于橋梁結(jié)構(gòu)橋墩的防護(hù)方式——在柱（墩）前增設(shè)擾流短柱見圖8（b）。短柱與房屋立柱的中心距為1 m，直徑為0.4 m，高1 m。短柱可以替房屋立柱抵擋洪水沖擊，改變立柱周圍的流場，從而改變結(jié)構(gòu)的受力特性。

3.2 計算結(jié)果分析

增加防護(hù)裝置后的模型與前文h=4 m進(jìn)行相同的邊界條件設(shè)置，得出底層柱柱頂位移和上游側(cè)中柱柱底表面應(yīng)力分別如圖9、圖10所示。

從圖9可知，結(jié)構(gòu)達(dá)到穩(wěn)定繞流階段后，未增加防護(hù)措施下的位移響應(yīng)值大于增設(shè)防護(hù)措施的位移響應(yīng)，分別為0.134 mm、0.107 mm、0.095 mm，三者最大位移響應(yīng)值接近，“過流堰”防護(hù)和短柱防護(hù)分別減小了穩(wěn)定繞流階段位移20%、29%。雖然防護(hù)裝置對最大位移響應(yīng)值影響不大，但其對降低穩(wěn)定繞流階段的穩(wěn)定位移值具有明顯效果。

從圖10中可以看出，未設(shè)防護(hù)裝置與過流堰防護(hù)所得柱底最大拉應(yīng)力基本相同，為2.10 MPa，短柱防護(hù)所得柱底最大拉應(yīng)力為1.85 MPa，較未防護(hù)情況減小13%。因瞬時沖擊拉應(yīng)力出現(xiàn)時間較短，且較混凝土極限應(yīng)力強度相近，因此認(rèn)為洪水來流時瞬時沖擊力并不會引起結(jié)構(gòu)出現(xiàn)致命損傷，主要以穩(wěn)定繞流階段的應(yīng)力進(jìn)行評估結(jié)構(gòu)安全性。未加防護(hù)措施、過流堰、短柱的在穩(wěn)定繞流階段的柱底拉應(yīng)力分別為1.67 MPa、1.61 MPa、1.34 MPa。根據(jù)防護(hù)裝置對柱底拉應(yīng)力的減弱程度，判斷其防護(hù)效率，過流堰和短柱分別減小柱底拉應(yīng)力4%、19%。短柱的防護(hù)能力明顯高于過流堰防護(hù)，而且短柱防護(hù)使柱底承受最大拉應(yīng)力低于混凝土極限抗拉強度，有效避免了混凝土開裂。

綜上所述，過流堰防護(hù)對結(jié)構(gòu)位移及應(yīng)力響應(yīng)影響較小，防護(hù)效果不明顯，不能避免混凝土開裂。而增設(shè)擾流短柱能夠有效減少結(jié)構(gòu)位移，減小混凝土表面應(yīng)力，有效避免出現(xiàn)拉裂縫，有很好的防護(hù)效果。

4 結(jié)束語

本文以某山區(qū)3層框架結(jié)構(gòu)作為研究對象，基于流固耦合分析理論，采用數(shù)值模擬的方法，對結(jié)構(gòu)在不同深度洪水沖擊時的安全性，和2種防護(hù)措施進(jìn)行了研究，得到結(jié)論：

（1）當(dāng)水深未達(dá)到底層樓面時，隨著洪水深度的增加，洪水沖擊力增大，結(jié)構(gòu)的位移和應(yīng)力隨之增加。但當(dāng)洪水漫過底層樓面時，由于樓面對水流的作用，結(jié)構(gòu)位移反而較小。但是結(jié)構(gòu)受到的洪水沖擊力和應(yīng)力響應(yīng)會急劇加大，威脅結(jié)構(gòu)安全。因此可以適當(dāng)增加首層層高，提高結(jié)構(gòu)的抗洪能力。

（2）當(dāng)水深未達(dá)到底層樓面時，洪水對房屋柱體的沖擊壓力沿高度方向呈倒“C”型分布，而且水深對壓力分布的最大值和最小值影響不大。當(dāng)洪水漫過樓面時，柱體在水面附近受到的壓力較大，柱體其他位置的壓力較小，分布較均勻。

（3）“過流堰”防護(hù)措施（門前低矮斜坡）可以顯著降低結(jié)構(gòu)穩(wěn)定繞流階段的穩(wěn)定位移，但應(yīng)力降低不明顯。短柱擾流防護(hù)措施能夠有效減少結(jié)構(gòu)位移及柱底拉應(yīng)力，降低混凝土開裂風(fēng)險。且短柱擾流防護(hù)造價低，施工簡單，不需作專項設(shè)計，可在山區(qū)低層房屋抗洪設(shè)計中做適當(dāng)推廣。

本文主要關(guān)注不同洪水深度對房屋結(jié)構(gòu)的作用。其他參數(shù)，如洪水流速、含沙量等也可能影響結(jié)構(gòu)的響應(yīng)，這些值得后續(xù)做進(jìn)一步的研究。

參考文獻(xiàn)

［1］ 中華人民共和國水利部.中國水旱災(zāi)害防御公報2020［R］.2021

［2］ 潘佳佳，曹志先，王協(xié)康，等.暴雨山洪水動力學(xué)模型及其簡化模型的比較研究［J］.四川大學(xué)學(xué)報（工程科學(xué)版），2012，44（S1）：77-82.

［3］ 肖詩云，王曉慶.洪水演進(jìn)模型及沖擊荷載數(shù)值分析［J］.工程力學(xué)，2010，27（9）：35-40.

［4］ Kundzewicz，Kanae，Seneviratne，et al. Flood risk and climate change： global and regional perspectives［J］. Hydrological Sciences Journal，2014，59（1）.

［5］ Alessio Radice，Laura Longoni，Monica Papini，et al. Generation of a Design Flood-Event Scenario for a Mountain River with Intense Sediment Transport［J］. Water，2016，8（12）.

［6］ 肖詩云，楊留娟，岳斌，等.不同開洞率鄉(xiāng)村建筑洪水水流荷載試驗［J］.沈陽建筑大學(xué)學(xué)報（自然科學(xué)版），2010，26（6）：1039-1045.

［7］ 于文，葛學(xué)禮，朱立新，等.山區(qū)鄉(xiāng)村房屋模型水流作用試驗［J］.北京交通大學(xué)學(xué)報，2011，35（1）：24-27.

［8］ 肖詩云，楊留娟，岳斌，等.山區(qū)鄉(xiāng)村房屋模型洪水沖擊試驗研究［J］.防災(zāi)減災(zāi)工程學(xué)報，2010，30（3）：235-240.

［9］ 鄭涵午，黃爾，路信，等.山區(qū)洪水對下游房屋沖擊特性的試驗研究［J］.中國農(nóng)村水利水電，2020（7）：41-45.

［10］ 楊留娟. 山區(qū)鄉(xiāng)村建筑洪水作用荷載試驗研究［D］.大連：大連理工大學(xué)，2010.

［11］ 楊哲豪，吳鋼鋒，張科鋒，等.潰壩洪水對構(gòu)筑物沖擊荷載的數(shù)值模擬［J］.長江科學(xué)院院報，2020，37（3）：45-50.

［12］ 范力陽，朱純，楊冰清，等.移動式防洪墻系統(tǒng)的三維有限元模擬［J］.水科學(xué)與工程技術(shù)，2021（1）：1-3.

［13］ 王友彪，姚昌榮，劉賽智，等.泥石流對橋墩沖擊力的試驗研究［J］.巖土力學(xué)，2019，40（2）：616-623.

［14］ 劉杰.潰壩洪水對并聯(lián)雙橋墩沖擊過程數(shù)值模擬［J］.人民珠江，2016，37（8）：6-9.

［15］ 張吉保.洪水漂浮物沖擊作用下橋梁結(jié)構(gòu)動力響應(yīng)研究［D］.哈爾濱：哈爾濱工業(yè)大學(xué)，2016.

［16］ 向琪芪，李亞東，魏凱，等.橋梁基礎(chǔ)沖刷研究綜述［J］.西南交通大學(xué)學(xué)報，2019，54（2）：235-248.

［17］ 郭良，丁留謙，孫東亞，等.中國山洪災(zāi)害防御關(guān)鍵技術(shù)［J］.水利學(xué)報，2018，49（9）：1123-1136.

［18］ 張志彤.山洪災(zāi)害防治措施與成效［J］.水利水電技術(shù)，2016，47（1）：1-5+11.

［19］ 侯極，劉興年，蔣北寒，等.山洪攜帶泥沙引發(fā)的山區(qū)大比降河流水深變化規(guī)律研究［J］.水利學(xué)報，2012，43（S2）：48-53.