• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    軟基水閘底板脫空動力學(xué)反演模型構(gòu)建與試驗驗證

    2021-01-14 05:14:32李火坤魏博文黃錦林
    農(nóng)業(yè)工程學(xué)報 2020年21期
    關(guān)鍵詞:水閘反演底板

    李火坤,余 杰,王 剛,魏博文,黃 偉,黃錦林

    軟基水閘底板脫空動力學(xué)反演模型構(gòu)建與試驗驗證

    李火坤1,余 杰1,王 剛1,魏博文1※,黃 偉1,黃錦林2,3

    (1. 南昌大學(xué)建筑工程學(xué)院,南昌 330031;2. 廣東省水利水電科學(xué)研究院,廣州 510635;3. 河口水利技術(shù)國家地方聯(lián)合工程實驗室,廣州 510635)

    軟基水閘極易因地基不均勻沉降、滲透變形等發(fā)生底板脫空現(xiàn)象,傳統(tǒng)的無損檢測方法難以在水下探測底板脫空,且不具實時性,其應(yīng)用還存在一定局限性。該研究基于軟基水閘室內(nèi)物理模型,提出了一種軟基水閘底板脫空動力學(xué)反演方法。首先,基于多參數(shù)變量的底板脫空范圍數(shù)學(xué)模型,改進(jìn)了反映軟基水閘底板脫空參數(shù)與水閘動力學(xué)參數(shù)(模態(tài)參數(shù))之間非線性映射關(guān)系的數(shù)學(xué)代理模型;其次,以水閘多測點多階頻率和振型變化率組合作為動力學(xué)敏感特征量,建立了軟基水閘底板脫空參數(shù)反演的目標(biāo)函數(shù),并基于遺傳算法對目標(biāo)函數(shù)進(jìn)行求解。最后,建立了軟基水閘室內(nèi)物理模型,并在模型中設(shè)置3種不同的底板脫空工況,采用軟基水閘底板脫空動力學(xué)反演方法對軟基水閘底板的脫空進(jìn)行反演識別與模型驗證。結(jié)果表明:基于動力學(xué)參數(shù)反演識別的水閘底板脫空區(qū)域趨勢和脫空面積與模型實際脫空情況吻合較好,3種工況下水閘底板脫空面積反演結(jié)果與模型實際脫空面積的相對誤差分別為7.47%、6.78%、6.90%,驗證了該方法的可靠性,可望為軟基水閘實際工程底板脫空隱患檢測提供一種新的思路。

    模型;試驗;反演;軟基水閘;底板脫空;動力學(xué)參數(shù)

    0 引 言

    水閘是水利工程中應(yīng)用廣泛而重要的控制性建筑物,中國建成流量超5 m3/s以上的水閘達(dá)10萬余座,在調(diào)控水資源、防汛抗洪、蓄水灌溉等方面發(fā)揮著重要作用,對推動國內(nèi)的農(nóng)業(yè)發(fā)展具有重要影響[1-5]。水閘經(jīng)過長久運(yùn)行,極易產(chǎn)生各類影響水閘安全運(yùn)行的病患,特別是位于平原地區(qū)或沿海地區(qū)的水閘,大部分建于軟基上,其存在軟土厚、含水量高、壓縮性高、承載力低的特點[6-7],閘基底部在水流作用下會出現(xiàn)不均勻沉降、滲透變形等問題[8-9],從而造成水閘底板脫空而導(dǎo)致水閘坍塌或閘底斷裂等事故[10-12],嚴(yán)重影響水閘正常功能的發(fā)揮。目前,在檢測水閘底板病患的方法中以探地雷達(dá)法應(yīng)用最為廣泛,戴呈祥等[13-14]利用探地雷達(dá)技術(shù)探測了軟基水閘閘基沖刷剝蝕程度、淘空區(qū)的位置和范圍,對其不同類型閘基隱患的雷達(dá)圖形特征進(jìn)行了分析,并指出雷達(dá)圖像的解譯和特征分析是判別閘基隱患類型和部位的關(guān)鍵;楊松華[15]利用探地雷達(dá)檢測技術(shù)對廈門馬鑾水閘閘基、護(hù)坦、海漫底板下閘基等工程隱患的位置及范圍進(jìn)行了探測,并根據(jù)地質(zhì)鉆探成果對檢測結(jié)果進(jìn)行了驗證,同時指出該方法易受外界干擾,應(yīng)用之前需對雷達(dá)模型進(jìn)行驗證;徐云乾等[16]應(yīng)用探地雷達(dá)技術(shù)檢測佛山市沙口水利樞紐引水閘底板的滲流隱患,給出了滲流隱患區(qū)的探地雷達(dá)圖像特征和隱患的定性識別方法;吳俊姿等[17-18]將探地雷達(dá)技術(shù)應(yīng)用于水閘脫空檢測之中,探討了探地雷達(dá)的天線選擇及測線布置,分析了水閘脫空狀況與雷達(dá)圖像之間的關(guān)系;劉金濤等[19]應(yīng)用探地雷達(dá)對水閘的底板進(jìn)行了檢測,有效地探測了水閘底板的厚度、滲漏通道情況及地基地層分布情況等。從上述探地雷達(dá)法的應(yīng)用特點來看,要求檢測設(shè)備(如雷達(dá)或聲發(fā)射探頭)能夠觸及被測結(jié)構(gòu),且前期需建立反射信號與各類結(jié)構(gòu)或地層物理特性之間的先驗?zāi)P?,以此修正實際探測圖像的解譯成果,因此在實施之前,還需進(jìn)行模型驗證。對于軟基水閘底板而言,其在正常運(yùn)行時都是位于水下,傳統(tǒng)的檢測設(shè)備難以觸及水閘底板,若要實現(xiàn)運(yùn)行期水閘底板脫空的實時檢測和診斷,還需要研究新的診斷方法。軟基水閘底板脫空會造成水閘整體約束邊界條件改變,從而使水閘整體模態(tài)參數(shù)發(fā)生變化,通過建立模態(tài)參數(shù)變化與水閘底板脫空的非線性映射關(guān)系可實現(xiàn)軟基水閘底板脫空診斷,基于模態(tài)參數(shù)的軟基水閘底板脫空診斷方法能夠很好地實現(xiàn)實時檢測。國內(nèi)外關(guān)于軟基水閘底板脫空動力學(xué)診斷方法的研究不多,陳鸚[20]在國內(nèi)外率先提出了基于水閘底板動力特性的變化的軟基水閘底板的脫空檢測思想,將水閘底板以下的軟基視為彈性地基,采用柔度和黏性阻尼系數(shù)作為底板脫空診斷的動力學(xué)敏感特征量,應(yīng)用脈沖錘擊法進(jìn)行室內(nèi)模態(tài)試驗以驗證該方法的可行性,并進(jìn)行了實際工程檢驗,結(jié)果表明該方法是可行的,但該方法需在特定的條件下進(jìn)行,即當(dāng)水閘底板處于無水時傳感器和激勵可觸及底板,對于正常運(yùn)行的水閘(水閘底板過水)的檢測還存在一定局限性。黃錦林等[21]基于數(shù)值模型,采用水閘頻率和振型變化率作為動力學(xué)敏感特征量,對軟基水閘的底板脫空進(jìn)行識別,進(jìn)一步驗證了基于動力學(xué)參數(shù)的底板脫空診斷的可行性,但該方法僅在數(shù)值模型上開展,在結(jié)合水閘實際物理模型或原型工程開展振動測試與模態(tài)辨識、響應(yīng)面數(shù)學(xué)代理模型的精準(zhǔn)表達(dá)等方面還需進(jìn)一步改進(jìn)。

    本文基于軟基水閘室內(nèi)物理模型,提出一種軟基水閘底板脫空動力學(xué)反演方法,基于多參數(shù)變量的底板脫空范圍數(shù)學(xué)模型,改進(jìn)反映軟基水閘底板脫空參數(shù)與水閘動力學(xué)參數(shù)(模態(tài)參數(shù))之間非線性映射關(guān)系的數(shù)學(xué)代理模型;以水閘多測點多階頻率和振型變化率組合作動力學(xué)敏感特征量,建立水閘底板脫空反演的目標(biāo)函數(shù),并基于遺傳算法對目標(biāo)函數(shù)進(jìn)行求解,以實現(xiàn)底板脫空的動力學(xué)反演。

    1 軟基水閘底板脫空動力學(xué)反演理論方法

    1.1脫空模式及脫空范圍的參數(shù)化數(shù)學(xué)模型

    以文獻(xiàn)[21]所提出的軟基水閘底板脫空模式分類為參考,本文以相鄰側(cè)型和相對側(cè)型脫空為例開展模型試驗研究,將水閘底板相鄰側(cè)型和相對側(cè)型脫空范圍分別以5個和10個脫空參數(shù)變量進(jìn)行數(shù)學(xué)描述,并確定基于多參數(shù)的相鄰側(cè)型和相對側(cè)型底板脫空范圍數(shù)學(xué)模型,其數(shù)學(xué)模型表達(dá)可參考文獻(xiàn)[21]。

    1.2 軟基水閘底板脫空反演的數(shù)學(xué)代理模型

    本文軟基水閘底板脫空的反演方法屬于模型修正方法[22-25]的一類,其中響應(yīng)面法是一種能夠通過構(gòu)造具有顯性的多項式來表達(dá)隱式的非線性函數(shù)的可靠手段,為準(zhǔn)確表達(dá)水閘模態(tài)參數(shù)與底板脫空參數(shù)之間的非線性映射關(guān)系,構(gòu)建響應(yīng)面數(shù)學(xué)代理模型并代替有限元模型,本文以底板脫空參數(shù)作為模型的自變量,選擇水閘模態(tài)參數(shù)(頻率、振型)作為因變量,采用多項式函數(shù)來表達(dá)模態(tài)參數(shù)與脫空參數(shù)之間非線性映射關(guān)系,其非線性映射關(guān)系可表示為

    式中()、()分別表示水閘的固有頻率、節(jié)點振型值;表示脫空參數(shù)值;表示脫空參數(shù)個數(shù)。可將方程式(1)和(2)用Taylor級數(shù)公式展開,在準(zhǔn)確表達(dá)水閘脫空參數(shù)和模態(tài)參數(shù)之間的非線性函數(shù)關(guān)系下充分考慮計算效率,本文采用三階多項式響應(yīng)面模型,相較于文獻(xiàn)[21]中的三階多項式響應(yīng)面模型,增加了考慮脫空參數(shù)之間相關(guān)性的二次交叉項,使得該響應(yīng)面模型能夠更加準(zhǔn)確地描述水閘底板脫空參數(shù)與模態(tài)參數(shù)之間的非線性數(shù)學(xué)關(guān)系,具體形式為

    響應(yīng)面數(shù)學(xué)代理模型確定基本方法為:采用拉丁超立方抽樣法[26-27]在脫空參數(shù)取值范圍內(nèi)來生成合理的脫空參數(shù)樣本d(2 000組以上),以此作為軟基水閘基準(zhǔn)有限元模型模態(tài)計算的脫空參數(shù)輸入,計算不同脫空組合下的模態(tài)參數(shù)并得到對應(yīng)的固有頻率和節(jié)點振型值,采用多元回歸分析方法對上述方程進(jìn)行求解并擬合待定系數(shù)和,從而確定數(shù)學(xué)代理模型。

    在擬合響應(yīng)面方程時,響應(yīng)面精度評價是檢驗響應(yīng)面模型是否符合要求的重要指標(biāo),通過響應(yīng)面模型與有限元模型之間的相對誤差來反映響應(yīng)面模型的回歸精度,以此衡量擬合的響應(yīng)面模型能否有效地替代有限元模型。本文通過有限元模型計算的及響應(yīng)面模型輸出的固有頻率f和節(jié)點振型值φ來計算響應(yīng)面模型的精度:

    式中表示響應(yīng)面模型的精度,響應(yīng)面模型精度控制在5‰以內(nèi)[28];RS表示響應(yīng)面模型模態(tài)參數(shù)的輸出值,F(xiàn)EM表示有限元模型模態(tài)參數(shù)的計算值。

    1.3 軟基水閘底板脫空反演的目標(biāo)函數(shù)

    基于數(shù)學(xué)代理模型和動力學(xué)敏感特征量,構(gòu)建底板脫空參數(shù)反演的最優(yōu)化數(shù)學(xué)模型,通過建立目標(biāo)函數(shù)*來反映水閘振動測試識別的動力學(xué)敏感特征量與代理模型(響應(yīng)面模型)計算得到動力學(xué)敏感特征量的相對偏差,將脫空參數(shù)的反演表達(dá)為最優(yōu)化求解問題。以軟基水閘多測點多階頻率和振型相對變化率組合作為動力學(xué)敏感特征量,構(gòu)造的目標(biāo)函數(shù)為:

    1.4 軟基水閘底板脫空反演基本流程

    綜合以上,本文軟基水閘底板脫空動力學(xué)反演流程如圖1所示。

    圖1 軟基水閘底板脫空動力學(xué)反演流程

    2 模型試驗驗證

    2.1 軟基水閘物理模型的制作

    為驗證本文方法的可靠性,以江西信江水閘工程為背景實例,制作一單孔軟基水閘室內(nèi)物理模型,模型比尺為1:10,模型基本尺寸為:水閘底板順?biāo)鞣较蜷L1.44 m,橫河向?qū)?.36 m,底板前后設(shè)有齒槽,水閘底板厚0.16 m,閘墩高1.6 m、厚0.16 m,水閘前、后部工作橋?qū)挿謩e為0.32和0.4 m,厚0.04 m,水閘模型示意圖如圖2a所示。為保證在進(jìn)行動力測試過程中振動波傳播介質(zhì)與原型水閘的一致性,水閘結(jié)構(gòu)模型采用與原型一致的鋼筋混凝土材料,軟基的材料模擬也與原型水閘地基材料基本保持一致,即采用細(xì)沙、礫石與黏土按一定界層填充壓實來模擬軟土地基,其中,軟土地基最底部細(xì)沙與礫石層厚0.3 m,往上為黏土層厚0.08 m,最上部為水閘底板凹槽處黏土層厚0.16 m。試驗制作完成的軟基水閘物理模型如圖2b所示。

    2.2 軟基水閘模型底板脫空的模擬及脫空工況設(shè)置

    基于水閘物理模型,采用人工掏空的方式對水閘模型底板進(jìn)行相鄰側(cè)型脫空和相對側(cè)型脫空的隨機(jī)模擬,并測量脫空縱深,如圖3所示;在底板脫空范圍的數(shù)學(xué)模型中,采用5個參數(shù)描述相鄰側(cè)脫空(即=5),采用10個參數(shù)描述相對側(cè)型脫空(即=10)。脫空模擬設(shè)置共分為3種工況:工況一和工況二均為水閘底板上游相鄰側(cè)型脫空,工況三為水閘底板相對側(cè)型脫空,不同脫空工況下的脫空參數(shù)(即脫空縱深)如表1所示。

    1.軟基 2.水閘底板 3.閘墩 4.工作橋

    注:d1~d5為水閘底板上游脫空參數(shù),d6~d10為水閘底板下游脫空參數(shù)。

    表1 水閘物理模型底板脫空工況及脫空參數(shù)

    2.3 軟基水閘模型振動測試及模態(tài)參數(shù)辨識

    2.3.1 軟基水閘模型振動測試

    水閘模型振動測試傳感器布置如圖4a所示,共布置20個水平向振動傳感器,振動傳感器采用BY-S07型高精度振動傳感器。傳感器布置方式如下:在水閘模型左右兩側(cè)沿閘墩水流方向依次對稱布置8個振動傳感器(H1~H4為右閘墩測點、H11~H14為左閘墩測點),在左右閘墩沿高度方向依次對稱布置12個振動傳感器(H5~H10為右閘墩測點、H15~H20為左閘墩測點)。水閘振動測試的測振物理量為速度,采樣頻率為500 Hz,采用DASP智能數(shù)據(jù)采集系統(tǒng)[30]進(jìn)行數(shù)據(jù)采集與預(yù)處理;本次模型振動測試以人工力棒瞬時脈沖激勵的方式對水閘進(jìn)行激勵并使其產(chǎn)生振動,如圖4b所示,典型工況(工況三)代表性測點振動響應(yīng)時程線及功率譜如圖5所示;根據(jù)所采集的各工況下的振動響應(yīng)時程進(jìn)行模態(tài)辨識即可獲取各脫空工況下的水閘模態(tài)參數(shù)。

    注:H1~H10為布置在水閘右側(cè)閘墩測點,H11~H20為布置在水閘左側(cè)閘墩測點。下同。

    2.3.2 各脫空工況下的水閘模型模態(tài)辨識結(jié)果

    根據(jù)水閘模型各底板脫空工況下測點的振動響應(yīng),采用基于奇異熵定階的隨機(jī)子空間法及穩(wěn)定圖法對水閘模型進(jìn)行模態(tài)參數(shù)識別。在人工脈沖激勵下,水閘模型共被激發(fā)出4階模態(tài),工況一識別得到的前4階頻率分別為:20.26、24.81、53.31和69.72 Hz;工況二識別得到前4階頻率分別為:20.24、24.80、52.20和69.58 Hz;工況三識別得到前4階頻率分別為:20.23、24.65、48.42和69.77 Hz。分析不同脫空工況下水閘模態(tài)頻率的變化趨勢可知:隨著水閘脫空區(qū)域的逐漸增大,水閘整體頻率逐漸變小,且對于該水閘模型,底板發(fā)生脫空對第三階頻率影響最大,該階頻率變化對底板脫空敏感性最強(qiáng),主要是第三階的振動在底板附近較大,而其他階次的振動在底板附近極小,當(dāng)水閘底板脫空使底板約束邊界條件發(fā)生改變時,第三階的振動變化較其他階次更大,即第三階頻率變化對底板脫空敏感性最強(qiáng)。

    圖5 典型測點振動響應(yīng)時程線和歸一化功率譜

    2.4 響應(yīng)面數(shù)學(xué)代理模型的建立

    本文響應(yīng)面數(shù)學(xué)代理模型的建立以有限元數(shù)值模型為基準(zhǔn),根據(jù)上述室內(nèi)物理模型對比建立相應(yīng)的軟基水閘ANSYS有限元模型。水閘有限元數(shù)值模型采用六面實體單元進(jìn)行網(wǎng)格劃分(圖6),共有20 036個單元與24 455個節(jié)點。在有限元數(shù)值模型中,采用溫克爾地基模型模擬軟基,定義Surface154單元附著于水閘底板(即底板與地基接觸面)并通過設(shè)置彈性地基剛度(Elastic foundation stiffness,EFS)值來模擬溫克爾地基。針對水閘數(shù)值模型的脫空部分,可將EFS值置為零來實現(xiàn)底板脫空的模擬。有限元模型中的材料參數(shù)通過與物理模型進(jìn)行對比修正獲得,如表2所示。

    圖6 軟基水閘有限元網(wǎng)格劃分

    表2 有限元模型材料參數(shù)

    水閘底板相鄰側(cè)型脫空和相對側(cè)型脫空分別采用5個、10個脫空控制參數(shù)對水閘底板脫空區(qū)域進(jìn)行數(shù)學(xué)描述。采用拉丁超立方抽樣法來隨機(jī)生成2 000組脫空控制參數(shù)樣本,其中工況一和工況二中5個脫空控制參數(shù)取樣范圍為0~0.65 m;工況三中上游側(cè)和下游側(cè)的10個脫空控制參數(shù)取樣范圍分別為0~0.65和0~0.4 m。將隨機(jī)生成的2 000組脫空控制參數(shù)以程序設(shè)計的方式輸入所建立的軟基水閘基準(zhǔn)有限元模型,并進(jìn)行有限元模態(tài)計算,提取各組對應(yīng)脫空狀態(tài)下的前四階頻率與振型,根據(jù)式(3)、(4)擬合響應(yīng)面方程的待定系數(shù)和,并將擬合后的響應(yīng)面模型替代有限元模型。經(jīng)擬合,3個工況下的響應(yīng)面方程均滿足精度要求(均在5‰以內(nèi)),限于篇幅,本文以工況三為典型代表給出頻率和典型測點振型響應(yīng)面精度,如圖7和8所示。

    圖7 工況三前4階頻率響應(yīng)面精度

    2.5 軟基水閘底板脫空反演結(jié)果與驗證

    根據(jù)圖1所示反演流程,得到水閘底板不同脫空工況下的脫空控制參數(shù)與模型脫空實際值對比如表3所示,反演識別脫空區(qū)域面積與實際脫空區(qū)域面積相對誤差如表4所示,不同脫空工況下脫空面積與趨勢反演識別結(jié)果與實際脫空結(jié)果對比如圖9所示。從反演識別的結(jié)果來看,本文所提出的方法所反演識別的底板脫空區(qū)域趨勢和脫空面積與模型實際脫空情況吻合較好,3種工況下水閘底板脫空面積反演結(jié)果與模型實際脫空面積的相對誤差分別為7.47%、6.78%、6.90%,相對誤差較小,驗證了本方法的可靠性。

    圖8 工況三典型測點振型響應(yīng)面精度

    表3 不同工況下脫空參數(shù)反演結(jié)果

    圖9 不同工況下脫空面積和趨勢識別結(jié)果與實際脫空值對比

    3 結(jié)論與建議

    1)改進(jìn)了反映軟基水閘底板脫空參數(shù)與水閘動力學(xué)參數(shù)之間非線性映射關(guān)系的響應(yīng)面數(shù)學(xué)代理模型,以水閘多測點多階頻率和振型變化率組合作為動力學(xué)敏感特征量,建立了水閘底板脫空反演的目標(biāo)函數(shù),并給出了脫空參數(shù)的反演流程。

    2)建立了軟基水閘室內(nèi)模型,對水閘底板相鄰側(cè)脫空和相對側(cè)脫空進(jìn)行了物理模擬,采用改進(jìn)的軟基水閘底板脫空動力學(xué)反演方法進(jìn)行了脫空參數(shù)的反演識別并與實際脫空情況進(jìn)行了對比驗證,結(jié)果表明,反演識別的水閘底板脫空區(qū)域趨勢和脫空面積與模型實際脫空情況吻合較好,3種工況下水閘底板脫空面積反演結(jié)果與模型實際脫空面積的相對誤差分別為7.47%、6.78%、6.90%,驗證了該方法的可靠性。

    頻率和各測點振型對水閘底板脫空參數(shù)的敏感性不同會對目標(biāo)函數(shù)中頻率和振型變化率的權(quán)重系數(shù)產(chǎn)生影響,如何從靈敏度角度分析各模態(tài)參數(shù)對底板脫空敏感性并確定其權(quán)重系數(shù),進(jìn)一步提高反演精度,是值得深入研究的。

    [1] 中華人民共和國水利部. 全國水利發(fā)展統(tǒng)計公報[M]. 北京:中國水利水電出版社, 2018.

    [2] 倪文進(jìn). 中國農(nóng)村水利發(fā)展?fàn)顩r與科技需求[J]. 農(nóng)業(yè)工程學(xué)報,2010,26(3):1-8.

    Ni Wenjin. Development and technology requirement of China rural water conservancy[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(3): 1-8. (in Chinese with English abstract)

    [3] 包加桐, 錢江, 張煒, 等. 基于多通道數(shù)據(jù)流在線相關(guān)分析及聚類的閘站工程安全監(jiān)測[J]. 農(nóng)業(yè)工程學(xué)報, 2019, 35(3):109-116.

    Bao Jiatong, Qian Jiang, Zhang Wei, et al. Safety monitoring of sluice-pump station project based on online correlation analysis and clustering of multichannel data streams[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(3): 109-116. (in Chinese with English abstract)

    [4] 汪恩良. 寒區(qū)水閘計算機(jī)輔助設(shè)計軟件系統(tǒng)開發(fā)研究[J]. 農(nóng)業(yè)工程學(xué)報,1999,15(2):239-239.

    Wang Enliang. Research and development of computer aided design software system for sluices in cold area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 1999, 15(2): 239-239. (in Chinese with English abstract)

    [5] 劉昉,谷欣玉,李文勝,等. 動水關(guān)閉的平面事故閘門體型優(yōu)化試驗研究[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(12):142-149.

    Liu Fang, Gu Xinyu, Li Wensheng, et al. Structural optimization of emergency plate gate for closure in moving water[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(12): 142-149. (in Chinese with English abstract)

    [6] 梁民陽,吳興龍. 浙東海塘上水閘病害成因分析及對策[J]. 中國農(nóng)村水利水電,2006(7):107-108.

    Liang Mingyang, Wu Xinglong. Cause and countermeasures of sluice disease in east sea of Zhejiang[J].China Rural Water and Hydropower, 2006(7): 107-108. (in Chinese with English abstract)

    [7] 張帆. 珠三角地區(qū)水閘軟土基礎(chǔ)沉降處理案例分析[J]. 西北水電,2007(2):50-51, 73.

    Zhang Fan. Case analysis on soft foundation settlement of sluice gates in delta region around Pear River[J].Northwest Hydropower, 2007(2): 50-51, 73. (in Chinese with English abstract)

    [8] 何金平,曹旭梅,李紹文,等. 基于安全監(jiān)測的水閘健康診斷體系研究[J]. 水利水運(yùn)工程學(xué)報,2018(5):1-7.

    He Jinping, Cao Xumei, Li Shaowen, et al. Analysis of sluice health diagnosis system based on safety monitoring[J]. Hydro-Science and Engineering, 2018(5): 1-7. (in Chinese with English abstract)

    [9] 鐘俊紀(jì). 水閘軟土地基沉降預(yù)測的主要方法[J]. 理論與現(xiàn)代化,2005(S1):131-132.

    Zhong Junji. Main methods of settlement prediction of sluice on soft foundation[J].Theory and Modernizatio, 2005(S1): 131-132. (in Chinese with English abstract)

    [10] 陳富強(qiáng),楊光華,黃致興,等. 某水閘樁基底板脫空案例分析及設(shè)計反思[J]. 建筑科學(xué),2015,31(S2):197-203.

    Chen Fuqiang, Yang Guanghua, Huang Zhixing, et al. Case analysis and design reflection of a sluice floor void[J]. Building Science, 2015, 31(S2): 197-203. (in Chinese with English abstract)

    [11] 陳仲策. 水閘基礎(chǔ)處理方法探討[J]. 中國防汛抗旱,2010(4):59-60,71.

    Chen Zhongce. Discussion on foundation treatment method of sluice[J]. China Flood & Drought Management, 2010(4): 59-60, 71. (in Chinese with English abstract)

    [12] 關(guān)錦榮. 沿海地區(qū)水閘破壞形式與加固對策[J]. 水利科技,2009(3):59-61.

    Guan Jinrong. Damaged forms of sluices in coastal areas and the strengthening countermeasures[J].Hydraulic Science and Technology, 2009(3): 59-61. (in Chinese with English abstract)

    [13] 戴呈祥,王士恩. 水閘閘基隱患探測雷達(dá)圖像特征分析[J]. 地球物理學(xué)進(jìn)展,2003,18(3):429-433.

    Dai Chenxiang, Wang Shien. Feature analysis on radar pictures for detecting hidden danger of floodgate base[J]. Progress in Geophysics, 2003, 18(3): 429-433. (in Chinese with English abstract)

    [14] 戴呈祥,王士恩. 水閘閘基隱患類型特征分析[J]. 工程地球物理學(xué)報,2004,1(4):353-357.

    Dai Chengxiang, Wang Shien. The character analysis of hidden danger types of the floodgate base[J]. Chinese Journal of Engineering Geophysics, 2004, 1(4): 353-357. (in Chinese with English abstract)

    [15] 楊松華. 地質(zhì)雷達(dá)檢測在廈門馬鑾水閘安全鑒定中的應(yīng)用[J]. 福建建筑,2011,156(6):73-75.

    Yang Songhua. The application of ground penetrating radar examination in Xiamen maluan sluice safe appraisal[J]. Fujian Architecture & Construction, 2011, 156(6): 73-75. (in Chinese with English abstract)

    [16] 徐云乾,楊文濱,袁明道,等. 水閘底板滲流隱患的探地雷達(dá)檢測[J]. 無損檢測,2017,39(9):83-86.

    Xu Yunqiang, Yang Wenbin, Yuan Mingdao, et al. The detection of seepage hazard in sluice gate bottom plate using GPR[J]. Nondestructive Testing, 2017, 39(9): 83-86. (in Chinese with English abstract)

    [17] 吳俊姿. 探地雷達(dá)在水閘脫空檢測中的應(yīng)用[J]. 水利技術(shù)監(jiān)督,2018(3):233-236.

    Wu Junzi. Application of ground penetrating radar in the detection of sluice voids[J]. Technical Supervision in Water Resources, 2018(3): 233-236. (in Chinese with English abstract)

    [18] 安鐸,陸新宇,李經(jīng)緯. 基于探地雷達(dá)技術(shù)的水閘脫空檢測方法研究[J]. 人民長江,2014,45(2):202-205.

    An Duo, Lu xinyu, Li jingwei. Research on detection method of sluice void based on ground penetrating radar technology[J].Yangtze River, 2014, 45(2): 202-205. (in Chinese with English abstract)

    [19] 劉金濤,張旭輝,潘展釗,等. 基于探地雷達(dá)的水閘水下結(jié)構(gòu)探測[J]. 無損檢測,2019,41(5):61-63.

    Liu Jingtao, Zhang Xuhui, Pan Zhanzhao, et al. Detection of underwater structures of a sluice based on GPR[J]. Nondestructive Testing, 2019, 41(5): 61-63. (in Chinese with English abstract)

    [20] 陳鸚. 無損動態(tài)檢測閘基隱患[J]. 大壩與安全,1994(1):40-46.

    Chen Ying. Non-destructive dynamic detection of hidden dangers of gate base[J]. Dam & Safety, 1994(1): 40-46. (in Chinese with English abstract)

    [21] 黃錦林,李火坤,鄧冰梅. 基于響應(yīng)面理論的閘基底板脫空區(qū)域識別方法[J]. 中國水利水電科學(xué)研究院學(xué)報,2018,16(4):249-256.

    Huang Jinlin, Li Huokun, Deng Bingmei. Identification methods of sluice floor pavement area based on the theory of response surface[J]. Journal of China Institute of Water Resources and Hydropower Research, 2018, 16(4): 249-256. (in Chinese with English abstract)

    [22] Girardi M, Padovani C, Pellegrini D, et al. Finite element model updating for structural applications[J]. Journal of Computational and Applied Mathematics, 2020(370): 112675.

    [23] Hofmeister B, Bruns M, Rolfes R. Finite element model updating using deterministic optimisation: A global pattern search approach[J]. Engineering Structures, 2019(195): 373-381.

    [24] Kang J, Zhang X, Cao H, et al. Research on multi-alternatives problem of finite element model updating based on IAFSA and kriging model[J]. Sensors, 2020, 20(15): 4274.

    [25] 費慶國,韓曉林,蘇鶴玲. 響應(yīng)面有限元模型修正的實現(xiàn)與應(yīng)用[J]. 振動、測試與診斷,2010,30(2):132-134.

    Fei Qingguo, Han Xiaolin, Su Heling. Response surface based finite element model updating and its application[J]. Journal of Vibration, Measurement & Diagnosis, 2010, 30(2): 132-134. (in Chinese with English abstract)

    [26] 方磊. 基于LHS抽樣的不確定性分析方法在概率安全評價中的應(yīng)用研究[D]. 合肥:中國科學(xué)技術(shù)大學(xué),2015.

    Fang Lei. Application Research of Uncertainty Analysis Method Based on LHS Sampling in Probabilistic Safety Valuation[D]. Hefei: University of Science and Technology of China, 2015. (in Chinese with English abstract)

    [27] Jirutitijaroen P, Singh C. Comparison of simulation methods for power system reliability indexes and their distributions[J]. IEEE Transactions on Power Systems, 2008, 23(2): 486-493.

    [28] Li H, Wang G, Wei B, et al. Dynamic inversion method for the material parameters of a high arch dam and its foundation[J]. Applied Mathematical Modelling, 2019(71): 60-76.

    [29] Li H, Zhong Y, Wei B, et al. Beat vibration mechanism of a sluice pier under high-speed flood discharge excitation [J]. Journal of Low Frequency Noise Vibration and Active Control, 2020, 39(1): 28-43.

    [30] 賀玉龍,張群,陳瑞,等. 成都地鐵誘發(fā)室內(nèi)結(jié)構(gòu)噪聲實測與分析[J]. 噪聲與振動控制,2020,40(2):202-206.

    He Yulong, Zhang Qun, Chen Rui, et al. Measurements and analyses of indoor ground-borne noise induced by Chengdu metro[J]. Noise and Vibration Control, 2020, 40(2): 202-206. (in Chinese with English abstract)

    Model construction of dynamic inversion and experimental verification for the void of sluice floor on the soft foundation

    Li Huokun1, Yu Jie1, Wang Gang1, Wei Bowen1※, Huang Wei1, Huang Jinlin2,3

    (1.330031,; 2.510635,; 3.510635,)

    As the uneven settlement and seepage deformation of the soft foundation, the floor void is prone to occurring in the sluices built on soft foundation. The traditional non-destructive testing methods are hardly difficult to detect the floor void of the sluices under water, which are non-real time, and their applications still have certain limitations. For the floor of sluice on the soft foundation, it has been under water during the operation period, and new diagnostic methods need to be studied to realize the real-time detection and diagnosis of the floor void of sluice on the soft foundation. The floor void of sluice on the soft foundation will cause the overall constraint boundary conditions of the sluice to change, thereby changing the overall modal parameters of the sluice. The diagnosis of the floor void of sluice on the soft foundation can be realized by establishing the nonlinear mapping relationship between the modal parameters and the floor void parameters of sluice. The method of diagnosing the floor void of sluice on the soft foundation based on modal parameters can realize real-time detection well. Therefore, a dynamic inversion method of the floor void of sluice on the soft foundation is proposed on the basis of the indoor physical model of the sluice on soft foundation in this research. First, based on mathematical model of multi-parameter variables of the floor void of sluice, a mathematical proxy model reflecting the nonlinear mapping relationship between the void parameters of the sluice floor and the dynamic parameters (modal parameters) of the sluice was improved. Second, the combination of the change rate of multi-order frequency and mode shape of multiple measurement points of the sluice was took as the dynamic sensitive characteristic quantity of the parameters inversion, and an objective function for the dynamic inversion of void parameters of sluice floor was established. Then the objective function is optimized through the genetic algorithm. Finally, a physical model of the indoor sluice on soft foundation had been made, and three different conditions of floor void of sluice on the soft foundation were set in the model. The prototype vibration test of the physical sluice on soft foundation on three conditions was carried out, and the corresponding modal parameters of sluice were identified based on the vibration response signal. And the dynamic inversion method of the floor void of sluice on the soft foundation proposed in the paper was used to invert the floor void of sluice and verified by the physical model. The results showed that the trend of floor void of sluice and the void area determined by the inversion method based on dynamic parameters of sluice were in good agreement with the actual condition of the physical model. With the three different conditions, the relative errors between the calculated void area and the actual void area are 7.47%, 6.78%, and 6.90%, respectively, which can confirm the reliability and rationality of the proposed method. The proposed method can provide a new idea for safety detection of the floor void of actual sluice on the soft foundation.

    models; test; inversion; sluice on soft foundation; floor void; dynamic parameters

    李火坤,余杰,王剛,等. 軟基水閘底板脫空動力學(xué)反演模型構(gòu)建與試驗驗證[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(21):145-153. doi:10.11975/j.issn.1002-6819.2020.21.018 http://www.tcsae.org

    Li Huokun, Yu Jie, Wang Gang, et al. Model construction of dynamic inversionand experimental verification for the void of sluice floor on the soft foundation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(21): 145-153. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.21.018 http://www.tcsae.org

    2020-06-30

    2020-10-30

    國家自然科學(xué)基金項目(52079061,51879126,51869011,51779115);江西省杰出青年基金項目(2018ACB21018,20192ACB21022);南昌大學(xué)研究生創(chuàng)新專項資金項目(CX2019117)

    李火坤,博士,教授,主要從事水工水力學(xué)與工程安全檢測方面研究。Email:lihuokun@ncu.edu.cn

    魏博文,博士,教授,主要從事水利工程安全監(jiān)測方向研究。Email:bwwei@ncu.edu.cn

    10.11975/j.issn.1002-6819.2020.21.018

    TV662

    A

    1002-6819(2020)-21-0145-09

    猜你喜歡
    水閘反演底板
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    2016年河南省已建成水閘數(shù)量
    全省已建成水閘數(shù)量
    軟土地基上超長底板節(jié)制閘沉降及應(yīng)力分析
    底板巷一巷兩用的回采工作面防火實踐
    底板隔水層破壞的力學(xué)模型及破壞判據(jù)
    河南省2014年已建成水閘數(shù)量
    疊前同步反演在港中油田的應(yīng)用
    大片电影免费在线观看免费| 欧美zozozo另类| 国产国拍精品亚洲av在线观看| 国产伦精品一区二区三区视频9| 日本午夜av视频| 精品人妻视频免费看| 精华霜和精华液先用哪个| 又大又黄又爽视频免费| 亚洲欧洲国产日韩| www.色视频.com| 免费电影在线观看免费观看| 亚洲怡红院男人天堂| 亚洲国产精品专区欧美| 国产精品国产三级国产av玫瑰| 午夜激情久久久久久久| 亚洲精品国产成人久久av| 黄色怎么调成土黄色| 亚洲va在线va天堂va国产| 国产亚洲精品久久久com| 亚洲欧美中文字幕日韩二区| 在线免费十八禁| 一级av片app| 免费黄频网站在线观看国产| 另类亚洲欧美激情| 亚洲无线观看免费| 亚洲欧美日韩东京热| 国产男女超爽视频在线观看| 精品人妻偷拍中文字幕| 日韩不卡一区二区三区视频在线| 精品午夜福利在线看| 久久影院123| 女人十人毛片免费观看3o分钟| 热99国产精品久久久久久7| 人妻系列 视频| 干丝袜人妻中文字幕| 久久精品国产鲁丝片午夜精品| 国产综合懂色| 在现免费观看毛片| 免费观看的影片在线观看| 国产亚洲一区二区精品| 亚洲精品一区蜜桃| 久热这里只有精品99| 日本熟妇午夜| 国产成人精品久久久久久| 亚洲av国产av综合av卡| 日韩不卡一区二区三区视频在线| 国产av不卡久久| 午夜亚洲福利在线播放| 亚洲欧美清纯卡通| 97精品久久久久久久久久精品| 国产91av在线免费观看| 久久久久久久国产电影| 亚洲av成人精品一区久久| 久久久欧美国产精品| 黄色配什么色好看| 亚洲av欧美aⅴ国产| 成人国产麻豆网| 国产男女内射视频| 69人妻影院| 国内精品宾馆在线| 免费在线观看成人毛片| 又爽又黄无遮挡网站| 久久久久久久亚洲中文字幕| 人人妻人人爽人人添夜夜欢视频 | 亚洲三级黄色毛片| 热re99久久精品国产66热6| 久久精品国产亚洲av天美| 国产在线男女| 国产精品一区二区在线观看99| 婷婷色麻豆天堂久久| 免费不卡的大黄色大毛片视频在线观看| 国产一区二区亚洲精品在线观看| 在线观看免费高清a一片| 亚洲欧美一区二区三区国产| 九九爱精品视频在线观看| av专区在线播放| 亚洲成人av在线免费| 久久99精品国语久久久| 国产成年人精品一区二区| 69人妻影院| 亚洲国产精品成人久久小说| 亚洲美女搞黄在线观看| 80岁老熟妇乱子伦牲交| 国产高清不卡午夜福利| av在线观看视频网站免费| 国产成人精品婷婷| 黄色怎么调成土黄色| 亚洲国产精品成人久久小说| 欧美高清性xxxxhd video| 可以在线观看毛片的网站| 久久久久久久久久久免费av| 日韩av不卡免费在线播放| 亚州av有码| 日本猛色少妇xxxxx猛交久久| 高清午夜精品一区二区三区| 综合色丁香网| 免费在线观看成人毛片| 亚洲av男天堂| 成年女人看的毛片在线观看| 丝袜喷水一区| 大香蕉97超碰在线| 爱豆传媒免费全集在线观看| 神马国产精品三级电影在线观看| 一个人看视频在线观看www免费| 嘟嘟电影网在线观看| 国产成人免费无遮挡视频| 亚洲熟女精品中文字幕| av天堂中文字幕网| 青春草视频在线免费观看| 内地一区二区视频在线| 国内精品宾馆在线| 又粗又硬又长又爽又黄的视频| 欧美区成人在线视频| 尤物成人国产欧美一区二区三区| 舔av片在线| freevideosex欧美| 老司机影院毛片| 99re6热这里在线精品视频| 另类亚洲欧美激情| 亚洲精品一区蜜桃| videossex国产| av在线播放精品| 日韩强制内射视频| 有码 亚洲区| 亚洲av福利一区| 久久久久久伊人网av| 九九久久精品国产亚洲av麻豆| 亚洲成人精品中文字幕电影| 久久久精品免费免费高清| 男女那种视频在线观看| 伊人久久国产一区二区| 欧美 日韩 精品 国产| 51国产日韩欧美| 国产久久久一区二区三区| 国产成人aa在线观看| 美女cb高潮喷水在线观看| 久久久久久久精品精品| 亚洲欧洲日产国产| 国产一区亚洲一区在线观看| 一级毛片电影观看| 18禁在线播放成人免费| 舔av片在线| 麻豆国产97在线/欧美| 国产精品国产三级国产av玫瑰| 亚洲av福利一区| 激情五月婷婷亚洲| 91在线精品国自产拍蜜月| 日韩成人av中文字幕在线观看| 午夜免费鲁丝| 免费观看的影片在线观看| 国产精品一二三区在线看| 国产亚洲一区二区精品| 超碰97精品在线观看| 国产精品久久久久久精品古装| av黄色大香蕉| 日韩一区二区视频免费看| kizo精华| 十八禁网站网址无遮挡 | 国产精品久久久久久av不卡| 国产成人免费观看mmmm| 一级毛片 在线播放| 日韩av免费高清视频| 国产视频首页在线观看| 亚洲欧美精品专区久久| 午夜日本视频在线| 亚洲伊人久久精品综合| 91aial.com中文字幕在线观看| 中文天堂在线官网| 深夜a级毛片| 欧美日韩视频精品一区| 亚洲美女视频黄频| 91久久精品国产一区二区三区| 18禁在线播放成人免费| 91久久精品国产一区二区成人| 91久久精品国产一区二区成人| 男人和女人高潮做爰伦理| 久久人人爽人人片av| 欧美成人午夜免费资源| 高清视频免费观看一区二区| 人体艺术视频欧美日本| 全区人妻精品视频| av天堂中文字幕网| 黄色配什么色好看| 内射极品少妇av片p| 三级男女做爰猛烈吃奶摸视频| av在线观看视频网站免费| 欧美bdsm另类| 国产精品久久久久久精品电影| 国产欧美另类精品又又久久亚洲欧美| 午夜激情久久久久久久| 免费大片18禁| 一区二区三区乱码不卡18| 99热全是精品| 亚洲在线观看片| 91狼人影院| 全区人妻精品视频| 免费看a级黄色片| 亚洲怡红院男人天堂| 91久久精品国产一区二区三区| 一区二区三区乱码不卡18| 51国产日韩欧美| 免费黄频网站在线观看国产| 尾随美女入室| 色网站视频免费| 亚洲精品自拍成人| 在线观看一区二区三区| 国产精品爽爽va在线观看网站| 18禁裸乳无遮挡免费网站照片| 丰满少妇做爰视频| 丝袜脚勾引网站| 又爽又黄a免费视频| 国产精品不卡视频一区二区| 九草在线视频观看| freevideosex欧美| 国产成年人精品一区二区| 校园人妻丝袜中文字幕| 亚洲精品日本国产第一区| 亚洲图色成人| 联通29元200g的流量卡| 欧美人与善性xxx| 欧美xxⅹ黑人| 日韩中字成人| 国产精品爽爽va在线观看网站| 热99国产精品久久久久久7| 99热6这里只有精品| 伊人久久国产一区二区| 精品少妇黑人巨大在线播放| 日韩av在线免费看完整版不卡| 久久久久精品久久久久真实原创| 日韩av免费高清视频| 亚洲va在线va天堂va国产| 国产精品一区www在线观看| 高清在线视频一区二区三区| 日日摸夜夜添夜夜添av毛片| 色吧在线观看| 久久久久久伊人网av| 狂野欧美激情性xxxx在线观看| 亚洲最大成人手机在线| 成人二区视频| 丝袜喷水一区| 国产精品一区二区三区四区免费观看| 欧美xxⅹ黑人| 国产成人免费观看mmmm| 久久久久九九精品影院| 内地一区二区视频在线| 免费看光身美女| 精品人妻视频免费看| 女人十人毛片免费观看3o分钟| 亚洲国产精品成人久久小说| 乱码一卡2卡4卡精品| 亚洲国产高清在线一区二区三| 大片电影免费在线观看免费| 亚洲精品影视一区二区三区av| 国产乱来视频区| 成人一区二区视频在线观看| 国产熟女欧美一区二区| 最新中文字幕久久久久| 国产精品女同一区二区软件| 天堂俺去俺来也www色官网| 可以在线观看毛片的网站| 精品熟女少妇av免费看| 亚洲精品乱久久久久久| 欧美成人a在线观看| 国产综合懂色| 最近中文字幕2019免费版| 日韩中字成人| 欧美激情在线99| videos熟女内射| 人妻制服诱惑在线中文字幕| 精华霜和精华液先用哪个| 国产成人免费无遮挡视频| 伊人久久国产一区二区| 亚洲色图av天堂| 精品久久久久久电影网| 王馨瑶露胸无遮挡在线观看| 精品亚洲乱码少妇综合久久| 99精国产麻豆久久婷婷| 国产欧美另类精品又又久久亚洲欧美| 97热精品久久久久久| 成年av动漫网址| 亚洲国产日韩一区二区| 少妇人妻 视频| 亚洲精品国产av成人精品| 别揉我奶头 嗯啊视频| 成人亚洲精品一区在线观看 | 免费黄频网站在线观看国产| 欧美最新免费一区二区三区| 久久久久久久久大av| 天美传媒精品一区二区| 欧美极品一区二区三区四区| 亚洲最大成人中文| 三级男女做爰猛烈吃奶摸视频| 国产成人91sexporn| 卡戴珊不雅视频在线播放| 99热这里只有精品一区| 亚洲国产高清在线一区二区三| 国产av不卡久久| 日韩免费高清中文字幕av| 99热这里只有精品一区| 成年版毛片免费区| 国产成人91sexporn| 狂野欧美激情性bbbbbb| 亚洲精华国产精华液的使用体验| 成人毛片60女人毛片免费| 久久ye,这里只有精品| 国产69精品久久久久777片| 亚洲精品一二三| 成人黄色视频免费在线看| 国产精品久久久久久av不卡| 国产日韩欧美亚洲二区| 日韩av不卡免费在线播放| 99九九线精品视频在线观看视频| 蜜桃亚洲精品一区二区三区| 国产黄片美女视频| 免费不卡的大黄色大毛片视频在线观看| 久久久久九九精品影院| 成人综合一区亚洲| 麻豆久久精品国产亚洲av| 日本猛色少妇xxxxx猛交久久| 精品人妻熟女av久视频| 可以在线观看毛片的网站| 午夜福利在线在线| 黄色怎么调成土黄色| 亚洲av不卡在线观看| 91精品一卡2卡3卡4卡| 少妇人妻一区二区三区视频| 亚洲高清免费不卡视频| videos熟女内射| 91午夜精品亚洲一区二区三区| 舔av片在线| 亚洲图色成人| 午夜福利网站1000一区二区三区| 听说在线观看完整版免费高清| 国产亚洲一区二区精品| 欧美最新免费一区二区三区| 亚洲av福利一区| 建设人人有责人人尽责人人享有的 | 只有这里有精品99| 国产免费视频播放在线视频| 又黄又爽又刺激的免费视频.| 日韩av在线免费看完整版不卡| 日本爱情动作片www.在线观看| 老师上课跳d突然被开到最大视频| 1000部很黄的大片| 搡老乐熟女国产| 亚洲精品中文字幕在线视频 | 最近中文字幕高清免费大全6| 久久久久久久亚洲中文字幕| 亚洲av欧美aⅴ国产| a级毛片免费高清观看在线播放| 国产一级毛片在线| 国内精品美女久久久久久| 一级二级三级毛片免费看| 中文在线观看免费www的网站| 国产精品一区二区性色av| 九色成人免费人妻av| 麻豆乱淫一区二区| 不卡视频在线观看欧美| xxx大片免费视频| 免费看日本二区| 全区人妻精品视频| 青春草亚洲视频在线观看| 国产久久久一区二区三区| 日韩一本色道免费dvd| 黄色怎么调成土黄色| 嫩草影院新地址| 精品人妻一区二区三区麻豆| 男的添女的下面高潮视频| 黄色怎么调成土黄色| 久久久a久久爽久久v久久| 午夜福利在线在线| 国产成人a∨麻豆精品| 国产日韩欧美在线精品| 欧美日韩在线观看h| 国产亚洲精品久久久com| 亚洲国产精品999| 99热网站在线观看| 丝袜脚勾引网站| 国产白丝娇喘喷水9色精品| 久久ye,这里只有精品| 高清在线视频一区二区三区| 只有这里有精品99| 热99国产精品久久久久久7| 久久鲁丝午夜福利片| 中文精品一卡2卡3卡4更新| 久久精品久久久久久久性| 成人亚洲精品一区在线观看 | 国产黄频视频在线观看| 国产成人a∨麻豆精品| av在线亚洲专区| 五月玫瑰六月丁香| 国产精品一二三区在线看| 免费观看av网站的网址| 青青草视频在线视频观看| 又黄又爽又刺激的免费视频.| av免费观看日本| 精品一区二区三区视频在线| 18+在线观看网站| 3wmmmm亚洲av在线观看| 乱系列少妇在线播放| 久久久久久久久久成人| 自拍欧美九色日韩亚洲蝌蚪91 | 777米奇影视久久| 国产欧美日韩一区二区三区在线 | eeuss影院久久| 2021少妇久久久久久久久久久| 免费看光身美女| 国产男女内射视频| 欧美3d第一页| 99视频精品全部免费 在线| 久久久久久久精品精品| 成人一区二区视频在线观看| 精品久久久久久电影网| 免费黄色在线免费观看| 成年版毛片免费区| 在线观看国产h片| 黄色欧美视频在线观看| 亚洲成人中文字幕在线播放| 黄色视频在线播放观看不卡| 国产成人精品福利久久| 国产av不卡久久| 丝瓜视频免费看黄片| 日产精品乱码卡一卡2卡三| 色播亚洲综合网| 欧美变态另类bdsm刘玥| 日韩大片免费观看网站| 干丝袜人妻中文字幕| 亚洲av日韩在线播放| 少妇猛男粗大的猛烈进出视频 | 国产亚洲午夜精品一区二区久久 | av免费在线看不卡| 女人被狂操c到高潮| 国产免费福利视频在线观看| 嫩草影院新地址| 亚洲av中文字字幕乱码综合| 免费av毛片视频| 极品教师在线视频| av在线观看视频网站免费| 亚洲成人久久爱视频| 免费观看av网站的网址| 男男h啪啪无遮挡| 亚洲国产日韩一区二区| 中国国产av一级| 亚洲成人一二三区av| 色网站视频免费| 一区二区三区四区激情视频| 99久久精品一区二区三区| 久久久久久久久大av| 亚洲精品成人av观看孕妇| 大话2 男鬼变身卡| 亚洲va在线va天堂va国产| 男人添女人高潮全过程视频| 99视频精品全部免费 在线| 一级毛片 在线播放| 在线免费十八禁| videos熟女内射| 你懂的网址亚洲精品在线观看| 久久久午夜欧美精品| 国产精品国产三级国产av玫瑰| 观看美女的网站| 在线免费十八禁| 成人午夜精彩视频在线观看| 国产伦在线观看视频一区| 麻豆久久精品国产亚洲av| 在线观看三级黄色| 18禁在线播放成人免费| 夜夜爽夜夜爽视频| 亚洲av成人精品一二三区| 亚洲精品第二区| 亚洲欧美日韩东京热| 成人国产av品久久久| 99精国产麻豆久久婷婷| 一区二区三区四区激情视频| 亚洲人成网站在线观看播放| 亚洲自偷自拍三级| 成人鲁丝片一二三区免费| 久久热精品热| 一本一本综合久久| 91狼人影院| 国产亚洲91精品色在线| 中文字幕人妻熟人妻熟丝袜美| 亚洲成人精品中文字幕电影| 午夜爱爱视频在线播放| 寂寞人妻少妇视频99o| 卡戴珊不雅视频在线播放| 人妻制服诱惑在线中文字幕| 免费大片18禁| 熟女av电影| 国产视频首页在线观看| 亚洲av中文字字幕乱码综合| 欧美精品国产亚洲| 两个人的视频大全免费| 国产av国产精品国产| 久久精品国产亚洲av涩爱| 国产高清国产精品国产三级 | 少妇的逼水好多| 美女国产视频在线观看| 九九爱精品视频在线观看| 日韩av不卡免费在线播放| 麻豆成人午夜福利视频| 国产午夜福利久久久久久| 男男h啪啪无遮挡| 黄色一级大片看看| 在线观看美女被高潮喷水网站| 人妻一区二区av| 我要看日韩黄色一级片| 久久久久久久亚洲中文字幕| 中国三级夫妇交换| 观看免费一级毛片| 亚洲国产av新网站| 中文字幕制服av| 国产一区二区在线观看日韩| 少妇 在线观看| 男女边摸边吃奶| 久久国内精品自在自线图片| 波野结衣二区三区在线| 性插视频无遮挡在线免费观看| 国产精品.久久久| 女人十人毛片免费观看3o分钟| 伊人久久精品亚洲午夜| 精品少妇久久久久久888优播| 亚洲在久久综合| 99九九线精品视频在线观看视频| 赤兔流量卡办理| 综合色av麻豆| 午夜视频国产福利| 2018国产大陆天天弄谢| 日韩视频在线欧美| 中文精品一卡2卡3卡4更新| 亚洲人成网站高清观看| 亚洲av.av天堂| 精品99又大又爽又粗少妇毛片| 国产精品无大码| 亚洲成人一二三区av| 街头女战士在线观看网站| 国产精品一二三区在线看| 亚洲真实伦在线观看| 99热这里只有是精品在线观看| 只有这里有精品99| 久久99热6这里只有精品| av播播在线观看一区| a级毛色黄片| 久久99精品国语久久久| 三级国产精品片| 亚洲自拍偷在线| 在线亚洲精品国产二区图片欧美 | 色吧在线观看| 免费黄频网站在线观看国产| 亚洲精品乱久久久久久| 又黄又爽又刺激的免费视频.| 18禁裸乳无遮挡动漫免费视频 | 亚洲国产精品专区欧美| 直男gayav资源| 色吧在线观看| 亚洲天堂国产精品一区在线| 中国三级夫妇交换| 观看免费一级毛片| 亚洲自偷自拍三级| 可以在线观看毛片的网站| 国产精品三级大全| 国产淫语在线视频| 亚洲av福利一区| 国产伦精品一区二区三区四那| 夫妻午夜视频| 国产成人午夜福利电影在线观看| 听说在线观看完整版免费高清| 97超碰精品成人国产| 舔av片在线| 尾随美女入室| 久久综合国产亚洲精品| 亚洲国产日韩一区二区| 国产免费视频播放在线视频| 久久女婷五月综合色啪小说 | av在线app专区| 在线观看三级黄色| 午夜免费观看性视频| av一本久久久久| 亚洲欧美日韩东京热| 一二三四中文在线观看免费高清| 免费av不卡在线播放| 国产亚洲av嫩草精品影院| 成年女人看的毛片在线观看| 简卡轻食公司| 亚洲精品,欧美精品| 久久久久久久国产电影| 日日撸夜夜添| 午夜精品一区二区三区免费看| 国产欧美另类精品又又久久亚洲欧美| 国产日韩欧美亚洲二区| 国产精品久久久久久精品电影| 国产成人精品一,二区| 国产成人福利小说| 中国美白少妇内射xxxbb| 一区二区三区四区激情视频| 夫妻性生交免费视频一级片| 99久久精品国产国产毛片| 交换朋友夫妻互换小说| 99热这里只有是精品在线观看| 日韩中字成人| 老女人水多毛片| 夫妻午夜视频| 日本熟妇午夜| 一区二区三区乱码不卡18| 永久网站在线| 国产69精品久久久久777片| 亚洲av成人精品一二三区| 成人美女网站在线观看视频| 国产精品爽爽va在线观看网站| 一个人看视频在线观看www免费| 王馨瑶露胸无遮挡在线观看| 成人毛片60女人毛片免费| 少妇熟女欧美另类| 久久精品久久久久久噜噜老黄| 高清在线视频一区二区三区| 国产黄a三级三级三级人|