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

    基于機器學習的瀝青路面壓實度質(zhì)量評估

    2024-05-21 00:19:19趙琪張健張智民陳鎮(zhèn)文劉澤佳周立成劉逸平
    濟南大學學報(自然科學版) 2024年3期
    關(guān)鍵詞:動力學模型壓實度機器學習

    趙琪 張健 張智民 陳鎮(zhèn)文 劉澤佳 周立成 劉逸平

    文章編號:1671-3559(2024)03-0331-10DOI:10.13349/j.cnki.jdxbn.20240318.001

    摘要: 為了提高瀝青路面壓實度預(yù)測的準確率,基于機器學習對瀝青路面進行壓實度質(zhì)量評估;通過對比實際工程中不同壓路機在瀝青路面上的振動規(guī)律,構(gòu)建振動壓路機-瀝青路面系統(tǒng)動力學模型,利用Simulink數(shù)值仿真軟件對所構(gòu)建的模型進行仿真,設(shè)計8種工況,對比壓實度計值、 壓實度控制值、 機械驅(qū)動功率、滾輪綜合剛度、 填筑體能量、 單位體積壓實功率6個壓實度質(zhì)量評估指標在各工況下的適用性差異;采用支持向量機、 邏輯回歸、 k最近鄰、 決策樹、 樸素貝葉斯法5種傳統(tǒng)的機器學習方法對各工況下的壓實度質(zhì)量評估指標樣本進行訓練,對比壓實度預(yù)測準確率;設(shè)計壓路機在不同碾壓次數(shù)時的碾壓路線,對比分別采用最優(yōu)壓實度質(zhì)量評估指標和單一壓實度質(zhì)量評估指標的壓實度預(yù)測準確率。結(jié)果表明:不同壓實度質(zhì)量評估指標在不同工況下的適用性不同,即使在同一種工況下,不同碾壓次數(shù)時的適用性也存在差異;采用最優(yōu)壓實度質(zhì)量評估指標代替單一壓實度質(zhì)量評估指標,壓實度預(yù)測準確率提高5.8%; 在5種傳統(tǒng)的機器學習方法中,樸素貝葉斯法預(yù)測最優(yōu)壓實度質(zhì)量評估指標類型的準確率最高,為96.22%。

    關(guān)鍵詞: 瀝青路面; 壓實度質(zhì)量評估; 機器學習; 壓實度; 動力學模型

    中圖分類號: U416.217

    文獻標志碼: A

    開放科學識別碼(OSID碼):

    Quality Evaluation on Compaction of

    Asphalt Pavements Based on Machine Learning

    ZHAO Qi1, ZHANG Jian2, ZHANG Zhimin2, CHEN Zhenwen2, LIU Zejia1, ZHOU Licheng1, LIU Yiping1

    (1. School of Civil Engineering and Transportation, South China University of Technology, Guangzhou 510641, Guangdong, China;

    2. Guangzhou Chengtou Road and Bridge Engineering Co., Ltd., Guangzhou 510199, Guangdong, China)

    Abstract: To improve accuracy of asphalt pavement compaction prediction, quality evaluation on compaction asphalt pavements was conducted on the basis of machine learning. By comparing vibration patterns of different rollers on asphalt pavements in practical engineering, a dynamical model of vibration roller-asphalt pavement system was constructed.? The constructed model was simulated by using Simulink numerical simulation software, and eight working conditions were designed to compare applicability differences of six quality evaluation indicators of compaction under different working conditions including compaction meter value, compaction control value, machine drive power, roller-integrated stiffness, filling energy, and unit volume compaction power. Five traditional machine learning methods, namely support vector machine, logistic regression, k-nearest neighbor, decision tree, and naive Bayes, were used to train quality evaluation indicator samples of compaction under different working conditions, and prediction accuracy of compaction was compared. Rolling routes of the roller with different rolling times were designed to compare predict accuracies of compaction using the

    收稿日期: 2023-03-14????????? 網(wǎng)絡(luò)首發(fā)時間:2024-03-19T16:00:19

    基金項目: 國家自然科學基金項目(11972162);廣東省自然科學基金項目(2023A1515012942)

    第一作者簡介: 趙琪(1998—),男,滿族,吉林四平人。碩士研究生,研究方向為瀝青路面壓實度的智能監(jiān)測。E-mail: 3487554692@qq.com。

    通信作者簡介: 劉澤佳(1975—),男,山東平度人。副教授,博士,碩士生導師,研究方向為固體力學。E-mail: zjliu@scut.edu.cn。

    網(wǎng)絡(luò)首發(fā)地址: https://link.cnki.net/urlid/37.1378.N.20240318.1842.002

    optimal quality evaluation indicators of compaction and the single quality evaluation indicators of compaction. The results show that the applicability of different quality evaluation indicators of compaction is different under each working condition. Even under the same working condition, the applicability with different rolling times is also different. Using the optimal quality evaluation indicators of compaction to evaluate indicators instead of the single quality evaluation indicators of compaction, the prediction accuracy of compaction is increased by 5.8%. Among the five traditional machine learning methods, naive Bayes method has the highest accuracy of predicting types of the optimal compaction quality evaluation indicators, which is 96.22%.

    Keywords: asphalt pavement; quality evaluation on compaction; machine learning; compaction; dynamical model

    瀝青路面施工工期較短,平整度好,行車舒適性好且噪聲小,同時養(yǎng)護方便。在瀝青路面建設(shè)施工過程中,壓實是瀝青路面施工過程中的最后一道工序,壓實度是壓實過程中嚴格控制的指標之一[1]。壓實度對瀝青路面施工質(zhì)量、 使用性能和使用壽命都有重要影響,在一定程度上決定了瀝青路面的抗水損害能力、 抗車轍能力和耐久性,因此壓實度質(zhì)量評估被高度重視。

    傳統(tǒng)的檢測壓實度的方法一般有灌砂法、 鉆芯取樣法等[2]。這些方法存在很多不足之處,例如:利用灌砂法和鉆芯取樣法進行壓實度的測定時,通常是對某個點進行檢測的,隨機性大,以點蓋面的做法缺乏代表性;檢測的結(jié)果具有滯后性,無法實時得到準確的壓實度測量結(jié)果,最終導致人力、 物力、 財力浪費;傳統(tǒng)的檢測壓實度的方法是有損的,對道路會造成一定程度的破壞。

    由于傳統(tǒng)壓實檢測技術(shù)存在諸多缺點,因此智能壓實檢測技術(shù)應(yīng)運而生[3]。智能壓實檢測技術(shù)是采用計算機技術(shù)、高精密傳感技術(shù)以及定位技術(shù)等高新技術(shù)集成的新型壓實檢測技術(shù)[4]。智能壓實檢測技術(shù)對于解決上述壓實質(zhì)量控制的不足具有非常顯著的優(yōu)勢,但是目前在智能壓實檢測領(lǐng)域,仍有一些問題需要解決: 1)大多數(shù)智能壓實檢測設(shè)備仍基于碾壓次數(shù)或單一壓實度質(zhì)量評估指標(簡稱單一指標)評價壓實度?;谀雺捍螖?shù)的方法忽略了其他因素,具有較大誤差,不宜采用;由于不同指標適用的工況不同,因此基于單一指標的方法對所有壓實工況均采用單一指標時也會產(chǎn)生較大誤差,甚至是錯誤的預(yù)測結(jié)論。 2)利用機器學習方法預(yù)測壓實度的方法改進了傳統(tǒng)方法的不足之處,但是需要進行大量的現(xiàn)場施工試驗從而建立數(shù)據(jù)庫,此過程需要耗費大量的人力和物力。

    為了提高壓實度預(yù)測準確率, 本文中對不同壓實度質(zhì)量評估指標在不同工況下的適用性進行探討, 并結(jié)合相應(yīng)工況選取最優(yōu)壓實度質(zhì)量評估指標(簡稱最優(yōu)指標), 構(gòu)建振動壓路機-瀝青路面系統(tǒng)動力學模型(簡稱本文模型), 結(jié)合MATLAB/Simulink軟件與機器學習, 對本文模型進行仿真, 設(shè)計8種不同類型的工況, 分別比較采用最優(yōu)指標預(yù)測和采用單一指標預(yù)測壓實度的準確率, 并結(jié)合路基壓實現(xiàn)場試驗, 進一步驗證最優(yōu)指標預(yù)測相對于單一指標預(yù)測的優(yōu)勢。

    濟南大學學報(自然科學版)第38卷

    第3期趙琪,等:基于機器學習的瀝青路面壓實度質(zhì)量評估

    1? 模型構(gòu)建

    1.1? 智能壓實檢測技術(shù)

    智能壓實檢測技術(shù)是用于土壤、 骨料和瀝青混合料壓實檢測的技術(shù)。該技術(shù)通過測量壓路機振動輪的加速度信號,計算得到相應(yīng)的壓實度質(zhì)量評估指標,以此反映填筑體的壓實程度[5],因此構(gòu)建振動加速度與壓實度的關(guān)系模型尤為重要。智能壓實檢測技術(shù)工作原理如圖1所示。加速度傳感器豎直置于振動輪中軸,以測量振動輪的加速度信號。隨著填筑體壓實度的增加,振動輪所受地面的反作用力相應(yīng)增加,加速度也增加。將所測加速度信號上傳至位于駕駛室內(nèi)的控制系統(tǒng),從而計算相應(yīng)的智能壓實測試值,最終得到實時預(yù)測的壓實度。根據(jù)智能壓實檢測技術(shù)工作原理,可以通過測量加速度信號評估填筑體的壓實狀態(tài)。

    1.2? 壓實度質(zhì)量評估指標

    僅采用加速度評價壓實度因所考慮影響因素有限,通常產(chǎn)生較大誤差。為了有效地提高預(yù)測準確率,在加速度的基礎(chǔ)上對壓實度質(zhì)量評估指標進行改進。指標的可靠性是直接影響智能壓實檢測技術(shù)適用性的關(guān)鍵因素。

    壓實度質(zhì)量評估指標最早出現(xiàn)于20世紀70年代,即基于對加速度信號進行傅里葉變換后,加速度頻域上2次諧波幅值與基頻幅值之比與壓實度之間的關(guān)系,以壓實度計值(compaction meter value, CMV)作為壓實度質(zhì)量評估指標,記為vcm。其他頻域上的振幅也對壓實度產(chǎn)生一定的影響,因此更多頻譜成分被引入畸變程度的計算中,以壓實度控制值(compaction control value, CCV)作為壓實度質(zhì)量評估指標,記為vcc[6]。vcc是vcm的擴展,不僅考慮了振動輪垂直加速度在頻域上的1、 2次諧波振幅,同時考慮了能夠反映諧波形狀的0.5、 1.5、 2.5、 3次諧波振幅。

    Rinehart等[7]通過建立壓實度和壓路機消耗能量的關(guān)系式,提出以機械驅(qū)動功率(machine drive power)作為壓實度質(zhì)量評估指標,記為Pmd,計算式[7]為

    Pmd=P-Wvsin α+z··g-(vl1+l2) ,(1)

    式中: P為壓路機總驅(qū)動功率; W為壓路機重力; v為壓路機行駛速度; α為壓路機的爬坡角; z··為振動輪振動加速度;g為重力加速度; l1、 l2為機械自身能量損耗率。

    Anderegg等[8]從振動力學的角度出發(fā),通過求解二自由度動力學模型的二階微分運動方程,提出以滾輪綜合剛度(roller-integrated stiffness)作為壓實度質(zhì)量評估指標,記為ks,計算式[8]為

    ks=Ω2m2+m0e0cos z

    ,(2)

    式中: Ω為壓路機激振頻率; m2為振動輪質(zhì)量; m0為振動輪中偏心塊質(zhì)量; e0為振動輪中偏心塊的偏心距; 為偏心力與振動輪位移之間的相位滯后角; z為振動輪位移。

    Krber等[9]采用作用力與位移的乘積表征能量,將2個連續(xù)的振動周期內(nèi)的抵抗力和振動輪位移的時間歷程進行積分,得到傳遞到填筑體的能量,提出以填筑體能量作為壓實度質(zhì)量評估指標,記為ω,計算式[9]為

    ω=∫2T[-m1z··+(m1+m2)g+m0e0Ω2]z·dt ,(3)

    式中: m1為機架質(zhì)量; z·為振動輪振動速度; t為振動輪激振時間; T為振動輪激振周期。

    Liu等[10]從壓路機做功的角度出發(fā),提出以單位體積壓實功率作為壓實度質(zhì)量評估指標,記為Pvol,計算式[10]為

    Pvol=2AQ+πF04ΩNbhv ,(4)

    式中: A為振動輪振幅; Q為振動輪豎向荷載; F0為振動輪激振力; N為振動輪碾壓次數(shù); b為振動輪寬度; h為碾壓層厚度。

    目前,國內(nèi)外已有10余種壓實度質(zhì)量評估指標應(yīng)用于實際工程,其中vcm、 vcc、 Pmd、 ks、 ω、 Pvol這6個壓實度質(zhì)量評估指標應(yīng)用效果較好,但是各指標在不同工況下的適用性和敏感性方面的研究還有所欠缺[11]。本文中對6個壓實度質(zhì)量評估指標在不同工況的適用性進行進一步研究,從而實現(xiàn)提高壓實度預(yù)測準確率的最終目標。

    1.3? 本文模型及其動力學方程

    路基的壓實工藝分為初壓、 復壓和終壓3個階段,而振動壓實主要集中在復壓階段,這也是本文中的重點研究階段。目前已有學者針對不同類型的振動壓路機提出相應(yīng)的數(shù)學模型[12],本文中通過對比實際工程中不同壓路機在瀝青路面上的振動規(guī)律,構(gòu)建本文模型,如圖2所示。

    m1—機架質(zhì)量; m2—振動輪質(zhì)量; m3—隨振瀝青質(zhì)量;

    k1、 c1—減震器剛度、 阻尼; k2、 c2—瀝青剛度、 阻尼;

    F0—振動輪激振力; ω0—振動輪激振角頻率;

    t—振動輪激振時間; x1、 x2、 x3—m1、 m2、 m3的自由度;

    Fs—地面的抵抗力。

    本文模型的構(gòu)建基于以下假設(shè): 1)只考慮豎直方向的振動,本文模型的上半部分代表壓路機的特性,下半部分代表瀝青材料的特性; 2)將壓路機的機架、 振動輪設(shè)為一定質(zhì)量的剛體,并且機架與振動輪之間用無質(zhì)量的彈簧-阻尼單元作為減震系統(tǒng)進行連接; 3)振動輪與地面之間緊密接觸,不存在跳振,即x1=x2,其中x1、 x2分別為m1、 m2的自由度。根據(jù)牛頓第二定律,本文模型的動力學方程為

    m1x··1+c1x·1+k1x1-c1x·2-k1x2=0 ,

    (m2+m3)x··2+(c1+c2)x·+(k1+k2)x2-c1x·1-

    k1x1=F0sinω0t ,

    F0=m0e0ω20 ,(5)

    式中: m3為隨振瀝青質(zhì)量; k1、 c1分別為減震器剛度、 阻尼; k2、 c2分別為瀝青剛度、 阻尼; ω0為振動輪的激振角頻率。式(5)的矩陣形式為

    MX··+CX·+KX=F(t) ,(6)

    式中:M=m10

    0m2+m3為質(zhì)量矩陣; C=

    c1-c1

    -c1c1+c2為阻尼矩陣; K=k1-k1

    -k1k1+k2為剛度矩陣; X、 X·、 X··分別為自由度位移、 速度、 加速度向量;

    F(t)為激振力矩陣。

    2? 模型仿真

    2.1? 狀態(tài)空間仿真

    采用MATLAB/Simulink軟件對本文模型進行機械振動分析仿真,較常用的方法有積分模塊仿真、 狀態(tài)空間仿真、 S函數(shù)模塊仿真3種方法[13]。由于狀態(tài)空間仿真有諸多優(yōu)點,如仿真過程簡單、 可讀性好、 容易

    操作等,因此本文中采用該方法進行仿真。

    狀態(tài)空間仿真采用狀態(tài)方程模塊State-Space,輸入為正弦信號,輸出為位移、 速度和加速度信號,并引入對信號進行快速傅里葉變換的Discrete模塊,從而有利于根據(jù)加速度時域信號求解vcm、 vcc, ΔuΔt為計算輸入對時間的變化率的微分模塊。為了便于觀察和分析,將仿真結(jié)果輸出至Workspace工作空間。本文模型仿真過程如圖3所示。

    State-Space模塊—狀態(tài)方程模塊;

    Discrete模塊—對信號進行快速傅里葉變換的模塊;

    ΔuΔt模塊—計算輸入對時間的變化率的微分模塊。

    2.2? 仿真參數(shù)設(shè)置

    在本文模型仿真的過程中,需要對各模塊以及仿真參數(shù)進行設(shè)置。參考文獻[14]設(shè)置本文模型的仿真參數(shù),如表1所示,其中m3設(shè)置為m2的30%[15]。

    設(shè)置2種振幅分別為0.25、 1.00 mm的信號波作為輸入,本文模型中振動輪的加速度-時間曲線仿真結(jié)果如圖4所示。由圖可知,輸出的2個振動輪加速度信號穩(wěn)定并且大致符合振幅的倍數(shù)關(guān)系,可以初步驗證仿真結(jié)果的可靠性。

    2.3? 工況設(shè)置

    根據(jù)文獻[16], k1、 c1根據(jù)相應(yīng)壓路機說明書確定范圍, k2、 c2根據(jù)瀝青材料級配確定范圍,并在此范圍內(nèi)梯度設(shè)置k2的數(shù)值。本文模型仿真參數(shù)值設(shè)置完成后,設(shè)置正弦波信號幅值和頻率,輸入仿真參數(shù)值并運行即可得到振動輪的加速度。再根據(jù)AC-25型瀝青路面壓實度-振動加速度-溫度關(guān)系擬合結(jié)果[14]計算瀝青路面壓實度,即

    振動壓路機-瀝青路面系統(tǒng)動力學模型中

    振動輪的加速度-時間曲線仿真結(jié)果

    K=0.826exp(0.009θ)z··-0.555θ+60.934

    ,(7)

    式中: K為瀝青路面壓實度; θ為瀝青路面攝氏溫度。

    通過給定溫度值得到相應(yīng)的壓實度值, 再結(jié)合壓路機的機械參數(shù)得到不同的壓實度質(zhì)量評估指標值, 并通過計算不同工況下壓實度質(zhì)量評估指標值與壓實度值之間的相關(guān)系數(shù)得到不同工況下的最優(yōu)指標。 在此過程中, 根據(jù)正弦波信號振幅的不同, 將工況設(shè)置為低、 高幅2種, 振幅分別為1、 2? mm; 將頻率設(shè)置為低、 高頻2種,頻率分別位30、 50 Hz;將溫度設(shè)置為低、 高溫2種,溫度分別為120、 150 ℃。將以上3種物理參數(shù)結(jié)合形成8種工況:低溫低頻低幅、 高溫低頻低幅、 低溫低頻高幅、 低溫高頻低幅、 低溫高頻高幅、 高溫低頻高幅、 高溫高頻低幅、 高溫高頻高幅,分別記為工況1、 2、 …、 8。通過梯度改變k2值,從而形成320個樣本。

    2.4? 仿真結(jié)果分析

    8種工況下壓實度與vcm、 vcc、 Pmd、 ks、 ω、 Pvol這6個壓實度質(zhì)量評估指標之間的關(guān)系如圖5所示。由圖可知:在工況1下,vcm、 vcc指標曲線突升和突降較明顯,與壓實度的相關(guān)性較差,而另外4種指標和壓實度的相關(guān)性較好;在工況7下, ω、 Pmd指標曲線在后期有突升和突降,而vcm、 vcc曲線相對平滑。由此可知,不同工況下,不同指標與壓實度的相關(guān)性不同,如果用不同工況下最優(yōu)指標評價壓實度,將能提升預(yù)測準確率。

    對不同指標和壓實度之間的相關(guān)性進行數(shù)值分析, 根據(jù)碾壓次數(shù)的不同, 將工況分為4個階段, 并且統(tǒng)計不同階段時不同指標與壓實度的相關(guān)系數(shù)ρ。ρ的絕對值越接近1, 說明指標與壓實度之間的相關(guān)性越強, 越適宜作為評價此工況下壓實度的最優(yōu)指標。 在統(tǒng)計過程中發(fā)現(xiàn): 3個能量類的指標中, ω與壓實度的ρ大于另外2個指標,因此將ω代表能量類指標參與對比分析。統(tǒng)計8種工況下不同壓實度質(zhì)量評估指標與壓實度的ρ,結(jié)果如表2所示。標注*的數(shù)值對應(yīng)的該工況下ρ最大,即該工況下壓實度質(zhì)量評估指標為最優(yōu)指標。由圖可知,最優(yōu)指標分布沒有規(guī)律,進一步驗證了不同工況下不同指標與壓實度的相關(guān)性不同的結(jié)論。

    統(tǒng)計8種工況下采用最優(yōu)指標與單一指標的相關(guān)系數(shù)ρ,針對不同工況下4種碾壓次數(shù)均選用單一指標預(yù)測,并取平均值計算相應(yīng)的ρ,每種工況下ρ最大的指標即該工況下的最優(yōu)指標。然后針對不同工況下4種碾壓次數(shù)均選用最優(yōu)指標預(yù)測,并取平均值計算相應(yīng)的ρ,最后計算不同工況下采用最優(yōu)指標與采用單一指標的壓實度預(yù)測準確率之間的差值,即采用最優(yōu)指標比單一指標預(yù)測提升的準確率,范圍為1%~22.5%,統(tǒng)計結(jié)果如表3所示。由表可知:對于不同工況下采用最優(yōu)指標比采用單一指標的壓實度預(yù)測準確率均有一定的提升;不同壓實度質(zhì)量評估指標在8種工況下的適用性有所不同,在不同工況下選取最優(yōu)指標預(yù)測壓實度,壓實度預(yù)測準確率平均提高9.575%。

    3? 最優(yōu)指標類型預(yù)測

    采用支持向量機(support vector machine,SVM)、 邏輯回歸(logistic regression)、 k最近鄰(k-nearest neighbor,KNN)、 決策樹(decision tree,DT)、 樸素貝葉斯(nave Bayes)法5種傳統(tǒng)的機器學習方法分別進行訓練學習。

    對8種工況下的320個樣本進行訓練學習。每個樣本對應(yīng)一組標簽,包括振動輪振幅A、 壓路機的激振頻率Ω、 瀝青路面攝氏溫度θ、 瀝青剛度k2、 瀝青阻尼c2、 振動輪振動加速度z··、 振動輪碾壓次數(shù)N。每個樣本同時對應(yīng)一組特征值為響應(yīng)的vcm、 vcc、 ks、 ω這4個指標中最優(yōu)指標的類型,將特征值與標簽一一對應(yīng)組合在一起,得到樣本集。然后將樣本集按照測試集與訓練集中樣本個數(shù)之比為3∶7的比例隨機分配。最終對4個指標進行分類學習,以實現(xiàn)根據(jù)已有樣本在不同工況下尋求最優(yōu)指標,并以成功尋求最優(yōu)指標的概率作為機器學習算法的壓實度預(yù)測準確率。

    5種傳統(tǒng)機器學習方法預(yù)測最優(yōu)指標類型準確率的統(tǒng)計結(jié)果如表4所示。 由表可知, 決策樹和樸素貝葉斯法的壓實度預(yù)測準確率較高, 其中樸素貝葉斯法的壓實度預(yù)測準確率最高, 可達96.22%, 因此在實際工程中可以使用樸素貝葉斯法尋求最優(yōu)指標。

    4? 路基壓實現(xiàn)場試驗

    4.1? 試驗條件與儀器

    4.1.1? 試驗條件

    為了對本文模型進行驗證, 開展路基壓實現(xiàn)場試驗, 施工前、 后現(xiàn)場試驗實景如圖6所示。 試驗路段面層采用AC-25型改性瀝青混凝土鋪筑, 厚度為4 cm。 考慮到與正常施工工序相吻合, 采取振

    動輪先靜壓2次后動壓10次的方案, 碾壓軌跡路線圖如圖7所示, 其中箭頭方向為碾壓軌跡方向, 采樣長度、 寬度分別為35、 2.2 m, 分為13個區(qū)域。 考慮到兩端位置處的機械啟停會產(chǎn)生一定的誤差, 啟停區(qū)長度取為4 m。 圖中圓圈的位置代表鉆芯取樣區(qū), 取樣區(qū)長度為2 m。 每個區(qū)域采樣3個點并計算壓實度的平均值。 2個圓圈之間的區(qū)域是壓路機需要換向的位置, 因此將該區(qū)域預(yù)留出來, 并且考慮到壓路機的車長,換向區(qū)長度取為3m。虛線

    位置是壓路機每次要換向時前輪軸所在位置。 壓路機按照碾壓路線行駛, 最終從起始位置駛出實驗區(qū)域。

    4.1.2? 試驗儀器

    試驗數(shù)據(jù)采集儀器包括: 1)加速度傳感器。 試驗時, 加速度傳感器垂直放置于振動軸中間位置, 如圖8所示。 在采集數(shù)據(jù)前, 為了防止加速度傳感器滑動, 先擦拭振動軸上的灰塵, 再用膠水將加速度傳感器固定于振動軸。 2)動態(tài)信號采集儀。 作為中間媒介, 動態(tài)信號采集儀一端與加速度傳感器相連, 另一端通過網(wǎng)線與計算機終端連接, 并把計算機終端放置在駕駛室內(nèi)以便于對壓實度進行實時監(jiān)控。 3)手持式測溫槍。

    振動軸中間位置的加速度傳感器

    4.2? 數(shù)據(jù)獲取與處理

    4.2.1? 數(shù)據(jù)獲取

    連續(xù)采樣10次動壓時加速度時域、頻域信號如圖9所示。由圖9(a)可知: 10次振動信號所含噪聲比例較?。挥蓤D9(b)可知:振動輪的振動頻率為50 Hz,與振動輪的理論振動頻率近似相等,間接驗證了信號采集的可靠性。利用測溫槍測得的數(shù)據(jù)包括:出料溫度為155 ℃;靜壓之前溫度為120 ℃;動壓過程全段路面不同測點溫度為100~110 ℃; 動壓

    過程內(nèi)部溫度為120~130 ℃。壓路機自始至終采用最高擋位振動,根據(jù)8種工況設(shè)置,將本次試驗工況類型設(shè)置為工況5。

    4.2.2? 數(shù)據(jù)處理

    根據(jù)圖9(a)計算10次動壓時加速度的幅值平均值,將5個區(qū)域鉆芯取樣所得的壓實度與計算所得的加速度的幅值平均值一一對應(yīng)。由于第6次動壓區(qū)域中有井蓋導致數(shù)據(jù)錯誤,因此對該數(shù)據(jù)進行合理估值,并對第1、 3、 5、 7、 9次動壓時的壓實度數(shù)據(jù)進行插值,最終得到10次動壓時加速度與壓實度的統(tǒng)計結(jié)果,如表5所示。由表可知,在10次動壓過程中,隨著瀝青路面壓實度的增大,振動輪的加速度逐漸增大,與振動壓實理論相吻合,間接驗證了此路基壓實現(xiàn)場試驗結(jié)果的可靠性。

    通過隨機插值,對加速度和壓實度進行樣本擴充,擴充到320個樣本。將壓路機的機械參數(shù)輸入到本文模型,其中除k2、 c2外都可以通過查找壓路機相應(yīng)型號的說明書得到。k2、 c2標定方法如下:先將壓路機放置在壓實度為100%的瀝青路面上振動并測量加速度,然后分別將預(yù)估的k2、 c2輸入到Simulink軟件,通過不斷迭代得到加速度對應(yīng)的k2、 c2。將c2作為定值,根據(jù)瀝青路面壓實度為100%時的k2值確定k2的變化范圍,從而完成對k2、 c2的標定。

    統(tǒng)計利用樸素貝葉斯算法訓練學習后,碾壓次數(shù)為10時采用最優(yōu)指標與單一指標的壓實度指標相關(guān)系數(shù)ρ,結(jié)果如表6所示。標注*的數(shù)值對應(yīng)的該碾壓次數(shù)時ρ最大,即該碾壓次數(shù)時的壓實度質(zhì)量評估指標為最優(yōu)指標,并以該最優(yōu)指標作為每次采用最優(yōu)指標的壓實度預(yù)測準確率。由表可知:當采用單一指標時,利用ω預(yù)測的壓實度預(yù)測準確率最佳,為74.31%;當采用最優(yōu)指標進行預(yù)測時,壓實度預(yù)測準確率為80.11%,比采用單一指標時的壓實度預(yù)測準確率提高了5.8%;在同一種工況下,不同碾壓次數(shù)時不同指標的適用性存在差異,對于不同工況選取最優(yōu)指標可以使壓實度預(yù)測準確率提高5.8%。

    5? 結(jié)論

    本文中通過對實際工程中瀝青路面上不同壓路機的振動規(guī)律進行對比,構(gòu)建振動壓路機-瀝青路面系統(tǒng)動力學模型,結(jié)合機器學習和Simulink模擬仿真,對瀝青路面壓實度進行質(zhì)量評估,并利用路基壓實現(xiàn)場試驗數(shù)據(jù)對采用最優(yōu)指標代替單一指標后壓實度預(yù)測準確率的可靠性進行驗證,得到以下主要結(jié)論:

    1)不同壓實度質(zhì)量評估指標在8種工況下的適用性存在差異,對于不同工況選取最優(yōu)指標可以使壓實度預(yù)測準確率提高9.575%。

    2)相對于支持向量機、 邏輯回歸、 k最近鄰、 決策樹,樸素貝葉斯法更適合預(yù)測不同工況下最優(yōu)指標類型,在預(yù)測最優(yōu)壓實度質(zhì)量評估指標類型時準確率最高,為96.22%。

    3)路基壓實現(xiàn)場試驗結(jié)果表明,即使在同一種工況下,相對于采用單一指標,不同碾壓次數(shù)時采用最優(yōu)指標預(yù)測能使壓實度預(yù)測準確率提高5.8%。

    參考文獻:

    [1]? WANG L B, 王含笑, 趙千, 等. 智能路面發(fā)展與展望[J]. 中國公路學報, 2019, 32(4): 50.

    [2]? 馬建, 孫守增, 芮海田, 等. 中國筑路機械學術(shù)研究綜述·2018[J]. 中國公路學報, 2018, 31(6): 1.

    [3]? 黃志福, 梁乃興, 趙毅, 等. 路面壓實度自動連續(xù)檢測技術(shù)[J]. 長安大學學報(自然科學版), 2015, 35(6): 24.

    [4]? 徐光輝, CHANG G. 智能壓實測量值的發(fā)展方向[J]. 筑路機械與施工機械化, 2018, 35(4): 19.

    [5]? 葉陽升,蔡德鉤,朱宏偉,等.基于振動能量的新型高速鐵路路基壓實連續(xù)檢測控制指標研究[J].鐵道學報,2020, 42(7): 127.

    [6]? 賈通. 瀝青路面智能壓實系統(tǒng)關(guān)鍵技術(shù)研究[D]. 南京: 東南大學, 2020: 8.

    [7]? RINEHART R V, MOONEY M A. Instrumentation of a roller compactor to monitor vibration behavior during earthwork compaction[J]. Automation in Construction, 2008, 17(2): 149.

    [8]? ANDEREGGR,KAUFMANNK.Intelligentcompactionwith vibratory rollers: feedback control systems in automatic compaction and compaction control[J]. Transportation Research Record: Journal of the Transportation Research Board, 2004, 1868(1): 131.

    [9]? KRBER W, FLOSS E H R, WALLRATH W. Dynamic soil stiffness as quality criterion for soil compaction[M]//Geotechnics for Roads, Rail Tracks, and Earth Structures. Lisse: Balkema, 2001: 196

    [10]? LIU D H, LIN M, LI S. Real-time quality monitoring and control of highway compaction[J]. Automation in Construction, 2016, 62: 118.

    [11]? 吳龍梁. 基于能量耗散的路基連續(xù)壓實控制技術(shù)研究[D]. 北京: 中國鐵道科學研究院, 2020: 20.

    [12]? PAULMICHL I, ADAM C, ADAM D. Parametric study of the compaction effect and the response of an oscillation roller[J]. Proceedings of the Institution of Civil Engineers: Geotechnical Engineering, 2020, 173(4): 285.

    [13]? 代然, 胡水英, 張前, 等. 基于MATLAB/Simulink的垂直振動壓路機動態(tài)性能的仿真研究[J]. 建設(shè)機械技術(shù)與管理, 2015, 28(6): 93.

    [14]? 趙毅, 楊臻, 梁乃興, 等. 瀝青混合料壓實度-振動加速度-溫度關(guān)系模型[J]. 建筑材料學報, 2021, 24(3): 621.

    [15]? 張世英. 關(guān)于振動壓路機設(shè)計中土的參數(shù)的確定方法[J]. 工程機械, 1998(1): 5.

    [16]? 陳博. 路基土壓實度實時檢測系統(tǒng)研究[D]. 西安: 長安大學, 2019: 34.

    (責任編輯:王? 耘)

    猜你喜歡
    動力學模型壓實度機器學習
    簡述灌砂法實測公路路基壓實度的重要性
    研究模糊控制下車輛的側(cè)傾穩(wěn)定性仿真分析
    基于機器學習的圖像特征提取技術(shù)在圖像版權(quán)保護中的應(yīng)用
    基于網(wǎng)絡(luò)搜索數(shù)據(jù)的平遙旅游客流量預(yù)測分析
    時代金融(2016年27期)2016-11-25 17:51:36
    前綴字母為特征在維吾爾語文本情感分類中的研究
    科教導刊(2016年26期)2016-11-15 20:19:33
    房地產(chǎn)投資系統(tǒng)動力學模型分析
    基于支持向量機的金融數(shù)據(jù)分析研究
    兩輪自平衡機器人的動力學模型的分析與建立
    四軸飛行器的結(jié)構(gòu)組成和動力學模型
    新型土壤固化劑在軟土治理中的應(yīng)用初探
    国产欧美日韩一区二区三| 一区二区三区精品91| 国产色视频综合| 精品熟女少妇八av免费久了| 亚洲人成伊人成综合网2020| 亚洲国产中文字幕在线视频| 欧美精品啪啪一区二区三区| 亚洲美女黄片视频| 亚洲精品粉嫩美女一区| 成人手机av| www.精华液| av电影中文网址| 国产精品综合久久久久久久免费| 啦啦啦 在线观看视频| 久久精品成人免费网站| 国产97色在线日韩免费| 国产亚洲精品一区二区www| 亚洲精品在线美女| 亚洲精品粉嫩美女一区| aaaaa片日本免费| 怎么达到女性高潮| 999久久久国产精品视频| 国产精品,欧美在线| 精品久久久久久久末码| 精品第一国产精品| 免费在线观看影片大全网站| 久久香蕉国产精品| 91国产中文字幕| 中文字幕精品亚洲无线码一区 | 一本综合久久免费| 97人妻精品一区二区三区麻豆 | 制服人妻中文乱码| 老司机在亚洲福利影院| 久久国产精品男人的天堂亚洲| 香蕉丝袜av| 在线看三级毛片| 波多野结衣巨乳人妻| 午夜福利高清视频| 黄色视频,在线免费观看| 99精品久久久久人妻精品| 亚洲欧美日韩高清在线视频| 国产99白浆流出| 在线免费观看的www视频| 久久精品91无色码中文字幕| 精品第一国产精品| 国产精品1区2区在线观看.| 日韩 欧美 亚洲 中文字幕| 满18在线观看网站| 久久伊人香网站| 亚洲av五月六月丁香网| 99re在线观看精品视频| www.www免费av| 精品午夜福利视频在线观看一区| 成人国产综合亚洲| 草草在线视频免费看| 大型av网站在线播放| 特大巨黑吊av在线直播 | 亚洲五月天丁香| www.精华液| 搞女人的毛片| 国产主播在线观看一区二区| 麻豆av在线久日| 丰满人妻熟妇乱又伦精品不卡| 曰老女人黄片| 亚洲男人天堂网一区| 久久久久九九精品影院| 国产激情偷乱视频一区二区| 国产亚洲精品第一综合不卡| 欧美在线一区亚洲| 正在播放国产对白刺激| 动漫黄色视频在线观看| 欧美性猛交黑人性爽| 精品免费久久久久久久清纯| 桃色一区二区三区在线观看| 欧美乱色亚洲激情| 久久婷婷人人爽人人干人人爱| 久久久国产精品麻豆| 久久亚洲真实| 国产成人啪精品午夜网站| 黄片小视频在线播放| 午夜福利在线观看吧| 日本三级黄在线观看| 侵犯人妻中文字幕一二三四区| 婷婷六月久久综合丁香| 这个男人来自地球电影免费观看| 黄色毛片三级朝国网站| 午夜福利18| 国产伦人伦偷精品视频| 亚洲第一电影网av| 悠悠久久av| 亚洲精品美女久久av网站| 好男人在线观看高清免费视频 | or卡值多少钱| 精品日产1卡2卡| 欧美日韩瑟瑟在线播放| 91麻豆精品激情在线观看国产| 久久久久九九精品影院| 亚洲av片天天在线观看| 免费看a级黄色片| 国内毛片毛片毛片毛片毛片| 亚洲男人的天堂狠狠| av免费在线观看网站| 免费av毛片视频| 这个男人来自地球电影免费观看| 成人精品一区二区免费| www日本黄色视频网| 日韩av在线大香蕉| 高潮久久久久久久久久久不卡| 国产精品综合久久久久久久免费| 国产精品精品国产色婷婷| 久久精品国产99精品国产亚洲性色| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产看品久久| 国产精品av久久久久免费| 亚洲成国产人片在线观看| 欧美精品啪啪一区二区三区| 黄片小视频在线播放| 午夜精品在线福利| 国产亚洲精品综合一区在线观看 | 免费高清在线观看日韩| 欧美一级毛片孕妇| 别揉我奶头~嗯~啊~动态视频| 午夜久久久在线观看| 欧美中文日本在线观看视频| 2021天堂中文幕一二区在线观 | 老司机深夜福利视频在线观看| 日韩免费av在线播放| 亚洲av成人不卡在线观看播放网| 亚洲专区中文字幕在线| 亚洲欧美日韩无卡精品| 香蕉av资源在线| 91国产中文字幕| 久久伊人香网站| 9191精品国产免费久久| 在线观看66精品国产| 2021天堂中文幕一二区在线观 | 国产精品一区二区免费欧美| 91在线观看av| 亚洲中文字幕一区二区三区有码在线看 | 亚洲成人国产一区在线观看| 免费观看人在逋| 黄色片一级片一级黄色片| 国产在线精品亚洲第一网站| 麻豆一二三区av精品| 亚洲 国产 在线| 岛国视频午夜一区免费看| 一本久久中文字幕| 熟女少妇亚洲综合色aaa.| 搞女人的毛片| 国产成人系列免费观看| 成人免费观看视频高清| 又大又爽又粗| 日本撒尿小便嘘嘘汇集6| 久久久精品国产亚洲av高清涩受| 亚洲精品国产一区二区精华液| 一区二区三区国产精品乱码| 婷婷精品国产亚洲av| 村上凉子中文字幕在线| 亚洲精品在线美女| 亚洲第一av免费看| 亚洲av电影在线进入| 欧美乱妇无乱码| 级片在线观看| 欧美日韩瑟瑟在线播放| 精品人妻1区二区| 美女高潮到喷水免费观看| 大香蕉久久成人网| 成年女人毛片免费观看观看9| 麻豆一二三区av精品| 男女午夜视频在线观看| 国产av又大| 成人精品一区二区免费| 久久精品国产清高在天天线| 成人国语在线视频| 国产一区二区三区视频了| 国产精品av久久久久免费| 欧洲精品卡2卡3卡4卡5卡区| 男人舔奶头视频| 69av精品久久久久久| 亚洲欧美日韩高清在线视频| 男女下面进入的视频免费午夜 | 少妇熟女aⅴ在线视频| 国产蜜桃级精品一区二区三区| 久9热在线精品视频| 亚洲一区二区三区不卡视频| 亚洲av成人av| 亚洲国产中文字幕在线视频| 午夜免费成人在线视频| 一夜夜www| 久久国产精品人妻蜜桃| 成人三级黄色视频| 中文字幕人妻熟女乱码| 国产精品1区2区在线观看.| 亚洲熟妇熟女久久| 久久久久国产一级毛片高清牌| 夜夜看夜夜爽夜夜摸| 精品久久久久久久毛片微露脸| 欧美日韩亚洲国产一区二区在线观看| 好男人电影高清在线观看| 国产成人精品久久二区二区91| 最新美女视频免费是黄的| 久久欧美精品欧美久久欧美| 免费人成视频x8x8入口观看| 亚洲在线自拍视频| av欧美777| 欧美成人一区二区免费高清观看 | 久久国产精品人妻蜜桃| 国产又色又爽无遮挡免费看| 热99re8久久精品国产| 国产精品99久久99久久久不卡| 19禁男女啪啪无遮挡网站| 久久久久久亚洲精品国产蜜桃av| 成人av一区二区三区在线看| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久久人妻蜜臀av| 99re在线观看精品视频| 日韩一卡2卡3卡4卡2021年| 欧美不卡视频在线免费观看 | 亚洲精品中文字幕一二三四区| 亚洲精品国产一区二区精华液| 一进一出抽搐gif免费好疼| 亚洲精品中文字幕一二三四区| 精品国产乱子伦一区二区三区| 国产在线观看jvid| 国产成人av教育| 麻豆成人av在线观看| 国产精品精品国产色婷婷| 制服人妻中文乱码| 国产主播在线观看一区二区| 国产精品亚洲av一区麻豆| 久久精品91蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 精品国内亚洲2022精品成人| 国产精品久久久久久精品电影 | 一a级毛片在线观看| 国产真实乱freesex| 亚洲一区中文字幕在线| or卡值多少钱| 久久国产精品人妻蜜桃| 亚洲片人在线观看| 欧美色欧美亚洲另类二区| 夜夜看夜夜爽夜夜摸| 男女之事视频高清在线观看| 欧美日本视频| 久久久国产精品麻豆| 久久午夜综合久久蜜桃| 国产真人三级小视频在线观看| 国产极品粉嫩免费观看在线| 亚洲电影在线观看av| 日本a在线网址| 成人欧美大片| 欧美黄色片欧美黄色片| 午夜久久久在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲色图 男人天堂 中文字幕| 中文字幕精品免费在线观看视频| 他把我摸到了高潮在线观看| 成人三级做爰电影| 妹子高潮喷水视频| www.www免费av| 欧美 亚洲 国产 日韩一| 老汉色av国产亚洲站长工具| 最新美女视频免费是黄的| 精华霜和精华液先用哪个| 欧美日韩一级在线毛片| 国产视频一区二区在线看| 国内毛片毛片毛片毛片毛片| 可以在线观看毛片的网站| 男女做爰动态图高潮gif福利片| 很黄的视频免费| 久久久久久免费高清国产稀缺| 日韩av在线大香蕉| 欧美精品啪啪一区二区三区| 亚洲性夜色夜夜综合| 久久国产精品人妻蜜桃| 国产精品久久久人人做人人爽| 国产亚洲欧美98| 亚洲精品中文字幕在线视频| 午夜免费激情av| 亚洲最大成人中文| 99国产综合亚洲精品| 最近最新中文字幕大全电影3 | 久久性视频一级片| 国产黄片美女视频| 免费搜索国产男女视频| 成人国产综合亚洲| 丝袜美腿诱惑在线| 亚洲av成人一区二区三| 视频区欧美日本亚洲| 最近最新中文字幕大全电影3 | av福利片在线| 不卡一级毛片| 91国产中文字幕| 国产精品一区二区精品视频观看| 亚洲国产精品成人综合色| 欧美日本视频| 一个人免费在线观看的高清视频| 免费一级毛片在线播放高清视频| 亚洲国产欧洲综合997久久, | 久久久水蜜桃国产精品网| 国产人伦9x9x在线观看| 国产欧美日韩一区二区三| 老司机福利观看| 久久人人精品亚洲av| 黄色片一级片一级黄色片| 老汉色∧v一级毛片| 一边摸一边抽搐一进一小说| 日日夜夜操网爽| 最好的美女福利视频网| bbb黄色大片| 一区二区三区激情视频| 俺也久久电影网| 亚洲自偷自拍图片 自拍| 一a级毛片在线观看| 一进一出抽搐动态| 日本撒尿小便嘘嘘汇集6| 在线观看www视频免费| 精品久久久久久久末码| 一二三四在线观看免费中文在| 亚洲七黄色美女视频| 欧美日本视频| 天天添夜夜摸| 91九色精品人成在线观看| 美女高潮到喷水免费观看| 亚洲欧美日韩高清在线视频| 亚洲精品美女久久久久99蜜臀| 色尼玛亚洲综合影院| 国产真人三级小视频在线观看| 桃色一区二区三区在线观看| 欧美日本亚洲视频在线播放| 欧美乱妇无乱码| 好看av亚洲va欧美ⅴa在| 伊人久久大香线蕉亚洲五| 狂野欧美激情性xxxx| 夜夜看夜夜爽夜夜摸| 90打野战视频偷拍视频| 真人一进一出gif抽搐免费| 搞女人的毛片| 欧美日韩中文字幕国产精品一区二区三区| 亚洲avbb在线观看| 黑人巨大精品欧美一区二区mp4| 国产一区二区在线av高清观看| 老司机福利观看| 久久久久久九九精品二区国产 | 日韩欧美免费精品| 久久久国产成人精品二区| 亚洲人成网站在线播放欧美日韩| 成人亚洲精品一区在线观看| 亚洲第一欧美日韩一区二区三区| 国产一区二区三区视频了| 日本免费a在线| 久久热在线av| 韩国精品一区二区三区| 日日干狠狠操夜夜爽| 91老司机精品| 91成人精品电影| 91国产中文字幕| 日本熟妇午夜| 国产成年人精品一区二区| 亚洲成国产人片在线观看| 999久久久精品免费观看国产| 精品免费久久久久久久清纯| 日日夜夜操网爽| 黄色女人牲交| 午夜激情福利司机影院| 一二三四在线观看免费中文在| 久久久久久久午夜电影| 欧美国产日韩亚洲一区| 波多野结衣高清无吗| 18美女黄网站色大片免费观看| 人人妻人人看人人澡| 婷婷六月久久综合丁香| 久久这里只有精品19| 精品不卡国产一区二区三区| 亚洲狠狠婷婷综合久久图片| 中文字幕最新亚洲高清| 韩国精品一区二区三区| 女警被强在线播放| av天堂在线播放| 国产日本99.免费观看| 女人高潮潮喷娇喘18禁视频| 久久草成人影院| 在线视频色国产色| 特大巨黑吊av在线直播 | 黄片播放在线免费| 啦啦啦 在线观看视频| 日日爽夜夜爽网站| 黄色视频,在线免费观看| 香蕉久久夜色| 狠狠狠狠99中文字幕| 欧美精品啪啪一区二区三区| 中亚洲国语对白在线视频| 麻豆成人av在线观看| 十八禁人妻一区二区| 亚洲国产高清在线一区二区三 | 婷婷精品国产亚洲av| 欧美成人一区二区免费高清观看 | 日韩精品中文字幕看吧| 色精品久久人妻99蜜桃| 在线观看www视频免费| 精品欧美国产一区二区三| 国产一区二区三区视频了| 99精品久久久久人妻精品| 91字幕亚洲| 久久久精品国产亚洲av高清涩受| 亚洲国产精品999在线| 久久草成人影院| 国产精品香港三级国产av潘金莲| 国产高清videossex| 午夜免费观看网址| 精品国内亚洲2022精品成人| 嫩草影院精品99| 熟妇人妻久久中文字幕3abv| 69av精品久久久久久| 国产精品免费视频内射| 午夜精品在线福利| 波多野结衣高清无吗| 欧美丝袜亚洲另类 | 成人国产一区最新在线观看| 欧美中文日本在线观看视频| 婷婷六月久久综合丁香| 欧美一级a爱片免费观看看 | 在线看三级毛片| 亚洲国产精品999在线| 亚洲男人的天堂狠狠| 国产蜜桃级精品一区二区三区| 最近最新中文字幕大全电影3 | 亚洲天堂国产精品一区在线| 99在线人妻在线中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 一边摸一边做爽爽视频免费| 男人操女人黄网站| 国产97色在线日韩免费| 在线视频色国产色| 国内久久婷婷六月综合欲色啪| 久久久久久久久久黄片| 一级毛片高清免费大全| 男女之事视频高清在线观看| 在线观看66精品国产| 成年人黄色毛片网站| 久久久久久人人人人人| 成人国语在线视频| 丰满的人妻完整版| 精品国内亚洲2022精品成人| 叶爱在线成人免费视频播放| 一卡2卡三卡四卡精品乱码亚洲| 最好的美女福利视频网| 午夜福利成人在线免费观看| 免费在线观看完整版高清| 中文字幕av电影在线播放| 99国产综合亚洲精品| 色老头精品视频在线观看| 久久精品成人免费网站| 国产高清视频在线播放一区| 女人被狂操c到高潮| 在线观看一区二区三区| 黑人欧美特级aaaaaa片| 色综合站精品国产| 中文字幕精品免费在线观看视频| 欧美成人一区二区免费高清观看 | 亚洲精品国产一区二区精华液| 午夜福利一区二区在线看| 黑人欧美特级aaaaaa片| 国产精品爽爽va在线观看网站 | videosex国产| 超碰成人久久| 国产成人影院久久av| 悠悠久久av| 99久久精品国产亚洲精品| www日本黄色视频网| 日韩欧美国产在线观看| 亚洲中文av在线| 久久精品亚洲精品国产色婷小说| 夜夜爽天天搞| 国产高清视频在线播放一区| 久久精品国产99精品国产亚洲性色| 午夜福利一区二区在线看| 国内精品久久久久精免费| or卡值多少钱| 美女国产高潮福利片在线看| 香蕉国产在线看| 国产精品影院久久| 亚洲avbb在线观看| 美女高潮喷水抽搐中文字幕| 91九色精品人成在线观看| 69av精品久久久久久| 欧美三级亚洲精品| 日韩一卡2卡3卡4卡2021年| 男女下面进入的视频免费午夜 | 午夜a级毛片| 人人澡人人妻人| 成人亚洲精品一区在线观看| 日韩精品中文字幕看吧| 1024手机看黄色片| 亚洲精品在线美女| 国产av又大| 亚洲电影在线观看av| 十八禁网站免费在线| tocl精华| 999久久久精品免费观看国产| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| 亚洲熟妇中文字幕五十中出| 日本一本二区三区精品| 色综合亚洲欧美另类图片| 中文字幕最新亚洲高清| 久久99热这里只有精品18| 黄色丝袜av网址大全| 18禁国产床啪视频网站| 中文在线观看免费www的网站 | 午夜福利在线在线| 在线观看日韩欧美| 老司机在亚洲福利影院| cao死你这个sao货| 欧美日韩精品网址| av视频在线观看入口| 在线十欧美十亚洲十日本专区| 最近最新中文字幕大全免费视频| 窝窝影院91人妻| 国产精品一区二区三区四区久久 | 麻豆成人av在线观看| 麻豆av在线久日| 亚洲国产高清在线一区二区三 | 久久久久免费精品人妻一区二区 | 国产成人欧美在线观看| www日本黄色视频网| 中出人妻视频一区二区| 一级a爱片免费观看的视频| 美女扒开内裤让男人捅视频| 欧美在线黄色| 午夜成年电影在线免费观看| 久久天堂一区二区三区四区| 国产高清videossex| 久久精品影院6| 亚洲av中文字字幕乱码综合 | 看片在线看免费视频| 哪里可以看免费的av片| 麻豆一二三区av精品| 91成年电影在线观看| 久久久久久久久免费视频了| 丁香欧美五月| 久久久国产成人精品二区| 日韩欧美一区视频在线观看| 国内久久婷婷六月综合欲色啪| 男女床上黄色一级片免费看| 一级毛片精品| 国产一区二区三区在线臀色熟女| 日本熟妇午夜| 麻豆av在线久日| 亚洲国产欧美日韩在线播放| 久久这里只有精品19| av天堂在线播放| 国产97色在线日韩免费| 欧美在线黄色| 国产片内射在线| 国产蜜桃级精品一区二区三区| 黄色丝袜av网址大全| 欧美午夜高清在线| 午夜久久久久精精品| 在线免费观看的www视频| 国产99久久九九免费精品| 国产精品爽爽va在线观看网站 | 午夜a级毛片| www.熟女人妻精品国产| 欧美黑人欧美精品刺激| 日本一本二区三区精品| 女生性感内裤真人,穿戴方法视频| netflix在线观看网站| 叶爱在线成人免费视频播放| 精品福利观看| 禁无遮挡网站| 成人国产一区最新在线观看| 午夜激情福利司机影院| 婷婷精品国产亚洲av在线| av片东京热男人的天堂| 午夜福利在线在线| 国产97色在线日韩免费| 国产在线观看jvid| 色综合欧美亚洲国产小说| 日韩 欧美 亚洲 中文字幕| 99热这里只有精品一区 | 又大又爽又粗| 亚洲 欧美 日韩 在线 免费| 亚洲成av人片免费观看| 身体一侧抽搐| 日韩有码中文字幕| 女人高潮潮喷娇喘18禁视频| 亚洲精品国产一区二区精华液| 国产免费av片在线观看野外av| 亚洲天堂国产精品一区在线| 免费高清视频大片| 在线国产一区二区在线| 视频区欧美日本亚洲| 中文字幕精品免费在线观看视频| 亚洲aⅴ乱码一区二区在线播放 | 身体一侧抽搐| 久久这里只有精品19| 国产精品香港三级国产av潘金莲| 制服人妻中文乱码| 黄色视频不卡| 欧美性长视频在线观看| 男女视频在线观看网站免费 | 国语自产精品视频在线第100页| 国产成人影院久久av| 每晚都被弄得嗷嗷叫到高潮| 日韩一卡2卡3卡4卡2021年| 国产真人三级小视频在线观看| 正在播放国产对白刺激| 99在线人妻在线中文字幕| 久久精品国产综合久久久| 免费人成视频x8x8入口观看| 久久久久久免费高清国产稀缺| 午夜激情福利司机影院| 亚洲片人在线观看| 黑人欧美特级aaaaaa片|