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

    基于智能駕駛電動汽車的多車-橋梁耦合系統(tǒng)動力學(xué)特性研究1)

    2022-10-05 07:21:26張宇王嘉偉李韶華任劍瑩
    力學(xué)學(xué)報(bào) 2022年9期
    關(guān)鍵詞:橋梁振動

    張宇 王嘉偉 李韶華 任劍瑩

    * (石家莊鐵道大學(xué)省部共建交通工程結(jié)構(gòu)力學(xué)行為與系統(tǒng)安全國家重點(diǎn)實(shí)驗(yàn)室,石家莊 050043)

    ? (石家莊鐵道大學(xué)工程力學(xué)系,石家莊 050043)

    引言

    面對全球氣候變暖、能源危機(jī)、空氣污染等現(xiàn)狀,迫于能源和環(huán)保問題的壓力,在新一輪科技革命下,電動汽車及智能駕駛受到了各國的高度重視[1-4].輪轂電機(jī)分布式獨(dú)立驅(qū)動電動汽車取消了傳動系統(tǒng),使其控制靈活、空間大、效率高,是目前電動汽車發(fā)展的大方向,得到了行業(yè)內(nèi)和學(xué)術(shù)界的極大關(guān)注[5-6].由于獨(dú)特的構(gòu)造特性,獨(dú)立驅(qū)動電動汽車的動力學(xué)研究也同樣面臨諸多新挑戰(zhàn): 輪轂電機(jī)直接安裝在車軸上,導(dǎo)致輪胎動載荷增加,影響車輛平順性;電機(jī)電磁力和轉(zhuǎn)矩波動會對車輪造成電機(jī)激勵,繼而加劇車輪振動;輪胎垂向動載會引起輪胎縱向力、橫向力及力矩的波動,影響輪胎地面附著特性.

    針對獨(dú)立驅(qū)動電動汽車的機(jī)電耦合振動特性,學(xué)者們開展了相關(guān)研究.Wang 等[7]研究發(fā)現(xiàn),當(dāng)垂向電磁力直接施加在車輪上,將導(dǎo)致輪胎負(fù)載顯著變化.Zhao 等[8]分析了輪轂電機(jī)電動車在道路影響和電機(jī)振動下,懸架阻尼對垂向特性的影響.Li 等[9]提出了一種主動懸架的多目標(biāo)優(yōu)化控制方法,可有效減少輪轂電機(jī)的機(jī)電耦合負(fù)效應(yīng).陳辛波等[10]設(shè)計(jì)了帶有吸振器的三自由度車模型,改善了電動汽車運(yùn)行的平順性.Shao 等[11]建立了1/4 汽車主動懸架與開關(guān)磁阻電機(jī)模型,數(shù)值模擬發(fā)現(xiàn)路面不平度與電機(jī)氣隙偏心率和不平衡殘余力間存在高度耦合.李哲等[12]提出了電磁主動懸架多目標(biāo)粒子群優(yōu)化方法,改善了電動汽車垂向動力學(xué)響應(yīng).馮桂珍等[13-14]考慮了電動汽車-路面系統(tǒng)機(jī)電耦合振動特性,研究了多種因素對電動汽車平順性的影響.李韶華等[15]考慮了電機(jī)電磁激勵、地基非線性和胎路多點(diǎn)接觸,建立了電動汽車-路面系統(tǒng)機(jī)電耦合動力學(xué)模型,分析了不同因素對汽車和路面響應(yīng)的影響.以上對電動汽車機(jī)電耦合振動的研究中,未涉及電動汽車在橋上行駛的工況,未考慮車-橋耦合振動對車輛動力學(xué)特性的影響.

    此外,在智能駕駛電動汽車的研究中,其重要關(guān)注點(diǎn)之一是協(xié)同式多車隊(duì)列控制,即車輛編隊(duì)行駛[16-19].現(xiàn)有關(guān)于智能駕駛車隊(duì)的研究多針對于控制策略方面,車輛編隊(duì)行駛時的動力學(xué)問題同樣值得關(guān)注,如多輛電動汽車編隊(duì)過橋時的動力學(xué)行為.

    關(guān)于車橋耦合作用問題,近年來一直是交通運(yùn)輸領(lǐng)域研究熱點(diǎn)之一[20].在應(yīng)用梁、板理論建立橋梁模型進(jìn)行處理時,伽遼金法是一種常用分析方法,可將控制方程離散得到常微分模態(tài)方程組[21-22].Sheng 等[23]考慮了橋梁幾何非線性效應(yīng),得到了車-橋耦合非線性振動模型,數(shù)值分析了多種因素對橋梁和車輛振動特性的影響.Li 等[24]應(yīng)用達(dá)朗貝爾原理和伽遼金法,考慮曲線梁“彎-扭耦合作用”[25],對三向耦合力作用下曲線梁橋振動響應(yīng)進(jìn)行了理論分析.Ma 等[26]采用三軸11 自由度HS20-44 卡車模型建立了車-橋耦合振動系統(tǒng),分析了15 座連續(xù)梁橋的動力放大系數(shù),發(fā)現(xiàn)運(yùn)動速度對共振有重要影響.鄧露等[27]建立了13 座常見中小跨徑公路混凝土簡支梁橋的實(shí)體有限元模型,結(jié)合能表征中國設(shè)計(jì)車輛荷載動力特性的三維車輛數(shù)值模型,研究了不同因素對動力沖擊系數(shù)的影響.Zhang 等[28]在研究車-橋相互作用問題時提出了一種非線性多彈簧輪胎力模型,基于此模型分析了路面不平順和車輛跳離路面時的系統(tǒng)動力學(xué)特性.當(dāng)前研究中,主要研究傳統(tǒng)汽車的車橋耦合作用問題,關(guān)于電動汽車與橋梁的動力相互作用分析,還未見報(bào)道.

    獨(dú)立驅(qū)動電動汽車簧下質(zhì)量大,車輪振動劇烈,相對于道路來說,與橋梁的動力相互作用更加明顯.本文在此背景下,考慮橋面不平順激勵、電機(jī)激勵及車輪耦合激勵的綜合影響,研究獨(dú)立驅(qū)動電動汽車-橋梁耦合動力學(xué)行為,建模時考慮了汽車輪胎和懸架的非線性剛度特性以及輪胎與路面的胎-路多點(diǎn)接觸關(guān)系[29].本文的創(chuàng)新點(diǎn)主要在于: (1)考慮了獨(dú)立驅(qū)動電動汽車和橋梁耦合振動的相互影響;(2)嘗試分析了橋面上多輛獨(dú)立驅(qū)動電動汽車同時行駛的情況(兩車和三車).所建模型有望為智能駕駛電動汽車與橋梁的耦合動力學(xué)研究及智能駕駛車隊(duì)過橋時的工況參數(shù)設(shè)定提供一定的理論參考.

    1 基于四分之一電動汽車模型的多車-橋梁耦合系統(tǒng)建模

    考慮車身和車輪的垂向振動,建立了二自由度四分之一車輛模型,輪轂電機(jī)與車輪固連.橋梁采用兩端簡支歐拉-伯努利梁模型,車、橋的接觸采用輪胎與橋面的多點(diǎn)接觸關(guān)系.考慮了兩輛參數(shù)相同的輪轂電機(jī)驅(qū)動電動汽車以間距l(xiāng)在跨度為L的橋上行駛,電車懸架具有剛度平方、立方非線性特性,輪胎具有剛度平方非線性特性.車-橋耦合系統(tǒng)分析模型如圖1 所示.

    考慮兩車自左向右行駛,應(yīng)用達(dá)朗貝爾原理和歐拉-伯努利梁理論[30-31],可推導(dǎo)出圖1 中獨(dú)立驅(qū)動電動汽車和橋面的垂向振動控制方程為

    圖1 多車-橋梁耦合系統(tǒng)模型Fig.1 Multi-vehicle-bridge coupling system model

    其中式(1)~式(4)分別表示前車和后車的動力學(xué)方程,z1和z3為非簧載質(zhì)量位移,z2和z4為簧載質(zhì)量位移,m1,m2分別為非簧載質(zhì)量、簧載質(zhì)量,md為電機(jī)及減速機(jī)構(gòu)質(zhì)量,k1,k2分別為輪胎和懸架剛度系數(shù),c1,c2分別為輪胎和懸架阻尼系數(shù),n為輪胎與路面的多點(diǎn)接觸數(shù)目,β1,β2和 β3分別為電車懸架剛度的平方、立方非線性和輪胎剛度的平方非線性系數(shù),F1d,F2d分別為電機(jī)垂向激勵的合力.E為梁材料的彈性模量,I為梁截面對中性軸的慣性矩,m為單位長度橋梁的質(zhì)量,F1,F2分別為前、后電動汽車與橋面的相互作用力,δ 為狄拉克函數(shù),r為橋面不平順激勵,x為汽車行駛方向的位移,x1t,x2t分別為前車和后車在t時刻所在的位置.w為橋梁的垂向振動位移,w1,w2分別為車橋耦合振動引起的前車和后車下方橋面二次位移激勵.

    本文中考慮輪胎與橋面的多點(diǎn)接觸關(guān)系,即假設(shè)輪胎模型與橋面的接觸是由多個接觸點(diǎn)組成的線接觸.線接觸縱向跨度為0.25 m,接觸點(diǎn)在縱向均勻分布,選取6 個等分點(diǎn)進(jìn)行計(jì)算,即取n=6.

    考慮前后兩電動汽車以速度v勻速行駛,間距為l,則式(5)中x1t,x2t可表示為

    橋面不平順幅值采用正弦激勵的形式,表示為

    其中,Ω=2πv/L0為路面激勵頻率;B0,L0表示路面不平順幅值及波長.

    文中電機(jī)選用開關(guān)磁阻電機(jī),其具有結(jié)構(gòu)簡單、啟動性好、效率高、調(diào)速范圍寬等優(yōu)點(diǎn),適合在電動汽車上應(yīng)用.參照文獻(xiàn)[13,32-33],電機(jī)垂向激勵合力形式可表述為

    其中,T=60/(aNr) 為周期,a=1500 r/min 為電機(jī)轉(zhuǎn)速,Nr=6為定子級數(shù),R=0.05 m為 定子半徑,Te=165 N·m為電磁轉(zhuǎn)矩,b=0.001 m 為定轉(zhuǎn)子極間最短距離,i=1 A為額定電流,Kθ=82.5 為L,θ 關(guān)系曲線斜率,Lmin=4.95 H為最小電感,r=0.047 m 為轉(zhuǎn)子半徑.

    式(5)中,前、后電動汽車與橋面的相互作用力F1,F2,即動態(tài)輪胎力可表述為

    與單車相比,兩車過橋存在前后車和橋梁間的振動耦合,受車距影響較大.當(dāng)研究更多車同時行駛過橋時,其建模方法與以上兩車的方法相同,只需將系統(tǒng)方程增加相應(yīng)個汽車的自由度,將式(9)和式(10) 所示的車橋相互作用力擴(kuò)展到多個并代入式(5)即可.因此本文方法同樣適用于更多車同時行駛過橋的情況.

    2 耦合模型的求解方法

    描述橋梁振動的動力學(xué)方程式(5)為偏微分方程,在計(jì)算時可采用假設(shè)模態(tài)法進(jìn)行處理.利用假設(shè)模態(tài)法,將式中的時間和空間項(xiàng)分離,轉(zhuǎn)化為常微分方程,進(jìn)而應(yīng)用伽遼金截?cái)喾椒ㄟM(jìn)行近似計(jì)算.橋梁垂向振動位移可表示為

    其中,qi(t) 表示橋梁垂向第i階模態(tài)的廣義坐標(biāo);Xi(x)表示第i階振型函數(shù).在簡支邊界條件下,梁的第i階振型可表示為

    其中L表示橋梁的長度.

    本文計(jì)算時考慮了橋梁的前5 階模態(tài)函數(shù)近似,即橋梁垂向振動位移為

    將式(13)代入式(5),應(yīng)用伽遼金方法將方程乘以相應(yīng)振型函數(shù)并沿梁長 [0,L] 進(jìn)行積分.考慮振型函數(shù)的正交性以及狄拉克函數(shù)的性質(zhì)

    通過代數(shù)運(yùn)算,可得到電動汽車-橋梁耦合系統(tǒng)受迫振動的時變動力學(xué)方程組,表示為

    其中,X為位移向量,M為質(zhì)量矩陣,C為阻尼矩陣,K為剛度矩陣,F為荷載向量.為了便于編程計(jì)算,由懸架和輪胎剛度非線性引起的非線性力項(xiàng)寫入F列陣中.

    由于電動汽車在橋梁上不斷運(yùn)動,式(15)中系數(shù)矩陣M,C,K,F的參數(shù)在不斷變化,動力學(xué)方程組為時變系數(shù)的二階微分方程組,本文采用Newmark-β方法編制程序求解.根據(jù)數(shù)值積分計(jì)算出的橋梁垂向位移時程曲線,由不同時刻兩車所處位置,可提取出汽車行駛過程中路面對車輛的二次位移激勵.

    在求得每個振動響應(yīng)量值后,可采用均方根值作為振動性能指標(biāo),響應(yīng)量均方根值表示如下

    其中xi為振動響應(yīng)在各時刻的數(shù)值.

    3 算例分析

    電動汽車參數(shù)選取為: 非簧載質(zhì)量m1=40 kg,簧載質(zhì)量m2=337.5 kg,減速機(jī)構(gòu)和電機(jī)質(zhì)量md=30 kg,輪胎剛度系數(shù)k1=250 000 kN/m,懸架剛度系數(shù)k2=196 00 kN/m,輪胎阻尼系數(shù)c1=375 N·s/m,懸架阻尼系數(shù)c2=1450 N·s/m,懸架平方非線性剛度系數(shù) β1=0.1,懸架立方非線性剛度系數(shù) β2=0.6,輪胎平方非線性剛度系數(shù) β3=0.01,兩車間距l(xiāng)=10 m,車速v=10 m/s.橋梁系統(tǒng)參數(shù)選取為: 橋長L=30 m,單位長度質(zhì)量m=1500 kg/m,抗彎剛度EI=6.72×108N·m2,不平順幅值B0=0.01 m,不平順波長L0=10 m.

    3.1 有效性驗(yàn)證

    首先,研究了橋梁振動選取不同模態(tài)截?cái)嚯A數(shù)的收斂性,分析了橋梁模態(tài)截?cái)嚯A數(shù)對橋梁動力響應(yīng)的影響.考慮橋面不平順幅值B0=0.01 m、不平順波長L0=10 m,兩個電動汽車以車速v=10 m/s、車距l(xiāng)=10 m 通過橋梁,分別選取前1~5 階模態(tài)數(shù),得到了不同時刻全橋不同位置的撓度曲線和橋梁跨中撓度時程曲線如圖2 和圖3 所示.表1 給出了選取不同階模態(tài)數(shù)時橋梁跨中撓度均方根值(root mean square value,RMSV)及相鄰階數(shù)均方根值差別.

    圖2 t =1,3 s 時,選取不同階模態(tài)時全橋撓度曲線Fig.2 The deflection curves of the bridge with different modes at t=1,3 s

    圖3 選取不同階模態(tài)時橋梁跨中撓度時程曲線Fig.3 The time-history curves of mid-span deflection of bridge with different modes

    表1 選取不同階模態(tài)時橋梁跨中撓度均方根值及相鄰階數(shù)均方根值差別Table 1 The root mean square values of mid-span deflection of bridge with different modes and relative differences of adjacent modes

    由圖2、圖3 及表1 可以看出: 當(dāng)模態(tài)截?cái)嚯A數(shù)選取2 以上時,不同時刻全橋的撓度曲線基本重合;不同階模態(tài)情況下,橋梁跨中的位移時程曲線及均方根值相差非常小.圖3 和表1 中,橋梁前1 和2 階模態(tài)數(shù)和前3 和4 階模態(tài)數(shù)跨中位移重合以及均方根值相等是由于偶數(shù)階模態(tài)在跨中的位移為0.實(shí)際上,車輛荷載在橋上運(yùn)行時間較短,主要激發(fā)起橋梁的低階振動模態(tài),因此本文選取前5 階模態(tài)進(jìn)行計(jì)算,能夠很好地滿足精度需求.

    3.2 電機(jī)質(zhì)量、電機(jī)激勵、非線性因素、車距、車速對橋面振動特性的影響

    圖4 和圖5 分別給出了考慮不同電機(jī)質(zhì)量時,前車(first vehicle,FV) 和后車(second vehicle,SV)的橋面二次位移激勵以及橋面跨中垂向位移響應(yīng)時程曲線.

    由圖4 和圖5 可以看出: 隨著電機(jī)質(zhì)量的增加,橋面的振動加劇,橋面二次位移激勵和橋跨中位移響應(yīng)都增大.當(dāng)電機(jī)質(zhì)量取15 kg 和30 kg 時,橋面二次激勵均方根值分別增加3.97% 和7.94%(前車)、3.85%和7.71%(后車).尋求輪轂電機(jī)的輕量化,可以在一定程度上提高對橋梁、道路的友好性.

    圖4 選取不同電機(jī)質(zhì)量時橋面二次位移激勵時程曲線Fig.4 The time history curves of bridge deck secondary excitation with different motor masses

    圖5 選取不同電機(jī)質(zhì)量時橋梁跨中撓度時程曲線Fig.5 The time-history curves of mid-span deflection of bridge with different motor masses

    為了研究電機(jī)工作狀態(tài)對橋梁振動的影響,圖6和圖7 給出了不考慮電機(jī)激勵和考慮電機(jī)激勵并且電機(jī)轉(zhuǎn)速不同時,前車和后車的橋面二次位移激勵以及橋面跨中垂向位移響應(yīng)的時程曲線.

    圖6 選取不同電機(jī)激勵時橋面二次位移激勵時程曲線Fig.6 The time history curves of bridge deck secondary excitation with different motor excitations

    圖7 選取不同電機(jī)激勵時橋梁跨中撓度時程曲線Fig.7 The time-history curves of mid-span deflection of bridge with different motor excitations

    從圖6 和圖7 中可以看出,不考慮電機(jī)激勵以及考慮電機(jī)激勵且當(dāng)電機(jī)轉(zhuǎn)速為1500,1000 r/min時對橋面二次位移和跨中撓度的影響不大.當(dāng)電機(jī)轉(zhuǎn)速為500 r/min 時,由于激勵周期改變引起振動曲線的波動,表明此時車橋耦合效應(yīng)更明顯.總體來說,幾種工況下橋面二次位移激勵均方根值相差很小,電機(jī)激勵對橋梁的振動特性影響較小.

    進(jìn)一步,考慮了電動車輛輪胎和懸架非線性剛度系數(shù)對橋面振動的影響.圖8 中給出了考慮和不考慮非線性因素時,橋梁跨中撓度時程曲線.從圖8中可以看出,考慮和不考慮電動汽車非線性因素時,橋梁振動位移曲線幾乎重合.

    圖8 考慮和不考慮非線性因素時橋梁跨中撓度時程曲線Fig.8 The time-history curves of mid-span deflection of bridge with consider and disregard nonlinearity

    表2 給出了車輛非線性剛度系數(shù)取不同值時的兩車橋面二次位移激勵均方根值,從表中可以看出,選取較大非線性系數(shù)時,兩車的橋面二次位移激勵均方根值仍差別不大.由于剛度非線性引起的相互作用力與垂向相對位移的2 次和3 次方成比例,當(dāng)電動汽車和橋梁垂向位移較小時,非線性因素引起的作用力更小,因此,文中汽車的懸架和輪胎剛度非線性對橋梁振動特性以及路面二次激勵影響很小.

    表2 非線性系數(shù)取不同值時橋面二次位移激勵均方根值Table 2 The root mean square values of bridge deck secondary excitation with different nonlinear coefficients

    圖9 和圖10 中分別給出了選取不同的兩車間距對橋梁跨中位移和路面二次位移激勵的影響情況.從圖9 和圖10 中可以看出,車距對橋梁振動和二次位移激勵有較大影響,隨著車距的減小,橋梁振動加劇,兩車的橋面二次位移激勵明顯增大.相比于15 m 的車距,車距選取10 m 和5 m 時,橋面二次激勵均方根值分別增加了25.6% 和53.4%(前車)、25.9%和51.1%(后車),橋梁跨中撓度峰值分別增加了13.35%和20.40%.

    圖10 選取不同車距時橋梁跨中撓度時程曲線Fig.10 The time-history curves of mid-span deflection of bridge with different vehicle distances

    圖11 進(jìn)一步給出了三車過橋情況下不同車距對橋梁跨中位移的影響情況.從圖11 中同樣可以看出,車距對多車過橋時振動特性的影響非常明顯,相比于15 m 的車距,車距選取10 m 和5 m 時,橋梁跨中撓度峰值分別增加了17.12%和44.46%.因此,當(dāng)多個車輛同時通過橋梁時,車距是影響系統(tǒng)振動特性的一個重要因素.近年來,隨著智能駕駛技術(shù)的發(fā)展,智能車隊(duì)通過橋梁的情況越來越多,在車隊(duì)過橋前應(yīng)根據(jù)多車-橋梁耦合動力學(xué)分析結(jié)果合理設(shè)定車距,以保證車輛行駛的平順性和橋梁的耐久性、安全性.

    圖11 三車過橋不同車距時橋梁跨中撓度時程曲線Fig.11 The time-history curves of mid-span deflection of bridge with different vehicle distances for three vehicles

    圖12 和圖13 中給出了電動汽車不同行駛車速對橋梁跨中位移和路面二次位移激勵的影響情況.

    從圖12 和圖13 中可以看出,車速對橋梁振動和二次位移影響較明顯,車速的提升,出現(xiàn)位移峰值的時間縮短,位移峰值并非隨車速的增加單調(diào)增大.相比于車速為10 m/s,當(dāng)車速取為15 m/s 和20 m/s時,橋面二次激勵均方根值分別增加1.9%和4.4%(前車)、4.7%和5.7%(后車).圖14 進(jìn)一步給出了三車過橋時不同車速對橋梁跨中位移的影響情況.從圖14 同樣可以看出,隨著車速提升,橋梁跨中撓度出現(xiàn)峰值的時間縮短,且峰值并非隨車速增加而單調(diào)增大,這與兩車情況類似.因此,當(dāng)多個車輛經(jīng)過橋梁時,車速對車橋耦合系統(tǒng)振動的影響也是不容忽視的,在智能駕駛車隊(duì)過橋前,應(yīng)根據(jù)多車-橋梁耦合動力學(xué)分析結(jié)果合理設(shè)定車速.

    圖12 選取不同車速時橋面二次位移激勵時程曲線Fig.12 The time history curves of bridge deck secondary excitation with different vehicle speeds

    圖13 選取不同車速時橋梁跨中撓度時程曲線Fig.13 The time-history curves of mid-span deflection of bridge with different vehicle speeds

    圖14 三車過橋不同車速時橋梁跨中撓度時程曲線Fig.14 The time-history curves of mid-span deflection of bridge with different vehicle speeds for three vehicles

    3.3 電機(jī)激勵、車速、橋面不平順激勵、車距、三重耦合激勵對車輛行駛平順性的影響

    針對于文中電動汽車-橋梁耦合系統(tǒng),進(jìn)一步分析了不同工況對電動汽車行駛平順性的影響情況.車輛的行駛平順性指標(biāo)主要為: 車身加速度、懸架動撓度、輪胎相對動載.本文研究的獨(dú)立驅(qū)動電動汽車,由于簧下質(zhì)量振動劇烈,補(bǔ)充了車輪加速度為評價指標(biāo)之一.

    類似地,第二輛電動汽車的車身垂向加速度指標(biāo)表示為車輪垂向加速度指標(biāo)表示為懸架動撓度指標(biāo)表示為z4-z3,輪胎相對動載表示為

    在后文研究中,定義指標(biāo)1 為車身垂向加速度,指標(biāo)2 為懸架動撓度,指標(biāo)3 為車輪垂向加速度,指標(biāo)4 為輪胎相對動載.圖15 給出了前車4 個平順性指標(biāo)的位移時程曲線.可以看出由于電機(jī)激勵的存在,導(dǎo)致車輪垂向加速度值增大.

    圖15 前車4 個平順性指標(biāo)時程曲線Fig.15 The time history curves of four ride comfort indexes of first vehicle

    類似3.2 節(jié),懸架和輪胎剛度非線性對電動汽車振動特性影響很小,因此在算例分析中并未給出結(jié)果.

    首先研究了電機(jī)激勵對電動汽車振動特性的影響情況.表3 給出了不同電機(jī)激勵工況下兩車的4 個平順性指標(biāo)均方根值,其中電機(jī)激勵(motor excitation,ME) 和無激勵(no motor excitation,NME)考慮電機(jī)質(zhì)量,無電機(jī)(no motor,NM)不考慮電機(jī)質(zhì)量和激勵.

    從表3 中可以看出: 電機(jī)激勵使車身垂向加速度增大明顯,這降低了車輛的乘坐舒適性;懸架動撓度值略有增加,這對懸架行程限度的要求變大;輪胎相對動載增大,略微增大了輪胎、地面附著能力并降低了道路友好性;電機(jī)激勵對車輪垂向加速度影響非常明顯,有無電機(jī)激勵,車輪垂向加速度差距顯著.綜合來看,電機(jī)激勵在一定程度上降低了車輛平順性和道路友好性.

    表3 不同電機(jī)激勵情況下平順性指標(biāo)均方根值Table 3 The root mean square values of ride comfort indexes with different motor excitations

    進(jìn)一步分析了車速對電動汽車振動特性的影響,表4 給出了不同車速工況下兩車的4 個平順性指標(biāo)均方根值.表4 中的數(shù)據(jù)表明: 車速增加使車身垂向加速度增大較明顯,降低了車輛的乘坐舒適性;車速對車輪垂向加速度影響不大,車輪垂向加速度主要受電機(jī)激勵影響;輪胎相對動載增大,增大了輪胎、地面附著能力并降低了道路友好性.可以看出,車速對車輛平順性和道路友好性也有較大影響.

    表4 不同車速情況下平順性指標(biāo)均方根值Table 4 The root mean square values of ride comfort indexes with different vehicle speeds

    橋面不平順激勵是影響汽車振動特性的主要激勵形式之一.表5 給出了考慮不同橋面不平整度時的電動汽車平順性指標(biāo)均方根值.由表5 中可以看出: 和車速工況類似,橋面不平順度對車輪垂向加速度影響不大,其主要受電機(jī)激勵影響;當(dāng)不平順幅值B0越大、不平順波長L0越小,即橋面越不平坦時,車身垂向加速度、懸架動撓度和輪胎相對動載的值越大;且當(dāng)B0越大時,隨著L0的改變,車身垂向加速度、懸架動撓度和輪胎相對動載的值變化越大.表5中數(shù)據(jù)表明橋面不平順激勵是影響車輛平順性和道路友好性的重要因素.

    表5 不同橋面不平順情況下平順性指標(biāo)均方根值Table 5 The root mean square values of ride comfort indexes with different bridge irregularities

    表6 給出了考慮不同車距時的電動汽車平順性指標(biāo)均方根值.由表6 中可以看出: 車距的改變對第一輛車4 個指標(biāo)的影響并不大;相較而言,車距的改變對第二輛車產(chǎn)生了較大的影響,除車輪垂向加速度外其余指標(biāo)均發(fā)生了較明顯變化.

    表6 不同車距情況下平順性指標(biāo)均方根值Table 6 The root mean square values of ride comfort indexes with different vehicle distances

    進(jìn)一步給出了不同車距下車身垂向加速度、懸架動撓度和輪胎相對動載的變化情況,如圖16 所示.從圖16 中可以清晰看出,車距的改變對前車的車身垂向加速度、懸架動撓度和輪胎相對動載幾乎無影響,而對后車的車身垂向加速度、懸架動撓度和輪胎相對動載產(chǎn)生明顯影響,其中對輪胎相對動載影響最顯著.在車距為7 m 和17 m 時,后車相對動載達(dá)到峰值,而在車距為10 m 時,后車相對動載與前車幾乎無變化.后車車身垂向加速度和懸架動撓度變化幅值不及相對動載,但趨勢類似,均在車距為5 m 和17 m 達(dá)到最小值,在車距為11 m 時達(dá)到峰值(略大于前車數(shù)值).車身垂向加速度、懸架動撓度和輪胎相對動載的變化會影響車輛乘坐舒適性、懸架行程限度、胎路附著能力以及道路友好性,因此在研究多車輛過橋問題時,車輛行駛間距的影響不容忽視.

    圖16 不同車距情況下兩車指標(biāo)均方根值Fig.16 The root mean square values of the ride comfort indexes with different vehicle distances

    進(jìn)一步,綜合考慮了電動汽車(electric vehicle,EV)平順性指標(biāo)(ride comfort index,RCI)在三重激勵耦合作用下相比于單一激勵的差別.三重激勵為橋面不平順激勵(bridge irregularity,BI)、電機(jī)激勵以及橋面二次位移激勵(secondary excitation,SE).由于車輪垂向加速度顯著受電機(jī)激勵控制,表7 中未列出其指標(biāo)值的變化.

    由表7 可以看出: 綜合考慮3 種激勵(橋面不平順激勵,電機(jī)激勵,橋面二次位移激勵)與考慮2 種激勵(橋面不平順激勵,橋面二次位移激勵)及只考慮單一激勵(橋面不平順激勵)的對比,考慮三重激勵時指標(biāo)幅值明顯增加;當(dāng)橋面越平坦時(即不平順幅值B0越小),平順性指標(biāo)增加越明顯,對電動汽車的影響越顯著.

    表7 三重激勵相較于橋面不平順單一激勵的平順性指標(biāo)均方根值對比Table 7 The comparison of root mean square values of ride comfort indexes between different excitations

    4 結(jié)論

    本文建立了獨(dú)立驅(qū)動電動汽車-橋梁耦合系統(tǒng)動力學(xué)模型,應(yīng)用假設(shè)模態(tài)法和伽遼金方法將描述橋梁振動的偏微分方程轉(zhuǎn)化為常微分方程并與車輛振動方程聯(lián)立,利用Newmark-β 法進(jìn)行數(shù)值積分求解.分析計(jì)算了多車-橋梁耦合系統(tǒng)的動力學(xué)響應(yīng),研究了不同因素對橋梁和電動汽車動力學(xué)響應(yīng)的影響情況,得出的主要結(jié)論如下.

    (1)電機(jī)質(zhì)量主要對橋面振動有影響,電機(jī)激勵對電動汽車加速度指標(biāo)影響較大.

    (2)在電動汽車垂向相對位移較小時,汽車剛度非線性因素對橋梁和汽車振動特性影響很小,為簡化計(jì)算可考慮線性模型建模.

    (3)當(dāng)考慮多車過橋時,車距和車速對橋梁振動特性和電動汽車平順性指標(biāo)影響都很大.在智能駕駛車隊(duì)過橋前,應(yīng)對多車-橋梁耦合系統(tǒng)進(jìn)行動態(tài)設(shè)計(jì),以合理設(shè)定車距和車速,從而保證車輛的平順性和橋梁的耐久性、安全性.

    (4)橋面不平順激勵是影響車輛振動特性的重要因素.當(dāng)綜合考慮耦合激勵影響時,橋面越平坦,電機(jī)激勵及橋面二次激勵對車輛平順性和道路友好性的影響越加顯著.當(dāng)汽車行駛在平坦橋面時,電機(jī)激勵及橋面二次激勵對獨(dú)立驅(qū)動電動汽車的影響不容忽視.

    猜你喜歡
    橋梁振動
    振動的思考
    噴水推進(jìn)高速艇尾部振動響應(yīng)分析
    This “Singing Highway”plays music
    手拉手 共搭愛的橋梁
    句子也需要橋梁
    振動攪拌 震動創(chuàng)新
    中國公路(2017年18期)2018-01-23 03:00:38
    中立型Emden-Fowler微分方程的振動性
    加固技術(shù)創(chuàng)新,為橋梁健康保駕護(hù)航
    中國公路(2017年11期)2017-07-31 17:56:30
    無人機(jī)在橋梁檢測中的應(yīng)用
    中國公路(2017年10期)2017-07-21 14:02:37
    高性能砼在橋梁中的應(yīng)用
    69精品国产乱码久久久| 侵犯人妻中文字幕一二三四区| 丁香六月欧美| 免费在线观看日本一区| xxx96com| 两个人看的免费小视频| 黄色丝袜av网址大全| √禁漫天堂资源中文www| 欧美大码av| 黄色a级毛片大全视频| 欧洲精品卡2卡3卡4卡5卡区| av天堂久久9| 欧美日韩亚洲国产一区二区在线观看| 夫妻午夜视频| 午夜免费鲁丝| 色播在线永久视频| 日日摸夜夜添夜夜添小说| 99国产精品99久久久久| 国产精品电影一区二区三区| 免费在线观看日本一区| 一区二区三区国产精品乱码| 窝窝影院91人妻| 怎么达到女性高潮| 淫妇啪啪啪对白视频| 精品一区二区三区av网在线观看| 精品久久久久久成人av| 欧美一级毛片孕妇| 少妇的丰满在线观看| 18禁黄网站禁片午夜丰满| 国产欧美日韩综合在线一区二区| 免费不卡黄色视频| 黄色片一级片一级黄色片| 99在线视频只有这里精品首页| 18美女黄网站色大片免费观看| 这个男人来自地球电影免费观看| 操出白浆在线播放| www.精华液| 高潮久久久久久久久久久不卡| 国产99久久九九免费精品| 很黄的视频免费| 国产精品成人在线| 成人黄色视频免费在线看| 精品无人区乱码1区二区| 亚洲成a人片在线一区二区| 国产亚洲欧美精品永久| 日韩一卡2卡3卡4卡2021年| 亚洲成a人片在线一区二区| 亚洲人成伊人成综合网2020| 精品无人区乱码1区二区| 两性夫妻黄色片| 在线观看日韩欧美| aaaaa片日本免费| 91麻豆av在线| 国产成人啪精品午夜网站| 欧美日韩亚洲国产一区二区在线观看| 正在播放国产对白刺激| 午夜免费激情av| 在线观看www视频免费| 午夜91福利影院| 在线观看舔阴道视频| 好看av亚洲va欧美ⅴa在| 久久青草综合色| 黄片播放在线免费| 久久国产亚洲av麻豆专区| 可以在线观看毛片的网站| 亚洲欧美日韩高清在线视频| 欧美最黄视频在线播放免费 | 亚洲欧美精品综合久久99| 最新在线观看一区二区三区| 欧美成人免费av一区二区三区| 国产高清videossex| 女人精品久久久久毛片| 亚洲成国产人片在线观看| 久久伊人香网站| 一级毛片女人18水好多| 成人精品一区二区免费| 欧美乱码精品一区二区三区| 高清欧美精品videossex| 9热在线视频观看99| 岛国在线观看网站| 国产精品98久久久久久宅男小说| 国产精品98久久久久久宅男小说| 亚洲精品中文字幕在线视频| 国产三级在线视频| 99国产极品粉嫩在线观看| 巨乳人妻的诱惑在线观看| 国产蜜桃级精品一区二区三区| 露出奶头的视频| 中文字幕人妻丝袜一区二区| 久久香蕉国产精品| 黑丝袜美女国产一区| 一个人观看的视频www高清免费观看 | 久久久久亚洲av毛片大全| 老司机靠b影院| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成网站在线播放欧美日韩| 侵犯人妻中文字幕一二三四区| 999久久久国产精品视频| ponron亚洲| 亚洲自拍偷在线| 另类亚洲欧美激情| 中文字幕另类日韩欧美亚洲嫩草| 一边摸一边抽搐一进一小说| 中文字幕人妻熟女乱码| 自拍欧美九色日韩亚洲蝌蚪91| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文日韩欧美视频| 一本大道久久a久久精品| 亚洲欧美激情在线| 午夜精品国产一区二区电影| 亚洲伊人色综图| 亚洲成人精品中文字幕电影 | 久久人妻熟女aⅴ| 国产一区二区三区在线臀色熟女 | 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人av教育| 精品午夜福利视频在线观看一区| 国产精品 欧美亚洲| 国产精品野战在线观看 | 久久久精品欧美日韩精品| 免费人成视频x8x8入口观看| 如日韩欧美国产精品一区二区三区| 欧美丝袜亚洲另类 | 在线播放国产精品三级| 50天的宝宝边吃奶边哭怎么回事| netflix在线观看网站| 99国产精品一区二区三区| 亚洲国产精品合色在线| 变态另类成人亚洲欧美熟女 | 久久国产乱子伦精品免费另类| 男女做爰动态图高潮gif福利片 | 日韩免费av在线播放| 在线视频色国产色| av视频免费观看在线观看| 91成年电影在线观看| 亚洲国产精品一区二区三区在线| 国产精品久久视频播放| 久99久视频精品免费| 成年女人毛片免费观看观看9| xxx96com| 亚洲国产精品合色在线| ponron亚洲| 日本欧美视频一区| 在线天堂中文资源库| 黄色视频不卡| 99riav亚洲国产免费| 午夜福利在线免费观看网站| 一区二区三区精品91| 日本vs欧美在线观看视频| 国产激情久久老熟女| 9191精品国产免费久久| 国产精品一区二区精品视频观看| 妹子高潮喷水视频| 最好的美女福利视频网| 在线免费观看的www视频| 男女下面插进去视频免费观看| 99热国产这里只有精品6| 精品国产乱码久久久久久男人| 欧美日韩av久久| av免费在线观看网站| 国产深夜福利视频在线观看| 在线观看免费视频网站a站| 99国产精品免费福利视频| 69av精品久久久久久| 女警被强在线播放| 伦理电影免费视频| 亚洲一区二区三区不卡视频| 动漫黄色视频在线观看| 日本欧美视频一区| 亚洲狠狠婷婷综合久久图片| 亚洲精品中文字幕一二三四区| 精品久久久久久久毛片微露脸| 欧美成人性av电影在线观看| 嫁个100分男人电影在线观看| 久久人妻熟女aⅴ| 一边摸一边抽搐一进一小说| 久久久久久久精品吃奶| 色婷婷av一区二区三区视频| 黑人猛操日本美女一级片| 亚洲国产精品合色在线| 午夜激情av网站| 伊人久久大香线蕉亚洲五| 亚洲精品久久成人aⅴ小说| 女生性感内裤真人,穿戴方法视频| 日本a在线网址| aaaaa片日本免费| 黄色a级毛片大全视频| 亚洲色图av天堂| 欧美一区二区精品小视频在线| 色在线成人网| 亚洲熟妇熟女久久| 91国产中文字幕| 丰满人妻熟妇乱又伦精品不卡| 国产伦一二天堂av在线观看| 男人舔女人的私密视频| 国产一卡二卡三卡精品| 国产精品久久视频播放| √禁漫天堂资源中文www| 在线十欧美十亚洲十日本专区| 午夜精品国产一区二区电影| 88av欧美| 精品福利永久在线观看| 真人一进一出gif抽搐免费| 天堂影院成人在线观看| 日韩精品中文字幕看吧| 亚洲自拍偷在线| videosex国产| 欧美色视频一区免费| 99久久综合精品五月天人人| 丝袜美足系列| 精品国产超薄肉色丝袜足j| 精品一品国产午夜福利视频| 天堂中文最新版在线下载| 波多野结衣一区麻豆| 精品第一国产精品| 国产成+人综合+亚洲专区| 国产精品综合久久久久久久免费 | 无限看片的www在线观看| 一进一出抽搐gif免费好疼 | 欧美另类亚洲清纯唯美| 在线观看舔阴道视频| 超碰97精品在线观看| 国产蜜桃级精品一区二区三区| 天堂中文最新版在线下载| 精品无人区乱码1区二区| 成人精品一区二区免费| 久久久久久人人人人人| 精品乱码久久久久久99久播| 免费人成视频x8x8入口观看| 精品午夜福利视频在线观看一区| 黄色怎么调成土黄色| 亚洲精品成人av观看孕妇| 精品一品国产午夜福利视频| 男女之事视频高清在线观看| 日韩一卡2卡3卡4卡2021年| 欧美日韩亚洲国产一区二区在线观看| 亚洲专区字幕在线| 自线自在国产av| 成人免费观看视频高清| 中文字幕人妻熟女乱码| 高潮久久久久久久久久久不卡| 久久精品亚洲熟妇少妇任你| 91成人精品电影| 日韩欧美免费精品| aaaaa片日本免费| 香蕉久久夜色| 亚洲精品一区av在线观看| 露出奶头的视频| 女人被狂操c到高潮| 好看av亚洲va欧美ⅴa在| 大型黄色视频在线免费观看| 91成年电影在线观看| x7x7x7水蜜桃| 亚洲精品av麻豆狂野| 黄色 视频免费看| 日韩欧美三级三区| 国产成人影院久久av| 成人18禁在线播放| 99国产精品一区二区三区| 天天影视国产精品| 亚洲av日韩精品久久久久久密| 他把我摸到了高潮在线观看| 在线观看午夜福利视频| 天堂影院成人在线观看| 成年版毛片免费区| 国产蜜桃级精品一区二区三区| 免费搜索国产男女视频| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 亚洲精品在线观看二区| 国产有黄有色有爽视频| 国产精品亚洲av一区麻豆| 真人做人爱边吃奶动态| 亚洲成人免费电影在线观看| 美女高潮喷水抽搐中文字幕| 一边摸一边抽搐一进一小说| 黑人欧美特级aaaaaa片| 亚洲va日本ⅴa欧美va伊人久久| 男人舔女人的私密视频| 亚洲av电影在线进入| 久久久国产成人免费| 日日爽夜夜爽网站| 国产av在哪里看| 俄罗斯特黄特色一大片| 久久久久久亚洲精品国产蜜桃av| 精品国产国语对白av| 村上凉子中文字幕在线| 欧美成人性av电影在线观看| 欧美日韩亚洲国产一区二区在线观看| 欧美乱色亚洲激情| 日本欧美视频一区| 老鸭窝网址在线观看| 女人被狂操c到高潮| 欧美日韩av久久| 日韩欧美一区二区三区在线观看| 在线观看午夜福利视频| 50天的宝宝边吃奶边哭怎么回事| 老熟妇仑乱视频hdxx| 成人免费观看视频高清| 欧美中文日本在线观看视频| 757午夜福利合集在线观看| 两个人免费观看高清视频| 中文欧美无线码| 精品一区二区三区四区五区乱码| 国产97色在线日韩免费| 欧美乱码精品一区二区三区| 亚洲伊人色综图| 99精国产麻豆久久婷婷| av福利片在线| 咕卡用的链子| 亚洲欧美日韩无卡精品| 男女下面插进去视频免费观看| 人妻久久中文字幕网| 水蜜桃什么品种好| 国产91精品成人一区二区三区| 国产色视频综合| 9191精品国产免费久久| 亚洲成人免费av在线播放| 交换朋友夫妻互换小说| 久久久久精品国产欧美久久久| 欧美成人午夜精品| 美女高潮喷水抽搐中文字幕| 人人妻人人澡人人看| 91在线观看av| 日本欧美视频一区| 啪啪无遮挡十八禁网站| 伦理电影免费视频| 国产成人欧美| 国产三级黄色录像| 国产97色在线日韩免费| 色老头精品视频在线观看| 亚洲国产精品999在线| 国产亚洲精品久久久久5区| 亚洲男人的天堂狠狠| 国产视频一区二区在线看| 国产一区二区在线av高清观看| 欧美中文日本在线观看视频| 亚洲一区中文字幕在线| 80岁老熟妇乱子伦牲交| 久9热在线精品视频| 热re99久久国产66热| 男女午夜视频在线观看| 欧美激情高清一区二区三区| 91麻豆av在线| 热re99久久精品国产66热6| 久久亚洲精品不卡| 日本欧美视频一区| 神马国产精品三级电影在线观看 | 日本a在线网址| 欧美日本中文国产一区发布| 男女下面进入的视频免费午夜 | cao死你这个sao货| 91国产中文字幕| 99国产精品99久久久久| 99re在线观看精品视频| 麻豆av在线久日| 亚洲精品av麻豆狂野| 色在线成人网| 淫秽高清视频在线观看| 成人三级做爰电影| 国产乱人伦免费视频| 天堂动漫精品| av片东京热男人的天堂| 亚洲专区字幕在线| 一区二区日韩欧美中文字幕| 国产精品av久久久久免费| 一级作爱视频免费观看| 欧美最黄视频在线播放免费 | 久久中文看片网| 精品一区二区三区四区五区乱码| 激情在线观看视频在线高清| 99国产精品免费福利视频| 99久久精品国产亚洲精品| 中文字幕精品免费在线观看视频| 精品一品国产午夜福利视频| 国产av又大| 久久精品国产综合久久久| 亚洲成人精品中文字幕电影 | 久久伊人香网站| 身体一侧抽搐| 丝袜美腿诱惑在线| 法律面前人人平等表现在哪些方面| 欧美日韩精品网址| 亚洲精品美女久久av网站| 两人在一起打扑克的视频| 亚洲欧美精品综合一区二区三区| 免费av毛片视频| 天天躁夜夜躁狠狠躁躁| 国产精品av久久久久免费| 伦理电影免费视频| 日韩欧美国产一区二区入口| 一级片免费观看大全| 中文字幕av电影在线播放| 日韩精品免费视频一区二区三区| 制服诱惑二区| 欧美精品啪啪一区二区三区| 美女高潮喷水抽搐中文字幕| 女人精品久久久久毛片| 欧美日韩av久久| 琪琪午夜伦伦电影理论片6080| 欧美成人免费av一区二区三区| 不卡一级毛片| 久久久精品欧美日韩精品| 亚洲国产中文字幕在线视频| 纯流量卡能插随身wifi吗| 在线看a的网站| 久久青草综合色| 免费看a级黄色片| 久久久久九九精品影院| 精品卡一卡二卡四卡免费| 国产精品久久久人人做人人爽| 高清毛片免费观看视频网站 | 亚洲少妇的诱惑av| 亚洲黑人精品在线| 国产精品野战在线观看 | 操出白浆在线播放| www日本在线高清视频| 99在线人妻在线中文字幕| 日本三级黄在线观看| 身体一侧抽搐| 母亲3免费完整高清在线观看| 91精品国产国语对白视频| 他把我摸到了高潮在线观看| 夜夜夜夜夜久久久久| 日韩精品中文字幕看吧| 国产亚洲欧美在线一区二区| 在线av久久热| 精品福利观看| 日韩视频一区二区在线观看| 日本三级黄在线观看| 最好的美女福利视频网| 中国美女看黄片| 亚洲欧美一区二区三区久久| 国产真人三级小视频在线观看| 夜夜爽天天搞| 久久影院123| 黑人巨大精品欧美一区二区蜜桃| 精品国产乱码久久久久久男人| 一级毛片高清免费大全| 欧美亚洲日本最大视频资源| 国产av一区二区精品久久| 免费搜索国产男女视频| 亚洲国产中文字幕在线视频| 手机成人av网站| a级毛片黄视频| 真人做人爱边吃奶动态| 另类亚洲欧美激情| 操美女的视频在线观看| 亚洲欧美激情在线| 女性生殖器流出的白浆| 这个男人来自地球电影免费观看| 国产熟女xx| 岛国视频午夜一区免费看| 亚洲成国产人片在线观看| 不卡一级毛片| 免费搜索国产男女视频| 麻豆一二三区av精品| 久久天堂一区二区三区四区| 女人被狂操c到高潮| 女人被躁到高潮嗷嗷叫费观| 日本免费a在线| 免费av中文字幕在线| 色在线成人网| 午夜精品久久久久久毛片777| 九色亚洲精品在线播放| 欧美日韩乱码在线| 看免费av毛片| 日本 av在线| aaaaa片日本免费| 亚洲精品中文字幕在线视频| 电影成人av| 国产熟女xx| 啦啦啦 在线观看视频| 国产高清国产精品国产三级| av中文乱码字幕在线| 在线观看午夜福利视频| 欧美久久黑人一区二区| 国产真人三级小视频在线观看| 久久久精品国产亚洲av高清涩受| 日本黄色视频三级网站网址| 亚洲七黄色美女视频| 看片在线看免费视频| 在线视频色国产色| 制服诱惑二区| 在线观看免费视频日本深夜| 久久人妻福利社区极品人妻图片| av国产精品久久久久影院| 日本 av在线| 中出人妻视频一区二区| 嫩草影视91久久| 一区在线观看完整版| 777久久人妻少妇嫩草av网站| 天天添夜夜摸| 怎么达到女性高潮| 宅男免费午夜| 成人18禁在线播放| 一级a爱片免费观看的视频| 免费在线观看亚洲国产| 中文亚洲av片在线观看爽| 久久精品国产亚洲av高清一级| 国产激情欧美一区二区| 美女高潮喷水抽搐中文字幕| 亚洲国产中文字幕在线视频| 波多野结衣av一区二区av| 视频区图区小说| 成人av一区二区三区在线看| 大香蕉久久成人网| av电影中文网址| 天天添夜夜摸| 成人精品一区二区免费| 午夜福利,免费看| 国产精品永久免费网站| 日本 av在线| 99热国产这里只有精品6| 久久精品91蜜桃| 桃红色精品国产亚洲av| 757午夜福利合集在线观看| 国产有黄有色有爽视频| 免费在线观看完整版高清| av福利片在线| 亚洲精品在线美女| 免费看十八禁软件| 一二三四在线观看免费中文在| 亚洲黑人精品在线| 色老头精品视频在线观看| 久久人妻福利社区极品人妻图片| 一级毛片精品| 国产一区二区三区综合在线观看| 国内久久婷婷六月综合欲色啪| 精品国产一区二区三区四区第35| 精品久久蜜臀av无| 在线观看一区二区三区激情| 97超级碰碰碰精品色视频在线观看| 一进一出抽搐gif免费好疼 | 久久久国产一区二区| 成人特级黄色片久久久久久久| 大陆偷拍与自拍| 精品久久蜜臀av无| 中文亚洲av片在线观看爽| 国产精品自产拍在线观看55亚洲| 美女扒开内裤让男人捅视频| 精品一区二区三区视频在线观看免费 | 亚洲精品国产色婷婷电影| 日韩免费av在线播放| 久久精品国产亚洲av高清一级| 国内久久婷婷六月综合欲色啪| 亚洲 国产 在线| 久久99一区二区三区| 制服人妻中文乱码| 欧美日韩黄片免| 99国产精品一区二区蜜桃av| 女人被狂操c到高潮| 一级a爱视频在线免费观看| 琪琪午夜伦伦电影理论片6080| 亚洲欧美精品综合久久99| 天天躁夜夜躁狠狠躁躁| 免费高清在线观看日韩| aaaaa片日本免费| 亚洲精品国产一区二区精华液| 最好的美女福利视频网| 久久久国产成人精品二区 | 欧美黑人欧美精品刺激| 久久久久久久久免费视频了| 无遮挡黄片免费观看| 欧美另类亚洲清纯唯美| 深夜精品福利| 在线观看免费高清a一片| 国产精品爽爽va在线观看网站 | 午夜免费观看网址| 最新美女视频免费是黄的| 深夜精品福利| 亚洲色图综合在线观看| 国产欧美日韩精品亚洲av| 看免费av毛片| 亚洲av日韩精品久久久久久密| 亚洲专区字幕在线| 欧美日韩乱码在线| 日韩精品中文字幕看吧| 亚洲黑人精品在线| 国产精品99久久99久久久不卡| 亚洲一码二码三码区别大吗| 高清欧美精品videossex| 亚洲七黄色美女视频| 国产精品一区二区精品视频观看| 如日韩欧美国产精品一区二区三区| 亚洲伊人色综图| 午夜久久久在线观看| 天堂动漫精品| 好看av亚洲va欧美ⅴa在| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品国产精品久久久不卡| 99国产精品99久久久久| 满18在线观看网站| 欧美成人午夜精品| 国产成人一区二区三区免费视频网站| 国产又爽黄色视频| 国内久久婷婷六月综合欲色啪| 啦啦啦 在线观看视频| 超碰成人久久| 久久精品亚洲av国产电影网| 大型黄色视频在线免费观看| 亚洲成国产人片在线观看| 欧美激情 高清一区二区三区| 久久 成人 亚洲| 精品福利永久在线观看| 亚洲成人久久性| 两人在一起打扑克的视频| 美女午夜性视频免费| 亚洲精品国产色婷婷电影| 99热只有精品国产| 久久久国产成人精品二区 | 国产一区在线观看成人免费| 国内久久婷婷六月综合欲色啪| 一区二区日韩欧美中文字幕| 精品第一国产精品|