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

    非勻速算法在橋梁動(dòng)態(tài)稱重中的應(yīng)用研究

    2020-11-18 10:40:58安家禾
    公路工程 2020年5期
    關(guān)鍵詞:計(jì)算誤差軸重單軸

    卓 穎, 安家禾, 張 斌, 朱 平, 趙 華

    (湖南大學(xué) 風(fēng)工程與橋梁工程湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 長沙 410082)

    0 引言

    車輛超載會(huì)加速道路基礎(chǔ)設(shè)施的老化,尤其會(huì)顯著縮短橋梁的使用壽命[1-2]。在一些極端情況下,超載車輛的重量甚至可能超過橋梁的承載能力,直接導(dǎo)致橋梁倒塌。因此,治理車輛超載問題對(duì)保護(hù)和維護(hù)現(xiàn)代運(yùn)輸系統(tǒng)變得越來越重要。

    為了獲取移動(dòng)車輛的軸重信息,通??刹捎渺o態(tài)稱重和動(dòng)態(tài)稱重(Weigh-In-Motion: WIM)這兩種技術(shù)。靜態(tài)稱重能夠精確獲得車輛軸重信息,但是經(jīng)濟(jì)和時(shí)間成本很高,并且會(huì)妨礙正常交通,因此是不切實(shí)際的。為了克服靜態(tài)稱重的限制,自20世紀(jì)60年代以來,基于WIM的路面稱重系統(tǒng)(Pavement Weigh-In-Motion: PWIM)發(fā)展了起來[3]。PWIM系統(tǒng)使用安裝在道路上的設(shè)備在正常交通條件下對(duì)公路車輛進(jìn)行稱重。

    Moses于1979年首先提出了橋梁動(dòng)態(tài)稱重(Bridge Weigh-In-Motion: BWIM)的概念。與PWIM系統(tǒng)不同,BWIM系統(tǒng)將一座橋作為稱重秤來計(jì)算車輛的重量。BWIM系統(tǒng)具有耐用、易安裝和維護(hù)、軸重識(shí)別精度高的優(yōu)點(diǎn)[4]。這些優(yōu)點(diǎn)使BWIM系統(tǒng)成為治理車輛超載問題一個(gè)有效工具。

    橋梁動(dòng)態(tài)稱重(BWIM)系統(tǒng)計(jì)算軸重的方法通??煞譃閮纱箢?,一類是靜態(tài)算法,主要計(jì)算移動(dòng)車輛的靜態(tài)軸重,以Moses影響線算法為主要代表;另一類是動(dòng)態(tài)算法,主要計(jì)算軸重的時(shí)程曲線,以移動(dòng)荷載識(shí)別(Moving Force Identification: MFI)算法為代表。

    Moses算法通過最小化實(shí)測(cè)橋梁響應(yīng)與理論橋梁響應(yīng)之間的差值來計(jì)算車輛軸重,而理論橋梁響應(yīng)是利用影響線計(jì)算出來的[4]。

    Moses算法現(xiàn)已被用來建立現(xiàn)代商業(yè)BWIM系統(tǒng)基本框架。20世紀(jì)80年代,PETERS[5]在澳大利亞開發(fā)了Axway系統(tǒng)。后來,PETERS[6]開發(fā)了一種更有效的系統(tǒng),稱為Culway,它使用一個(gè)涵洞作為稱重秤。之所以用涵洞而不是橋梁來進(jìn)行稱重,是因?yàn)檐囕v和涵洞之間相互作用所產(chǎn)生的動(dòng)力效應(yīng)可以更快地被周圍的土壤所抑制。在歐洲,COST 323行動(dòng)和WAVE(Weighing in Motion of Axles and Vehicles for Europe)項(xiàng)目在1990年代末實(shí)施[7-8]。這些項(xiàng)目極大地提高了BWIM技術(shù)的準(zhǔn)確性,并促使了著名商業(yè)BWIM系統(tǒng)——SiWIM系統(tǒng)的開發(fā)。SiWIM系統(tǒng)由ZAG[9]在2005年基于Moses理論所研發(fā)出來的。該系統(tǒng)采用一種修正的彎矩影響線,使得單軸和總軸重計(jì)算更為精確。該系統(tǒng)在板橋上的軸重計(jì)算結(jié)果精度是可以接受的。

    MCNULTY[10]和O′BRIEN[10]在2003年提出只有標(biāo)定影響線盡可能地與實(shí)際一致,BWIM算法才能得到理想的精度。ZHAO[11]等在2015年提出了一個(gè)用于軸重識(shí)別的修正Moses算法,也證實(shí)了獲得與實(shí)際一致的標(biāo)定影響線對(duì)于BWIM系統(tǒng)的成功應(yīng)用十分重要。

    LANSDELL A[12]等在2017年提出了一種考慮車輛以分勻速行駛的情況的算法,并運(yùn)用實(shí)橋試驗(yàn)證實(shí)了該算法的有效性。

    本文的研究主要基于Moses的影響線算法,針對(duì)其局限性,提出了適用范圍更廣的非勻速算法,并利用MATLAB軟件進(jìn)行數(shù)值模擬,驗(yàn)證了其有效性。

    1 Moses影響線算法

    Moses算法主要由兩部分組成。第一部分需要對(duì)測(cè)試橋梁進(jìn)行標(biāo)定試驗(yàn),用已知軸重和軸距信息的車輛去標(biāo)定橋梁,得到實(shí)際的橋梁影響線。第二部分是讓實(shí)際車輛在橋上行駛,獲得實(shí)測(cè)應(yīng)變數(shù)據(jù),再用標(biāo)定試驗(yàn)得到的計(jì)算影響線來進(jìn)行計(jì)算,從而獲得車輛的軸重和軸距信息。

    1.1 影響線計(jì)算公式推導(dǎo)

    當(dāng)車輛通過橋梁時(shí),跨中的縱橋向總彎矩,可以表示成時(shí)間的函數(shù),并且等于每根梁的彎矩之和。為了簡化起見,將整座橋視為一個(gè)整體,且每根梁具有相同的彈性模量E與截面模量Z。故橋梁某一截面的總彎矩Mk,在時(shí)間步k可以由式(1)表示。

    (1)

    在整個(gè)標(biāo)定期間,假定標(biāo)定車是勻速行駛。車速由布置在橋面板下的兩排FAD(Free of Axle Detector: FAD)傳感器所采集的數(shù)據(jù)得到。已知一輛標(biāo)定車的N個(gè)軸的軸重P1,P2,…,PN,當(dāng)車過橋時(shí),在時(shí)間步k時(shí),由車產(chǎn)生的跨中彎矩Mk即為:

    (2)

    (3)

    (4)

    基于最小二乘法,理論應(yīng)變與實(shí)測(cè)應(yīng)變的誤差函數(shù)可定義為:

    (5)

    為了使E最小,對(duì)E中IR項(xiàng)求偏導(dǎo)并令其偏導(dǎo)數(shù)為0,可得:

    (6)

    測(cè)量總時(shí)間步為K,則從時(shí)間步R=1到K,共可得到K-CN個(gè)方程,將這些方程聯(lián)立便可以解出各個(gè)位置的影響線豎標(biāo)值。

    1.2 軸重計(jì)算公式推導(dǎo)

    (7)

    (8)

    則對(duì)應(yīng)時(shí)刻的理論應(yīng)變?yōu)椋?/p>

    (9)

    按照1.1小節(jié)的計(jì)算公式,可由標(biāo)定試驗(yàn)數(shù)據(jù)計(jì)算出影響線矩陣[I]。由式(9),時(shí)間步k從1到K有K個(gè)理論應(yīng)變:

    (10)

    其中,K是采集的應(yīng)變數(shù)目,式(10)可簡化為:

    [εt]K×1=[IL]K×N[P]N×1

    (11)

    其中,[εt]是理論應(yīng)變的列向量;[IL]是影響線縱坐標(biāo)的矩陣;[P]是待求的軸重列向量。

    參考式(5),構(gòu)造應(yīng)變誤差函數(shù),

    (12)

    將式(12)寫成矩陣形式,有:

    (13)

    將式(11)代入,對(duì)E中[P]求偏導(dǎo)并令其偏導(dǎo)數(shù)為0,可得:

    (14)

    從式(14)就可以求出各軸的軸重,將其求和便可得到總軸重GVW。

    (15)

    1.3 算法的局限性

    盡管Moses理論在BWIM中已經(jīng)應(yīng)用得很廣泛,但是該算法仍然存在缺點(diǎn)。Moses算法中假定了車輛以勻速行駛,但車輛在實(shí)際行駛時(shí),其車速可能會(huì)發(fā)生變化,速度變化也會(huì)引起軸重計(jì)算誤差。當(dāng)車速發(fā)生較大的變化時(shí),Moses算法便不適用了。因此,引入非勻速算法就顯得很重要。

    回顧Moses算法的計(jì)算流程,由式(3)可知Ci為與Di相應(yīng)的采集點(diǎn)數(shù),且與車速相關(guān)。當(dāng)車輛以勻速行駛時(shí),車輛位移隨著時(shí)間均勻變化,影響線可以分成相等的若干段,每一段的距離便是v/f,而Ci保持不變。而當(dāng)車速在不斷變化時(shí),車輛位移與時(shí)間的關(guān)系發(fā)生了變化,影響線矩陣也需要重新進(jìn)行分段處理才能與車的位置對(duì)應(yīng),Ci也會(huì)隨之變化。在軸重計(jì)算中,式(11)中的[IL]矩陣與車軸的位置有關(guān),而軸距Di在該矩陣中是用Ci間接表示的。也就說,當(dāng)車速在不斷變化時(shí),Di需要根據(jù)不同的車速換算成新的Ci然后才能構(gòu)建[IL]矩陣,使計(jì)算變得十分繁瑣。若Ci與影響線分段不變,而車速在變,則[IL]矩陣中的元素?cái)?shù)值會(huì)與實(shí)際不一致,為軸重計(jì)算帶來誤差。

    2 非勻速算法

    非勻速算法的計(jì)算建立在Moses所提出的BWIM算法上,同樣是基于最小二乘法,通過使跨中理論應(yīng)變與實(shí)測(cè)應(yīng)變的誤差最小來求軸重。在非勻速算法中,不再需要假定車輛以勻速行駛。

    2.1 影響線與軸重的計(jì)算

    從車輛第一個(gè)軸上橋開始測(cè)量,到最后一個(gè)軸離橋停止測(cè)量,車輛行駛的總距離為X=L+DN(L為橋長)。以0.025 m為一個(gè)區(qū)間,可將X分為K=X/0.025個(gè)區(qū)間,各區(qū)間端點(diǎn)編號(hào)依次為0,1……K。

    為了克服Moses算法在處理軸距時(shí)的缺陷,用距離分段代替時(shí)間分段,減少轉(zhuǎn)換步驟,為后續(xù)考慮變速情況提供便利。即將式(3)換成式(16)。

    (16)

    其中,Ci是與Di相對(duì)應(yīng)的分段數(shù)。

    將式(16)代入式(6),此時(shí)R表示車輛第一軸所在的區(qū)間端點(diǎn)號(hào),同樣從端點(diǎn)1到K可以構(gòu)建K-CN個(gè)方程,聯(lián)立之后便可求得影響線豎標(biāo)值。

    至于軸重計(jì)算,關(guān)鍵也在于Ci的變化。將式(16)代入式(8)~式(10),再用最小二乘法使得理論應(yīng)變與實(shí)測(cè)應(yīng)變的差值平方和最小,即可求得車輛軸重。

    該算法其實(shí)是用距離間隔代替時(shí)間間隔。在原始Moses算法中,由于儀器的采集頻率不變,因此時(shí)間間隔是恒定的。而在非勻速算法中,使距離間隔保持不變,然后根據(jù)車速和車輛行駛的距離算出時(shí)間。當(dāng)車速不變時(shí),距離間隔與時(shí)間間隔都是恒定的,與原Moses算法等效。而當(dāng)車速改變時(shí),距離間隔不變,時(shí)間間隔會(huì)變化。這樣處理可以在考慮車速變化時(shí),減少計(jì)算步驟。除此之外更為重要的是,距離,速度與時(shí)間這三者的關(guān)系可以單獨(dú)定義,而不影響軸重計(jì)算中[IL]矩陣的構(gòu)建。

    2.2 速度假定與軸距計(jì)算

    Moses算法假定車輛以勻速行駛,其計(jì)算車速是通過先求出每一根軸經(jīng)過兩排FAD傳感器的平均速度,再將這些速度平均所得到。為了充分利用試驗(yàn)數(shù)據(jù),現(xiàn)引入新的速度計(jì)算方法。一輛N軸車,每一個(gè)軸經(jīng)過一排FAD傳感器都會(huì)產(chǎn)生一個(gè)信號(hào),而在橋梁上裝有兩排FAD傳感器,所以一輛N軸車行駛一趟一共可以產(chǎn)生2N個(gè)信號(hào)。2N個(gè)信號(hào)就意味著有2N個(gè)時(shí)間點(diǎn),2N-1段距離,那么便可以得到車輛在2N-1個(gè)區(qū)間的平均速度。本文假定車輛在這2N-1個(gè)區(qū)間內(nèi)分別以各個(gè)區(qū)間的平均速度行駛。除此之外還可以假定車輛以其他方式運(yùn)動(dòng),采用不同的假定會(huì)有不同的計(jì)算結(jié)果。本文僅討論車輛以分段勻速的方式行駛。

    在實(shí)際的軸重計(jì)算中,移動(dòng)車輛的軸距是未知的,需要利用傳感器的信號(hào)先對(duì)其進(jìn)行計(jì)算,這樣才能順利構(gòu)建軸重計(jì)算矩陣。

    假定車輛的相鄰2個(gè)軸在通過兩排FAD傳感器時(shí)做勻加速運(yùn)動(dòng),前軸經(jīng)過2個(gè)FAD傳感器的時(shí)間分別記為t1,1,t2,1,對(duì)應(yīng)的車輛瞬時(shí)速度記為V1,V3;后軸經(jīng)過2個(gè)FAD傳感器的時(shí)間分別記為t1,2,t2,2,對(duì)應(yīng)的車輛瞬時(shí)速度記為V2,V4;2個(gè)車軸的距離記為D,2排傳感器之間的距離記為S。根據(jù)D與S的大小關(guān)系,可能出現(xiàn)2種情況,如圖1和圖2所示。

    圖1 軸距計(jì)算圖(S>D)Figure 1 Wheelbase calculation (S>D)

    圖2 軸距計(jì)算圖(D>S)Figure 2 Wheelbase calculation (D>S)

    由圖1可知,當(dāng)S>D時(shí),車輛后軸經(jīng)過第1排FAD傳感器時(shí),前軸還未到達(dá)第2排傳感器,此時(shí)前軸與第2排FAD傳感器的距離為S-D;由圖2可知,當(dāng)D>S時(shí),車輛前軸經(jīng)過第2排傳感器時(shí),后軸還未到達(dá)第1排FAD傳感器,此時(shí)后軸距離第1排FAD傳感器的距離為D-S。

    前軸在2排FAD傳感器之間的平均速度為:

    (17)

    Δt1=t2,1-t1,1

    (18)

    后軸在兩排FAD傳感器之間的平均速度為:

    (19)

    Δt2=t2,2-t1,2

    (20)

    由此,可以計(jì)算車輛的加速度a。

    (21)

    (22)

    (23)

    故可求出V1和V4:

    (24)

    (25)

    由第1排FAD傳感器算出來的軸距D1為:

    (26)

    由第2排FAD傳感器算出來的軸距D2為:

    (27)

    相鄰兩軸的軸距D取D1和D2的平均值,

    (28)

    當(dāng)車輛做勻速運(yùn)動(dòng)時(shí),即加速度為0時(shí),上述計(jì)算公式依然適用,且與原來的計(jì)算方法等價(jià)。

    3 數(shù)值模擬

    數(shù)值模擬不考慮車橋的振動(dòng)及相互作用,僅考慮數(shù)值上的計(jì)算。計(jì)算橋梁為一跨長度為20 m的簡支梁,傳感器布置如圖3所示。應(yīng)變計(jì)始終位于橋梁跨中,用于采集橋梁跨中的應(yīng)變數(shù)據(jù),在計(jì)算中為了減少轉(zhuǎn)換步驟,實(shí)際采用的是跨中彎矩,標(biāo)明應(yīng)變計(jì)主要是示意其位置。2排FAD傳感器之間距離為S,F(xiàn)AD與最近橋梁支點(diǎn)的距離為S1,F(xiàn)AD2與最近橋梁支點(diǎn)的距離為S2,且S1+S+S2=20 m。車輛行進(jìn)方向與x軸正方向一致。

    圖3 傳感器布置圖Figure 3 Sensor layout

    計(jì)算選用了3軸車和5軸車這2種車型,因?yàn)檫@兩者都同時(shí)具有單軸和組軸,比較具有代表性。車輛的軸距及軸重信息如表1所示。

    表1 車輛軸距軸重信息Table 1 Vehicle wheelbase and axle load information軸數(shù)軸距/m軸重/kNA1-A2A2-A3A3-A4A4-A5總長A1A2A3A4A5總軸重34.00 1.25 ——5.25 60 120 125 ——305 54.30 1.35 11.20 1.30 18.15 50 70 70 80 80 350

    數(shù)據(jù)采樣頻率f為512 Hz,利用MATLAB模擬車輛過橋過程來產(chǎn)生計(jì)算數(shù)據(jù)。計(jì)算所用的影響線為橋梁的理論影響線。

    3.1 勻速運(yùn)動(dòng)

    在討論車輛做勻加速運(yùn)動(dòng)前,首先要對(duì)計(jì)算數(shù)據(jù)段的選取做一個(gè)說明。一般來說,橋梁動(dòng)態(tài)稱重試驗(yàn)會(huì)布置2排FAD傳感器,用于檢測(cè)車軸位置和計(jì)算車速。所以車輛在這2排FAD傳感器之外的速度是未知的,也就是說車輛的運(yùn)動(dòng)可以有無限種可能,因此傳感器區(qū)間之外的數(shù)據(jù)不宜用來進(jìn)行計(jì)算。計(jì)算所采用的數(shù)據(jù)區(qū)間,是從車輛第1軸到達(dá)FAD1開始,到車輛最后1軸離開FAD2為止。

    車輛初始速度為10 m/s,從車輛第1軸到達(dá)FAD1開始加速,加速度分別取0.5、1.0、1.5、2.0 、2.5、3.0 m/s2。用2種方法進(jìn)行計(jì)算,第1種,假定車輛全程做勻速運(yùn)動(dòng),用Moses算法求解,用“全程勻速”表示;第2種,假定車輛做分段勻速運(yùn)動(dòng),用非勻速算法求解,用“分段勻速”表示。

    3軸車的軸重計(jì)算結(jié)果如圖4~圖6所示。由圖可知,2種方法計(jì)算的單軸和組軸誤差都非常小,保持在1%以內(nèi),但是單軸的誤差都較大,誤差最大時(shí)都超過了6%。

    圖4 3軸車總軸重計(jì)算誤差折線圖Figure 4 The calculation error line chart of the gross vehicle weight of 3-axle vehicle

    圖5 3軸車單軸重(A1)計(jì)算誤差折線圖Figure 5 The calculation error line chart of the single axle (A1) weight of 3-axle vehicle

    圖6 3軸車組軸重(A2+A3)計(jì)算誤差折線圖Figure 6 The calculation error line chart of the group of axles (A2+A3) weight of 3-axle vehicle

    5軸車的軸重計(jì)算結(jié)果如圖7~圖10所示。由圖可知,就總軸重而言,“全程勻速”的計(jì)算結(jié)果更好,其最大誤差比“分段勻速”的最大誤差小了近1%。對(duì)于單軸A1,“分段勻速”的計(jì)算結(jié)果明顯更好,其最大誤差基本在10%以內(nèi),而“全程勻速”的最大誤差約20%,是“分段勻速”的2倍。對(duì)于A2和A3組成的組軸,2種方法的計(jì)算結(jié)果相近,兩者最大誤差均在2%以內(nèi)。對(duì)于A4和A5組成的組軸,“分段勻速”最大誤差約為2%,而“全程勻速”最大誤差約為4%,兩者有2%的差距。

    圖7 5軸車總軸重計(jì)算誤差折線圖Figure 7 The calculation error line chart of the gross vehicle weight of 5-axle vehicle

    圖8 5軸車單軸重(A1)計(jì)算誤差折線圖Figure 8 The calculation error line chart of the single axle (A1) weight of 5-axle vehicle

    圖9 5軸車組軸重(A2+A3)計(jì)算誤差折線圖Figure 9 The calculation error line chart of the group of axles (A2+A3) weight of 5-axle vehicle

    圖10 5軸車組軸重(A4+A5)計(jì)算誤差折線圖Figure 10 The calculation error line chart of the group of axles (A4+A5) weight of 5-axle vehicle

    綜上,“分段勻速”在單軸和組軸的軸重計(jì)算上具有較高的精度,總軸重的計(jì)算精度略低,但也在可接受范圍之內(nèi);“全程勻速”的總軸重計(jì)算精度很高,組軸精度略低,但也在可接受范圍之內(nèi),而單軸的計(jì)算很不穩(wěn)定,可能出現(xiàn)較大的誤差。

    3.2 隨機(jī)變速運(yùn)動(dòng)

    車輛在實(shí)際運(yùn)動(dòng)中可能多次變速,為了模擬更為真實(shí)的車速情況,現(xiàn)在假設(shè)車輛以10 m/s的速度行駛,當(dāng)車輛第1軸到達(dá)FAD1時(shí)以隨機(jī)加速度a1做勻加速運(yùn)動(dòng),直到到達(dá)FAD2,然后再以隨機(jī)加速度a2做勻加速運(yùn)動(dòng),直到最后一軸經(jīng)過FAD2。加速度a1和a2由MATLAB的隨機(jī)函數(shù)生成,兩者均在0~3內(nèi)取值,單位為m/s2。車輛每行駛1趟,便重新生成一次加速度,共生成100趟。計(jì)算采用“全程勻速”和“分段勻速”2種方法來進(jìn)行對(duì)比。計(jì)算選用車型仍為3.1小節(jié)中的3軸車和5軸車,且每1種車型的計(jì)算原始數(shù)據(jù)均相同。

    3軸車的計(jì)算結(jié)果匯總?cè)绫?所示,總軸重,單軸重和組軸重的誤差箱線圖如圖11~圖13所示。

    表2 隨機(jī)變速3軸車軸重計(jì)算結(jié)果比較Table 2 Comparison of axle weight calculation results of 3-axle vehicle with random variable speed項(xiàng)目A1A2A3總軸重單軸組軸數(shù)量100100100100100100全程勻速均值3.16-17.5816.060.293.16-0.41標(biāo)準(zhǔn)差2.038.027.370.192.030.28分段勻速均值-3.41-0.601.29-0.38-3.410.36標(biāo)準(zhǔn)差2.202.972.680.232.200.28注: 均值與標(biāo)準(zhǔn)差的數(shù)據(jù)為百分?jǐn)?shù)分子; “單軸”即為“A1”; “組軸”為“A2+A3”。

    由表2所示,2種方法在總軸重,單軸A1和組軸上的計(jì)算精度非常接近,無明顯區(qū)別。但是組軸中的單軸A2和A3的計(jì)算精度有明顯差別,“全程勻速”的計(jì)算誤差非常大,誤差均值超過了15%,而且標(biāo)準(zhǔn)差也在一個(gè)較高水平,而“分段勻速”的誤差均值在1.5%以內(nèi),標(biāo)準(zhǔn)差也很小。這說明“分段勻速”的軸重計(jì)算更為穩(wěn)定一些。

    由圖11~圖13可知,2種方法在總軸重,單軸A1和組軸上的計(jì)算精度非常接近,無明顯區(qū)別,與表2中的數(shù)據(jù)一致。

    圖11 隨機(jī)變速3軸車總軸重計(jì)算誤差箱線圖Figure 11 Calculation error boxplot of the gross vehicle weight of 3-axle vehicle with random variable speed

    圖12 隨機(jī)變速3軸車單軸重(A1)計(jì)算誤差箱線圖Figure 12 Calculation error boxplot of the single axle (A1) weight of 3-axle vehicle with random variable speed

    圖13 隨機(jī)變速3軸車組軸重計(jì)算誤差箱線圖Figure 13 Calculation error boxplot of the group of axles weight of 3-axle vehicle with random variable speed

    5軸車的計(jì)算結(jié)果匯總?cè)绫?所示,總軸重、單軸重和組軸重的誤差箱線圖如圖14~圖16所示。

    圖14 隨機(jī)變速5軸車總軸重計(jì)算誤差箱線圖Figure 14 Calculation error boxplot of the gross vehicle weight of 5-axle vehicle with random variable speed

    圖15 隨機(jī)變速5軸車單軸重(A1)計(jì)算誤差箱線圖Figure 15 Calculation error boxplot of the single axle (A1) weight of 5-axle vehicle with random variable speed

    圖16 隨機(jī)變速5軸車組軸重計(jì)算誤差箱線圖Figure 16 Calculation error boxplot of the group of axles weight of 5-axle vehicle with random variable speed

    由表3所示,“全程勻速”在總軸重計(jì)算上有略微的優(yōu)勢(shì)。除總軸重這一項(xiàng)外,“全程勻速”的標(biāo)準(zhǔn)差均高于“分段勻速”,說明“分段勻速”的計(jì)算結(jié)果更穩(wěn)定。兩者的組軸計(jì)算結(jié)果相近,但是“分段勻速”的單軸計(jì)算精度要顯著高于“全程勻速”。

    表3 隨機(jī)變速5軸車軸重計(jì)算結(jié)果比較Table 3 Comparison of axle weight calculation results of 5-axle vehicle with random variable speed項(xiàng)目A1A2A3A4A5總軸重單軸組軸數(shù)量100100100100100100100200全程勻速均值10.38-29.4430.98-55.0050.670.8010.38-0.70標(biāo)準(zhǔn)差6.2516.3317.1026.4724.410.566.251.67分段勻速均值-5.4225.43-26.91-36.5034.43-1.54-5.42-0.89標(biāo)準(zhǔn)差2.1112.9214.1815.4114.160.782.110.68注: 均值與標(biāo)準(zhǔn)差的數(shù)據(jù)為百分?jǐn)?shù)分子; “單軸”即為“A1”; “組軸”為“A2+A3”和“A3+A4”。

    由圖14~圖16可知,“全程勻速”的總軸重計(jì)算精度要比“分段勻速”略高,但兩者的總軸重誤差均小于3.0%,在可接受范圍內(nèi)?!胺侄蝿蛩佟钡膯屋S重計(jì)算誤差基本可以控制在10%以內(nèi),要顯著高于“全程勻速”。對(duì)于組軸重的計(jì)算精度,綜上,“全程勻速”的總軸重計(jì)算結(jié)果更好,而“分段勻速”的單軸計(jì)算結(jié)果更好,兩者的組軸計(jì)算結(jié)果相近??紤]到兩者的總軸重和組軸的計(jì)算精度相差不大,而“分段勻速”在單軸精度上具有明顯優(yōu)勢(shì),因此可認(rèn)為在車輛以非勻速行駛時(shí),非勻速算法要比原始Moses算法更好。

    4 結(jié)論

    a.Moses算法只考慮車輛以勻速行使的情況,當(dāng)車輛做非勻速運(yùn)動(dòng)時(shí),Moses算法難以將車速的變化考慮進(jìn)去,而非勻速算法可以有效解決這個(gè)問題。

    b.數(shù)值模擬的計(jì)算結(jié)果表明,當(dāng)考慮車速變化時(shí),非勻速算法的計(jì)算精度整體在一個(gè)較高的水平,尤其是在單軸計(jì)算上有顯著優(yōu)勢(shì);Moses算法在總軸重計(jì)算上具有一定的優(yōu)勢(shì),組軸精度在可接受范圍之內(nèi),但是其單軸的精度很差,計(jì)算結(jié)果非常不穩(wěn)定。

    c.在非勻速算法中,通過對(duì)FAD信號(hào)的合理處理,獲得了車輛在各個(gè)區(qū)間的平均速度,并假定車輛在各個(gè)區(qū)間以其平均速度行駛。數(shù)值模擬分析證實(shí)了該假設(shè)的有效性。

    猜你喜歡
    計(jì)算誤差軸重單軸
    炭黑填充天然橡膠超彈性本構(gòu)方程的適用性分析
    單軸壓縮條件下巖石峰后第Ⅱ種類型應(yīng)力——應(yīng)變曲線的新解釋
    20t軸重米軌轉(zhuǎn)向架裝用搖枕的研制
    CFRP-鋼復(fù)合板的單軸拉伸力學(xué)性能
    水尺計(jì)重中密度測(cè)量與計(jì)算誤差分析及相關(guān)問題的思考
    水尺計(jì)重中密度測(cè)量與計(jì)算誤差分析及相關(guān)問題的思考
    單軸應(yīng)變Si NMOS電流模型研究
    斜單軸跟蹤式光伏組件的安裝傾角優(yōu)化設(shè)計(jì)
    32.5t軸重貨車作用下重載鐵路軌道的合理剛度
    強(qiáng)度折減法中折減參數(shù)對(duì)邊坡穩(wěn)定性計(jì)算誤差影響研究
    狂野欧美白嫩少妇大欣赏| 精品欧美国产一区二区三| 欧美丝袜亚洲另类 | 国产在视频线在精品| 内射极品少妇av片p| 一区二区三区高清视频在线| 深夜精品福利| 在现免费观看毛片| 中文字幕久久专区| ponron亚洲| 久久久久久久久大av| 免费看美女性在线毛片视频| 男人狂女人下面高潮的视频| 国产亚洲精品综合一区在线观看| 好男人电影高清在线观看| 亚洲av免费在线观看| 亚洲精品影视一区二区三区av| 又爽又黄无遮挡网站| 欧美最黄视频在线播放免费| 动漫黄色视频在线观看| 我的女老师完整版在线观看| 亚洲av免费高清在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 女生性感内裤真人,穿戴方法视频| 禁无遮挡网站| 啦啦啦观看免费观看视频高清| 国产亚洲精品久久久久久毛片| av欧美777| 精品免费久久久久久久清纯| 亚洲美女搞黄在线观看 | 欧美高清成人免费视频www| 欧美一区二区精品小视频在线| 亚洲18禁久久av| 十八禁人妻一区二区| 亚洲精品乱码久久久v下载方式| 欧美日韩瑟瑟在线播放| 深夜a级毛片| 欧美一区二区亚洲| 如何舔出高潮| 国产精品1区2区在线观看.| 中文字幕人成人乱码亚洲影| 亚洲五月天丁香| 一卡2卡三卡四卡精品乱码亚洲| 搞女人的毛片| 久久精品久久久久久噜噜老黄 | 亚洲avbb在线观看| 在线a可以看的网站| 搡女人真爽免费视频火全软件 | 国产伦精品一区二区三区视频9| 日本a在线网址| 国产成人欧美在线观看| 成年版毛片免费区| 色综合亚洲欧美另类图片| 亚洲aⅴ乱码一区二区在线播放| 色视频www国产| 色5月婷婷丁香| 国产三级中文精品| 亚洲性夜色夜夜综合| x7x7x7水蜜桃| 黄色配什么色好看| 精品国产三级普通话版| 一边摸一边抽搐一进一小说| 级片在线观看| 悠悠久久av| 99在线人妻在线中文字幕| 亚洲av免费高清在线观看| 欧美一级a爱片免费观看看| 不卡一级毛片| 一个人免费在线观看的高清视频| 老司机福利观看| 国内精品一区二区在线观看| 亚洲国产精品999在线| 免费观看人在逋| 欧美丝袜亚洲另类 | 亚洲成人中文字幕在线播放| 男女视频在线观看网站免费| 观看免费一级毛片| 嫩草影视91久久| 精品久久久久久久久av| 一个人观看的视频www高清免费观看| 99热这里只有精品一区| 国产av一区在线观看免费| 久久久久久久亚洲中文字幕 | 麻豆国产97在线/欧美| 精品熟女少妇八av免费久了| 91九色精品人成在线观看| 免费看a级黄色片| 午夜亚洲福利在线播放| 日本在线视频免费播放| 悠悠久久av| 丰满乱子伦码专区| 如何舔出高潮| 亚洲经典国产精华液单 | 免费在线观看日本一区| av中文乱码字幕在线| 俄罗斯特黄特色一大片| 亚洲精品乱码久久久v下载方式| 精品不卡国产一区二区三区| 久久久久久国产a免费观看| 国产中年淑女户外野战色| 亚洲午夜理论影院| 草草在线视频免费看| 日韩欧美精品免费久久 | 亚洲最大成人av| 在线看三级毛片| 亚洲va日本ⅴa欧美va伊人久久| 在线观看一区二区三区| 久久午夜亚洲精品久久| 色5月婷婷丁香| 午夜福利欧美成人| 国产成人影院久久av| 真实男女啪啪啪动态图| av中文乱码字幕在线| 精品免费久久久久久久清纯| 国产 一区 欧美 日韩| 欧美日韩黄片免| 有码 亚洲区| 免费在线观看亚洲国产| 天天一区二区日本电影三级| 亚洲av免费在线观看| 日韩有码中文字幕| 人妻制服诱惑在线中文字幕| 88av欧美| 日本黄色片子视频| 全区人妻精品视频| 听说在线观看完整版免费高清| 国产v大片淫在线免费观看| 亚洲专区国产一区二区| 国产三级黄色录像| 99热这里只有精品一区| 中文字幕av成人在线电影| 偷拍熟女少妇极品色| 91麻豆av在线| av在线蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 欧美中文日本在线观看视频| 一区二区三区激情视频| 亚洲av一区综合| 欧美另类亚洲清纯唯美| 女同久久另类99精品国产91| 国产精品自产拍在线观看55亚洲| 直男gayav资源| 成年女人看的毛片在线观看| 少妇裸体淫交视频免费看高清| 色综合欧美亚洲国产小说| 在线观看66精品国产| 天堂√8在线中文| 日本撒尿小便嘘嘘汇集6| 天天一区二区日本电影三级| 有码 亚洲区| 久久久国产成人精品二区| 国产乱人伦免费视频| 久久中文看片网| 亚洲成人久久爱视频| 国产探花极品一区二区| 亚洲av成人av| 亚洲av电影不卡..在线观看| 国产精品综合久久久久久久免费| 一个人观看的视频www高清免费观看| 国产成人a区在线观看| 每晚都被弄得嗷嗷叫到高潮| 69人妻影院| 尤物成人国产欧美一区二区三区| 欧美又色又爽又黄视频| www.www免费av| 国产伦精品一区二区三区四那| 久久人人精品亚洲av| 一区二区三区四区激情视频 | 中文字幕熟女人妻在线| 精品人妻一区二区三区麻豆 | 色综合婷婷激情| 男女下面进入的视频免费午夜| 18禁裸乳无遮挡免费网站照片| 成人欧美大片| 在线观看一区二区三区| 99在线视频只有这里精品首页| 国产午夜精品论理片| 日本 av在线| 男人和女人高潮做爰伦理| 在线国产一区二区在线| 99热只有精品国产| 成年免费大片在线观看| 国产一区二区三区视频了| 亚洲熟妇熟女久久| 亚洲精品乱码久久久v下载方式| 12—13女人毛片做爰片一| 一夜夜www| 乱码一卡2卡4卡精品| 日韩人妻高清精品专区| 国产精品久久视频播放| 成人永久免费在线观看视频| 日韩欧美一区二区三区在线观看| aaaaa片日本免费| 在线观看66精品国产| 久久久久久久午夜电影| 精品人妻一区二区三区麻豆 | 亚洲无线观看免费| 波多野结衣高清作品| 国产亚洲欧美98| 在线观看免费视频日本深夜| 午夜亚洲福利在线播放| 麻豆一二三区av精品| 国内揄拍国产精品人妻在线| 国产成人av教育| 黄色一级大片看看| 露出奶头的视频| 蜜桃久久精品国产亚洲av| 国产爱豆传媒在线观看| 嫩草影院入口| 久久久久免费精品人妻一区二区| 长腿黑丝高跟| 尤物成人国产欧美一区二区三区| 97人妻精品一区二区三区麻豆| 给我免费播放毛片高清在线观看| 久久久久国产精品人妻aⅴ院| 美女高潮喷水抽搐中文字幕| 午夜福利成人在线免费观看| 免费看a级黄色片| 午夜福利视频1000在线观看| 国内精品一区二区在线观看| 最近视频中文字幕2019在线8| 乱码一卡2卡4卡精品| 日韩欧美国产一区二区入口| 美女xxoo啪啪120秒动态图 | 国产一区二区在线av高清观看| 啪啪无遮挡十八禁网站| 免费看日本二区| 18+在线观看网站| 久久欧美精品欧美久久欧美| 天堂动漫精品| 1024手机看黄色片| 日本一本二区三区精品| 两个人的视频大全免费| 九九久久精品国产亚洲av麻豆| 我要看日韩黄色一级片| 国内精品久久久久久久电影| 人人妻人人看人人澡| 日韩免费av在线播放| 丰满乱子伦码专区| 亚洲国产色片| 99久久精品一区二区三区| 又粗又爽又猛毛片免费看| 国产精品三级大全| 国产视频内射| 免费一级毛片在线播放高清视频| 日韩精品青青久久久久久| av视频在线观看入口| 日日摸夜夜添夜夜添小说| av欧美777| 很黄的视频免费| 国产精品久久久久久久电影| 亚洲精品影视一区二区三区av| 国产亚洲精品av在线| 丰满人妻一区二区三区视频av| 欧美丝袜亚洲另类 | www.999成人在线观看| 国产精品久久久久久久久免 | 十八禁人妻一区二区| 国产伦一二天堂av在线观看| 久久久久久久午夜电影| 国产探花在线观看一区二区| 欧美黄色片欧美黄色片| 一本精品99久久精品77| 午夜福利欧美成人| 免费在线观看成人毛片| 国产精品综合久久久久久久免费| 男人和女人高潮做爰伦理| 欧美不卡视频在线免费观看| 国产伦人伦偷精品视频| 91在线观看av| 欧美性猛交黑人性爽| 一本精品99久久精品77| 国产免费av片在线观看野外av| 色视频www国产| 一进一出抽搐gif免费好疼| 最近最新免费中文字幕在线| 人妻制服诱惑在线中文字幕| 免费观看人在逋| 国产伦精品一区二区三区四那| 国产91精品成人一区二区三区| 精品日产1卡2卡| 国产三级中文精品| 国产淫片久久久久久久久 | а√天堂www在线а√下载| 少妇的逼水好多| 亚洲国产欧洲综合997久久,| 久久久久久久亚洲中文字幕 | 99久久精品一区二区三区| 欧美午夜高清在线| 国产av一区在线观看免费| 亚洲成人中文字幕在线播放| 国产精品伦人一区二区| 97碰自拍视频| bbb黄色大片| 成人特级黄色片久久久久久久| 色播亚洲综合网| 国产午夜精品论理片| 久久人人爽人人爽人人片va | 精品久久久久久久人妻蜜臀av| 少妇被粗大猛烈的视频| 1024手机看黄色片| 国产精品人妻久久久久久| 最新中文字幕久久久久| 成熟少妇高潮喷水视频| 3wmmmm亚洲av在线观看| 老司机午夜福利在线观看视频| 中文字幕av成人在线电影| 成年人黄色毛片网站| 久久久久国产精品人妻aⅴ院| 亚洲,欧美精品.| 非洲黑人性xxxx精品又粗又长| 久久国产精品影院| 天天躁日日操中文字幕| 国语自产精品视频在线第100页| av在线天堂中文字幕| 天堂av国产一区二区熟女人妻| 丰满人妻一区二区三区视频av| 亚洲精品乱码久久久v下载方式| 国产真实伦视频高清在线观看 | 国产伦在线观看视频一区| 欧美性猛交黑人性爽| 脱女人内裤的视频| 国产一区二区在线av高清观看| 久久久久久久亚洲中文字幕 | 人人妻人人看人人澡| 九色成人免费人妻av| 国产综合懂色| 美女高潮喷水抽搐中文字幕| 国产成人a区在线观看| 香蕉av资源在线| 欧洲精品卡2卡3卡4卡5卡区| 色噜噜av男人的天堂激情| 国产精品伦人一区二区| 欧美在线黄色| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品sss在线观看| 99热这里只有是精品在线观看 | 日韩欧美在线二视频| 一级a爱片免费观看的视频| 亚洲欧美激情综合另类| 亚洲乱码一区二区免费版| bbb黄色大片| 欧美一区二区精品小视频在线| 日韩亚洲欧美综合| 一进一出抽搐gif免费好疼| 国产精品嫩草影院av在线观看 | 日韩欧美国产一区二区入口| aaaaa片日本免费| 男人舔奶头视频| 国产美女午夜福利| 午夜福利在线在线| 欧美色欧美亚洲另类二区| 国产精品综合久久久久久久免费| 69av精品久久久久久| 波多野结衣巨乳人妻| 色吧在线观看| 国产精品永久免费网站| 一进一出抽搐gif免费好疼| 性色avwww在线观看| 搞女人的毛片| 哪里可以看免费的av片| 国产又黄又爽又无遮挡在线| 亚洲欧美日韩无卡精品| 久久九九热精品免费| 欧美精品啪啪一区二区三区| 久久久久久国产a免费观看| 国产aⅴ精品一区二区三区波| 免费看a级黄色片| 真人一进一出gif抽搐免费| 国产精品乱码一区二三区的特点| 国产精品日韩av在线免费观看| 精品人妻一区二区三区麻豆 | 能在线免费观看的黄片| 国产精品98久久久久久宅男小说| 国产精品三级大全| 精品日产1卡2卡| 好男人电影高清在线观看| 十八禁人妻一区二区| 99国产精品一区二区三区| 天堂√8在线中文| 日韩欧美在线二视频| 国产野战对白在线观看| 一本久久中文字幕| 亚洲成人久久性| 国产伦人伦偷精品视频| 最近最新中文字幕大全电影3| 中文字幕av成人在线电影| 一级毛片久久久久久久久女| 久久国产精品人妻蜜桃| av中文乱码字幕在线| 一区二区三区免费毛片| 午夜免费男女啪啪视频观看 | 9191精品国产免费久久| 国产私拍福利视频在线观看| 欧美黄色淫秽网站| 国产精品亚洲一级av第二区| 国产黄色小视频在线观看| 国产精品一区二区免费欧美| 亚洲欧美日韩无卡精品| 最近中文字幕高清免费大全6 | 男人狂女人下面高潮的视频| 久久国产精品影院| 欧美性感艳星| 国产精品电影一区二区三区| 亚洲av电影在线进入| 99久久99久久久精品蜜桃| 国产精品精品国产色婷婷| 两性午夜刺激爽爽歪歪视频在线观看| 毛片一级片免费看久久久久 | 波多野结衣高清无吗| 国内精品久久久久久久电影| 精华霜和精华液先用哪个| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产亚洲av香蕉五月| 亚洲天堂国产精品一区在线| 精品久久久久久久久av| 欧美日韩国产亚洲二区| 欧美性感艳星| 成人精品一区二区免费| 一二三四社区在线视频社区8| 高清在线国产一区| 午夜a级毛片| 亚洲成a人片在线一区二区| 精品人妻偷拍中文字幕| 国产亚洲精品久久久com| av福利片在线观看| 女人十人毛片免费观看3o分钟| 91九色精品人成在线观看| 变态另类丝袜制服| 亚洲人成网站高清观看| 亚洲精品成人久久久久久| aaaaa片日本免费| 色哟哟哟哟哟哟| av视频在线观看入口| 伊人久久精品亚洲午夜| 在线观看午夜福利视频| 麻豆久久精品国产亚洲av| 亚洲人成电影免费在线| 少妇人妻一区二区三区视频| 国产av一区在线观看免费| 天堂动漫精品| 日日摸夜夜添夜夜添av毛片 | 91麻豆av在线| 宅男免费午夜| 波多野结衣巨乳人妻| 超碰av人人做人人爽久久| 99久久无色码亚洲精品果冻| 黄色一级大片看看| 免费观看的影片在线观看| 午夜福利高清视频| 在线看三级毛片| 少妇被粗大猛烈的视频| 亚洲人成网站高清观看| 我要看日韩黄色一级片| 成年版毛片免费区| 亚洲欧美日韩高清专用| 人人妻,人人澡人人爽秒播| 久久草成人影院| 亚洲综合色惰| 麻豆一二三区av精品| 天美传媒精品一区二区| 全区人妻精品视频| 亚洲人成伊人成综合网2020| 最新中文字幕久久久久| 人妻丰满熟妇av一区二区三区| 亚洲av日韩精品久久久久久密| 久久精品国产亚洲av涩爱 | 精品人妻偷拍中文字幕| 婷婷丁香在线五月| 欧美黑人欧美精品刺激| 色综合婷婷激情| 欧美激情在线99| 精品不卡国产一区二区三区| 三级毛片av免费| eeuss影院久久| 午夜久久久久精精品| 国产av不卡久久| 国产成+人综合+亚洲专区| 丝袜美腿在线中文| 麻豆久久精品国产亚洲av| 桃色一区二区三区在线观看| 久久99热这里只有精品18| 免费人成在线观看视频色| 亚洲av一区综合| 日本成人三级电影网站| 少妇人妻精品综合一区二区 | 成人鲁丝片一二三区免费| 久久婷婷人人爽人人干人人爱| 欧美最黄视频在线播放免费| 日本a在线网址| 久久精品91蜜桃| 国产日本99.免费观看| 久久这里只有精品中国| 99久久精品一区二区三区| 69人妻影院| 动漫黄色视频在线观看| 欧美黄色淫秽网站| 精品乱码久久久久久99久播| 日韩欧美在线乱码| 午夜福利在线观看吧| 久久热精品热| 欧美另类亚洲清纯唯美| 国产成人福利小说| 伊人久久精品亚洲午夜| 一本一本综合久久| 97人妻精品一区二区三区麻豆| 久久人妻av系列| 黄色视频,在线免费观看| av天堂中文字幕网| 国产人妻一区二区三区在| 亚洲国产精品sss在线观看| av在线老鸭窝| 我要搜黄色片| 精品一区二区三区视频在线观看免费| 97超级碰碰碰精品色视频在线观看| 国产免费av片在线观看野外av| 看片在线看免费视频| 69人妻影院| 国内精品美女久久久久久| 97超级碰碰碰精品色视频在线观看| 久久久成人免费电影| 中文字幕高清在线视频| 一区二区三区免费毛片| 国产精品影院久久| 欧美国产日韩亚洲一区| 女生性感内裤真人,穿戴方法视频| 黄色配什么色好看| 久99久视频精品免费| 亚洲精品成人久久久久久| 怎么达到女性高潮| 女人十人毛片免费观看3o分钟| 免费在线观看影片大全网站| 成熟少妇高潮喷水视频| 亚洲精品粉嫩美女一区| 亚洲国产精品合色在线| 日韩有码中文字幕| 九九久久精品国产亚洲av麻豆| 国产在线男女| 中亚洲国语对白在线视频| 国产三级中文精品| 午夜久久久久精精品| 99国产综合亚洲精品| 亚洲中文字幕日韩| 国产成人欧美在线观看| 久久精品影院6| 国产精品国产高清国产av| 欧美日韩福利视频一区二区| 国产精品av视频在线免费观看| 国产av麻豆久久久久久久| 久久亚洲真实| 人人妻,人人澡人人爽秒播| 中文字幕人妻熟人妻熟丝袜美| 黄色视频,在线免费观看| 国产精品久久久久久久电影| 久久精品国产99精品国产亚洲性色| 亚洲,欧美精品.| 一个人看视频在线观看www免费| 免费观看的影片在线观看| 国产乱人视频| 性欧美人与动物交配| 午夜a级毛片| 国产aⅴ精品一区二区三区波| 国产成人福利小说| 国语自产精品视频在线第100页| 非洲黑人性xxxx精品又粗又长| 亚洲av.av天堂| 国产午夜精品论理片| 欧美乱妇无乱码| 少妇被粗大猛烈的视频| 少妇的逼水好多| 高清毛片免费观看视频网站| 美女大奶头视频| 国产欧美日韩精品一区二区| 国产精品久久视频播放| 亚洲美女黄片视频| 国产探花在线观看一区二区| 毛片一级片免费看久久久久 | 国产精品久久久久久精品电影| 91九色精品人成在线观看| 亚洲最大成人手机在线| 波野结衣二区三区在线| 欧美日韩福利视频一区二区| 国产免费av片在线观看野外av| 色5月婷婷丁香| 1024手机看黄色片| 草草在线视频免费看| 好看av亚洲va欧美ⅴa在| 真人一进一出gif抽搐免费| 一区福利在线观看| 人人妻人人看人人澡| 午夜福利视频1000在线观看| 日韩中字成人| 国产精品人妻久久久久久| 国产亚洲精品久久久com| 久久精品久久久久久噜噜老黄 | 久久国产乱子免费精品| 国产av麻豆久久久久久久| 97热精品久久久久久| 成人三级黄色视频| 国产在线男女| 久久九九热精品免费| 亚洲片人在线观看| 国产欧美日韩精品一区二区| 日韩高清综合在线| 3wmmmm亚洲av在线观看| 99久久99久久久精品蜜桃| 国产精品嫩草影院av在线观看 | 色综合站精品国产| 亚洲国产高清在线一区二区三| 欧美日韩亚洲国产一区二区在线观看| 亚洲最大成人av| 亚洲精品在线美女| 日本三级黄在线观看| 日韩欧美免费精品| 久久热精品热|