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

    考慮梁彎扭耦合對(duì)鋼軌橫向振動(dòng)特性的影響

    2023-11-03 03:43:20豆銀玲王紹華
    關(guān)鍵詞:膠墊共振頻率扣件

    豆銀玲 ,韋 凱 ,曹 勇 ,王紹華 ,亓 偉 ,王 平

    (1.西南交通大學(xué)高速鐵路線路工程教育部重點(diǎn)實(shí)驗(yàn)室,四川 成都 610031;2.成都工業(yè)職業(yè)技術(shù)學(xué)院,四川 成都610031;3.西南交通大學(xué)土木工程學(xué)院,四川 成都 610031;4.奧克蘭大學(xué)機(jī)械學(xué)院,奧克蘭 1010)

    隨著我國城市軌道交通規(guī)模的擴(kuò)大以及運(yùn)營速度的不斷增加,軌道交通所引發(fā)的環(huán)境振動(dòng)污染問題也日益顯著.因此,不少學(xué)者通常采用數(shù)值仿真的方法對(duì)軌道系統(tǒng)的振動(dòng)和噪聲進(jìn)行預(yù)測[1-3].而建立合理的鋼軌-扣件耦合模型和選取符合實(shí)際運(yùn)營條件的軌道結(jié)構(gòu)動(dòng)參數(shù)對(duì)準(zhǔn)確地分析鋼軌振動(dòng)問題極為重要.

    在以往研究中,長直無限長鋼軌主要采用經(jīng)典Timoshenko 梁模擬.由于鋼軌橫截面的剪切中心與質(zhì)心不重合,在橫向荷載作用下,鋼軌橫截面質(zhì)心偏移將引起鋼軌橫向彎曲與扭轉(zhuǎn)的耦合變形,即彎扭耦合變形[4].而經(jīng)典的Timoshenko 梁不能考慮鋼軌彎曲和扭轉(zhuǎn)變形的耦合作用.此外,車輛在軌道上運(yùn)行時(shí),鋼軌會(huì)承受來自橫向、垂向與扭轉(zhuǎn)荷載的共同作用.Vincent 等[5]指出垂向和橫向?qū)壍老到y(tǒng)輻射噪聲的相對(duì)貢獻(xiàn)量主要取決于車輪與鋼軌的接觸位置以及各個(gè)波沿鋼軌的衰減特征.如果接觸點(diǎn)明顯偏移,由橫向振動(dòng)引起的輻射噪聲可能會(huì)達(dá)到甚至超過垂向的貢獻(xiàn).Kostovasilis 等[6]提出忽略扭轉(zhuǎn)與橫向的耦合作用將會(huì)低估鋼軌的橫向振動(dòng)響應(yīng).通過引入經(jīng)驗(yàn)系數(shù)能近似考慮垂橫向之間的交叉耦合作用,但該系數(shù)不具有普遍性.Timoshenko 等[7]首次提出了薄壁梁的彎扭耦合振動(dòng)理論.隨后該理論得到廣泛應(yīng)用,如航空航天、汽車和民用建筑行業(yè)等.Li 等[8]利用解析法分析軸向加載的復(fù)合薄壁Timoshenko 梁在各種集中和分布荷載作用下的彎曲與扭轉(zhuǎn)耦合動(dòng)力響應(yīng).易強(qiáng)[9]采用考慮彎扭耦合的Timoshenko 梁模型描述高速鐵路無砟軌道結(jié)構(gòu),通過對(duì)比現(xiàn)場錘擊得到的鋼軌橫向一階pinnedpinned 共振頻率與理論計(jì)算的結(jié)果,發(fā)現(xiàn)兩者較為吻合.農(nóng)興中等[10]通過測試地鐵常用減振軌道鋼軌的橫向加速度導(dǎo)納,發(fā)現(xiàn)鋼軌的橫向振動(dòng)特性比垂向振動(dòng)更為復(fù)雜,且鋼軌的橫向振動(dòng)包括橫向彎曲波和扭轉(zhuǎn)波.綜上可知,如果不考慮鋼軌橫向彎曲與扭轉(zhuǎn)的耦合作用,將可能影響實(shí)際橫向彈性波在鋼軌中的傳播特性.

    另一方面,扣件膠墊動(dòng)參數(shù)也是影響軌道動(dòng)力特性的重要結(jié)構(gòu)參數(shù).為了更準(zhǔn)確地預(yù)測鋼軌的橫向振動(dòng)響應(yīng),需考慮實(shí)際運(yùn)行條件下列車附加的垂向預(yù)壓對(duì)扣件膠墊等高分子材料的動(dòng)力性能的影響[11].Thompson 等[12]通過測試膠墊扭轉(zhuǎn)剛度發(fā)現(xiàn),由于扭轉(zhuǎn)剛度測試頻率范圍有限,可根據(jù)均質(zhì)材料特性用膠墊垂向剛度估計(jì)扭轉(zhuǎn)剛度.因此,垂向預(yù)壓力變化將引起扭轉(zhuǎn)剛度產(chǎn)生相應(yīng)的黏彈性動(dòng)力性能,而采用本文提出的考慮彎扭耦合的Timoshenko可分析扣件膠墊垂向預(yù)壓特性對(duì)軌道橫向動(dòng)力響應(yīng)特性的影響.

    鑒于此,為準(zhǔn)確預(yù)測彈性波在鋼軌中的傳播特性,亟需建立合理的鋼軌-扣件耦合模型和考慮扣件膠墊垂向預(yù)壓特性的影響.為此,本文基于波譜-辛混合法建立了考慮梁彎扭耦合的整體道床無砟軌道空間無限長模型,通過對(duì)比分析未考慮鋼軌橫向質(zhì)心偏移與采用經(jīng)典Timoshenko 梁模型的計(jì)算結(jié)果,驗(yàn)證本文建立模型的正確性.在此基礎(chǔ)上,探究了不同預(yù)壓下扣件膠墊動(dòng)態(tài)黏彈性力學(xué)性能對(duì)鋼軌橫向速度導(dǎo)納的影響規(guī)律.最后,通過對(duì)比分析實(shí)測與理論計(jì)算的不同預(yù)壓下鋼軌橫向一階彎曲共振頻率,從而闡明考慮梁彎扭耦合和扣件膠墊動(dòng)態(tài)黏彈性力學(xué)性能的必要性.

    1 不同垂向預(yù)壓下鋼軌橫向錘擊試驗(yàn)

    為了獲取不同垂向預(yù)壓力下鋼軌橫向一階彎曲共振頻率,以地鐵常用的DZIII 型鋼軌-扣件系統(tǒng)為測試對(duì)象,通過千斤頂在垂向施加不同的預(yù)壓力,利用力錘橫向錘擊鋼軌軌頭得到不同預(yù)壓力下鋼軌橫向加速度導(dǎo)納.

    1.1 試驗(yàn)設(shè)備與測試工況

    測試裝置主要采用東方所 INV3018CT 型采集儀、力錘(靈敏度為 0.197 mV/N)以及加速度傳感器,設(shè)置采樣頻率為2 000 Hz,傳感器布置及錘擊位置如圖1 所示.測試工況主要依據(jù)我國規(guī)范[13]中測試車輛準(zhǔn)靜態(tài)荷載作用下鋼軌扣件彈性墊板動(dòng)剛度的加載范圍選取,由于現(xiàn)有裝置的加載限制,本次試驗(yàn)主要選取30、40 、50 kN 作為施加在鋼軌上的垂向預(yù)壓力.實(shí)際運(yùn)行中列車作用在扣件膠墊上的預(yù)壓一般在0~50 kN 變化[14].在本次試驗(yàn)中,將千斤頂?shù)念A(yù)壓力取為30~50 kN,使該預(yù)壓力對(duì)鋼軌橫向振動(dòng)特性有較明顯的影響.

    圖1 鋼軌橫向加速度導(dǎo)納測試裝置Fig.1 Schematic of testing device of lateral acceleration admittance of rail

    試驗(yàn)?zāi)P椭胁捎冒惭b狀態(tài)下的單跨鋼軌,通過扣件系統(tǒng)安裝在混凝土基礎(chǔ)上,上部通過閉環(huán)鋼板頂住千斤頂.通過在理論模型中考慮梁的彎扭耦合特性和扣件膠墊采用空間模型的方法,模擬改變垂向預(yù)壓力對(duì)鋼軌橫向振動(dòng)特性的影響.試驗(yàn)中通過力錘橫向錘擊鋼軌軌頭,在上部千斤頂約束下,由于鋼軌橫截面質(zhì)心與剪切中心不重合使得鋼軌產(chǎn)生彎扭耦合作用.

    1.2 不同垂向預(yù)壓下鋼軌橫向加速度導(dǎo)納

    采用千斤頂給鋼軌施加垂向預(yù)壓力時(shí),測試的軌道結(jié)構(gòu)與周圍的鋼板框架形成整體振動(dòng),將使得測試結(jié)構(gòu)的剛度增大.為了消除外圍鋼板結(jié)構(gòu)的參振影響,采用本文裝置進(jìn)行試驗(yàn)時(shí),在千斤頂與鋼軌中間放置軟橡膠墊等可提供良好彈性的結(jié)構(gòu),且保證結(jié)構(gòu)具有一定厚度以達(dá)到較好的隔振效果.獲取加速度導(dǎo)納時(shí),確保有5 次有效錘擊結(jié)果,最終結(jié)果取其平均值.試驗(yàn)有效分析頻率的下限值由加速度導(dǎo)納的相干系數(shù)決定,本節(jié)分析段內(nèi)的相干系數(shù)均在0.8 以上.本次試驗(yàn)中錘擊力為瞬時(shí)定點(diǎn)動(dòng)荷載,通過傅里葉變換到頻域?yàn)楹喼C荷載,與下文理論部分所用的頻域激勵(lì)荷載一致.圖2 給出了垂向30 kN預(yù)壓下鋼軌橫向加速度導(dǎo)納5 次錘擊結(jié)果.觀察圖2 可知,30 kN 預(yù)壓下鋼軌橫向一階彎曲共振頻率為132.8 Hz.由于在40、50 kN 預(yù)壓下鋼軌5 次力錘敲擊均具有良好的測試結(jié)果.因此,僅展示5 次測試結(jié)果的平均值.垂向預(yù)壓力分別為30、40、50 kN 時(shí),錘擊得到的鋼軌橫向加速度導(dǎo)納平均值如圖3 所示.根據(jù)圖3 中多次力錘敲擊獲得的鋼軌橫向加速度導(dǎo)納平均值可知,在30、 40、50 kN垂向預(yù)壓下,鋼軌橫向一階彎曲共振頻率分別為132.8、140.6、146.5 Hz.

    圖2 垂向30 kN 預(yù)壓下鋼軌橫向加速度導(dǎo)納Fig.2 Lateral acceleration admittance of rail under vertical preload of 30 kN

    圖3 不同預(yù)壓下鋼軌橫向加速度導(dǎo)納測試結(jié)果平均值Fig.3 Average test results of lateral acceleration admittance of rail under different preloads

    2 考慮梁彎扭耦合的整體道床無砟軌道空間無限長模型

    為準(zhǔn)確高效求解鋼軌振動(dòng)特性,本節(jié)采用考慮彎扭耦合的Timoshenko 梁模型模擬鋼軌,采用波譜單元法(SEM)和辛數(shù)學(xué)方法[14]建立整體道床無砟軌道空間無限長模型.

    2.1 考慮彎扭耦合的Timoshenko 梁的波譜單元?jiǎng)觿偠染仃?/h3>

    首先建立考慮彎扭耦合的Timoshenko 梁動(dòng)力學(xué)運(yùn)動(dòng)方程:

    式中: δ(?) 為狄拉克函數(shù);m為單位長度質(zhì)量;Fy(x,t)、Fz(x,t)、T(x,t)分別為時(shí)刻t作用于軌梁上縱向x位置處的橫向力、垂向力和扭轉(zhuǎn)力矩;uy、uz分別為垂向和橫向位移;θx、θy、θz分別為扭轉(zhuǎn)轉(zhuǎn)角、橫向轉(zhuǎn)角和垂向轉(zhuǎn)角;Ky和Kz分別為橫向和垂向剪切修正因子;Iy、Iz和J分別為橫向、垂向截面慣性矩和截面扭轉(zhuǎn)常數(shù);A為橫截面面積;ρ為密度;I0為轉(zhuǎn)動(dòng)慣量;y0為鋼軌質(zhì)心至剪心的垂向距離;E、G分別為曲線梁的彈性模量和剪切模量;Kyp、Kzp和Kxr分別為扣件膠墊橫向、垂向和扭轉(zhuǎn)剛度;xrj為第j個(gè)扣件支點(diǎn)坐標(biāo);Nr為研究范圍內(nèi)扣件支點(diǎn)個(gè)數(shù).

    為分析鋼軌的振動(dòng)特性,首先將鋼軌運(yùn)動(dòng)方程進(jìn)行傅里葉變換到頻域,不考慮外荷載作用,得到頻域內(nèi)鋼軌的振動(dòng)方程.鋼軌自由波解存在以下形式:

    式中:U(?) 為時(shí)域位移;U1(?) 為頻域內(nèi)位移;U2(?)為波數(shù)域位移;ω為角頻率;k為波數(shù).

    由于鋼軌運(yùn)動(dòng)位移之間存在耦合關(guān)系,可令鋼軌頻域位移為

    式中:b、c為未知量;R1、T1、T2均為波模態(tài)振幅特征系數(shù).

    采用波數(shù)法求解自由梁的波數(shù)k,將式(6)、(7)代入式(1)~(5)可得鋼軌自由波動(dòng)方程為

    以梁的橫向運(yùn)動(dòng)為例,通過整理式(10)~(12)可得橫向彎曲波和扭轉(zhuǎn)波的6 個(gè)波數(shù),其特征方程為

    式中:X11=-mω2+GAKzk2;X12=-ikGAKz;X13=my0ω2;X21= -X12;X22=-ω2ρIy+EIyk2+GAKz;X23=0;X31=X13;X32=X23;X33=-ω2ρI0+GJk2.

    式(13)可寫為

    令X=0,可得自由梁的特征方程為

    式中:e1~e4為系數(shù).

    應(yīng)用范盛金公式[14]對(duì)式(15)進(jìn)行求解,可求解得到k的6 個(gè)取值k1~k6.根據(jù)波數(shù)ki可計(jì)算式(13)中對(duì)應(yīng)的模態(tài)系數(shù)T1、T2.

    將求得的波數(shù)和對(duì)應(yīng)的模態(tài)系數(shù)代入式(7),便可得到自由梁頻域-波數(shù)域位移.同理,可得垂向位移和垂向轉(zhuǎn)角的表達(dá)式.結(jié)合Timoshenko 梁的位移邊界條件可消除位移中的系數(shù)b和c,得到頻域內(nèi)的位移.

    求得自由梁波數(shù)后,便可根據(jù)鋼軌頻域運(yùn)動(dòng)方程求解梁波譜動(dòng)剛度矩陣,進(jìn)而可利用動(dòng)剛度矩陣法求解梁的位移.為準(zhǔn)確方便獲得波譜動(dòng)剛度矩陣,采用變分法對(duì)梁頻域-波數(shù)域位移進(jìn)行變分.對(duì)橫向和扭轉(zhuǎn)運(yùn)動(dòng)方程進(jìn)行變分,可得其弱積分形式為

    式中:l為扣件間距.

    根據(jù)材料力學(xué)中梁內(nèi)部荷載與位移的關(guān)系,通過分步積分,可得降階后的變分方程,進(jìn)而將消除系數(shù)的位移表達(dá)式代入式(16),可得結(jié)構(gòu)橫向和扭轉(zhuǎn)位移的譜單元方程為

    式中:Dc(ω)=GAKz(-KcEzKc+iKcEzT1)-EIyT1KcEz×KcT1+GAKz(T1EzT1+iT1EzKc)-GJT2KcEzKcT2-ω2m×Ez-ω2(ρIyT1EzT1-my0EzT2+ρI0T2EzT2-my0T2Ez)=,Kc=diag[k1k2···k6],Ez=(x,ω)ez(x,ω)×dz=,Eci j如式(18),ez(x,ω)=diag[e-ik1xe-ik2x···e-ik6x] ;Hc為根據(jù)位移邊界條件得到的用于消除未知系數(shù)c的矩陣.

    式中:a,b=1,2,···6.

    同理,可用同樣的方法得到縱向和垂向運(yùn)動(dòng)方程的譜動(dòng)剛度矩陣Su.

    2.2 整體道床軌道空間無限長模型的建立

    本節(jié)基于上述推導(dǎo)的考慮梁彎扭耦合的波譜單元?jiǎng)觿偠染仃?,建立考慮梁彎扭耦合的整體道床無砟軌道空間無限長模型.由于整體道床及其下部基礎(chǔ)剛度很大,整體道床無砟軌道可視為周期離散點(diǎn)支承的鋼軌-扣件系統(tǒng),如圖4 所示,圖中:Cyp為扣件膠墊阻尼;I為截面慣性矩,因此,能夠應(yīng)用辛數(shù)學(xué)方法[15]等求解該類軌道結(jié)構(gòu)的動(dòng)力響應(yīng).鋼軌-扣件周期子系統(tǒng)的動(dòng)剛度矩陣Kt如式(19),由扣件系統(tǒng)動(dòng)剛度矩陣Kp(式(20)與鋼軌波譜單元?jiǎng)觿偠染仃嘖g(式(21))構(gòu)成,其中,單個(gè)子結(jié)構(gòu)的鋼軌動(dòng)剛度矩陣是由單元1 和單元2 動(dòng)剛度矩陣集成的,如圖4 中子結(jié)構(gòu)所示.

    圖4 鋼軌-扣件系統(tǒng)的整體結(jié)構(gòu)與周期子結(jié)構(gòu)示意Fig.4 Overall structure and periodic substructure of a rail-fastener system

    式中:Crp和Krp分別為膠墊的阻尼矩陣和剛度矩陣;分別為結(jié)構(gòu)單元1、2 的鋼軌剛度矩陣,=Su(ω)Pv+Sc(ω)Pl,Pv和Pl分別為使譜剛度矩陣Su和Sc擴(kuò)展到整體矩陣的轉(zhuǎn)換矩陣,式同.

    式(21)中不考慮扣件膠墊的質(zhì)量,扣件膠墊作用于鋼軌的物理模型如圖5.圖中:C為質(zhì)心;S為剪切中心;ωF為激勵(lì)圓頻率;Czp為橫向阻尼系數(shù);hr2為軌底距離;Kzp為橫向剛度;br2為底寬的1/2.

    圖5 鋼軌橫截面示意Fig.5 Schematic of rail cross section

    由于采用辛數(shù)學(xué)理論[15]求解周期結(jié)構(gòu)的方法是成熟的,這里將不再贅述.

    2.3 模型正確性驗(yàn)證

    本節(jié)采用與現(xiàn)有文獻(xiàn)計(jì)算結(jié)果進(jìn)行對(duì)比來驗(yàn)證上述提出的模型和計(jì)算方法的正確性.根據(jù)文獻(xiàn)[3]中的參數(shù),采用波譜-辛混合法計(jì)算有砟軌道結(jié)構(gòu)的鋼軌速度導(dǎo)納.由于文獻(xiàn)[3]中采用的梁未考慮彎扭耦合,應(yīng)用上述模型計(jì)算時(shí),將y0取為無窮小,即10-7,再與文獻(xiàn)[3]計(jì)算結(jié)果進(jìn)行對(duì)比.

    本文與文獻(xiàn)[3]關(guān)于鋼軌跨中垂向速度導(dǎo)納的計(jì)算結(jié)果對(duì)比如圖6 所示.

    圖6 鋼軌跨中垂向速度導(dǎo)納Fig.6 Vertical velocity admittance of a rail at the mid-span

    由圖6 可知,本文采用的波譜-辛混合法計(jì)算的結(jié)果和文獻(xiàn)中采用解析法所得結(jié)果基本一致,驗(yàn)證了本文采用的模型和方法的可靠性.

    3 膠墊垂向預(yù)壓特性對(duì)鋼軌橫向振動(dòng)特性的影響規(guī)律

    分析考慮梁彎扭耦合與否對(duì)鋼軌橫向和扭轉(zhuǎn)固有頻率的影響規(guī)律,并與傳統(tǒng)有限單元法(FEM)的計(jì)算結(jié)果進(jìn)行對(duì)比分析,進(jìn)一步驗(yàn)證本文建立模型的正確性;采用文獻(xiàn)[16]測得的不同預(yù)壓下扣件膠墊動(dòng)參數(shù),計(jì)算分析其對(duì)鋼軌橫向速度導(dǎo)納的影響規(guī)律,并對(duì)比分析仿真的與試驗(yàn)獲得的不同預(yù)壓下鋼軌橫向一階彎曲共振頻率.

    3.1 考慮梁彎扭耦合對(duì)Timoshenko 梁橫向和扭轉(zhuǎn)固有頻率的影響分析

    為便于應(yīng)用有限元軟件進(jìn)行模態(tài)分析,將鋼軌視為一根兩端簡支支撐的Timoshenko 梁,選取總長為5.9 m 的梁進(jìn)行計(jì)算分析,該長度已足夠有效計(jì)算鋼軌前20 階固有頻率.其幾何和材料屬性可參考文獻(xiàn)[2]給出的T60 軌的相關(guān)參數(shù),軌道結(jié)構(gòu)其他參數(shù)如表1 所示.

    表1 軌道結(jié)構(gòu)參數(shù)Tab.1 Structural parameters of rail

    為了獲得簡支梁固有頻率,首先,對(duì)式(21)獲得的動(dòng)剛度矩陣施加相應(yīng)的簡支梁邊界條件使其行列式為0,然后,采用Wittrick-Williams 算法[17]搜索其根值,繼而得到準(zhǔn)確的固有頻率,如表2 所示.采用FEM 進(jìn)行計(jì)算時(shí),Beam188 梁單元采用鋼軌實(shí)體截面,由于采用不同的單元數(shù)對(duì)計(jì)算結(jié)果有較大的影響,這里給出了單元數(shù)n分別取20、100、1 000 時(shí)簡支梁的固有頻率(表2).

    表2 Timoshenko 簡支梁橫向彎曲和扭轉(zhuǎn)固有頻率Tab.2 Lateral bending and torsional natural frequencies of a simply-supported Timoshenko beamHz

    從表2 中可以看出:FEM 采用單元數(shù)為1 000時(shí)與SEM 采用一個(gè)單元計(jì)算得到的固有頻率相近.為了提高計(jì)算效率,有限單元數(shù)為100 時(shí)便可得到相對(duì)準(zhǔn)確的結(jié)果.考慮梁的彎扭耦合振動(dòng)主要使得鋼軌3 階以上模態(tài)的固有頻率增大,而其中一些偶數(shù)階模態(tài)的固有頻率減小.

    3.2 考慮梁彎扭耦合對(duì)鋼軌橫向振動(dòng)特性的影響分析

    扣件膠墊采用分?jǐn)?shù)階Zener 模型模擬[16].為分析考慮Timoshenko 梁彎扭耦合對(duì)鋼軌振動(dòng)特性的影響,圖7 給出了橫向單位簡諧激勵(lì)荷載下考慮梁彎扭耦合與否的鋼軌橫向速度導(dǎo)納和扭轉(zhuǎn)速度交叉導(dǎo)納.圖中實(shí)線為梁彎扭解耦條件(情況1)下的橫向和扭轉(zhuǎn)速度導(dǎo)納,虛線為考慮梁彎扭耦合變形(情況2)的結(jié)果.

    圖7 考慮梁彎扭耦合對(duì)鋼軌橫向和扭轉(zhuǎn)速度導(dǎo)納的影響Fig.7 Influence of vertical preload dependence of fastener pads on lateral velocity admittance of rail

    考慮梁彎扭耦合對(duì)鋼軌橫向和扭轉(zhuǎn)速度導(dǎo)納有相似的影響,使其第1 個(gè)峰(圖中記為fl)的頻率從110.1 Hz 增加到139.6 Hz,增加了大約29.6 Hz,即鋼軌相對(duì)扣件子系統(tǒng)的橫向彎曲共振頻率,峰幅值略微增加.同時(shí),從表2 中可以看出,當(dāng)鋼軌固有頻率接近fl時(shí),考慮梁彎扭耦合的鋼軌固有頻率大于梁彎扭解耦時(shí)對(duì)應(yīng)的值.然而,對(duì)于橫向速度導(dǎo)納,其第2 個(gè)峰值(圖中記為flp1)的頻率減小了大約47.0 Hz,即鋼軌橫向一階pinned-pinned 共振頻率,且該pinned-pinned 模態(tài)出現(xiàn)在扭轉(zhuǎn)速度導(dǎo)納曲線上,這是由于考慮了梁的彎扭耦合作用.同時(shí),從橫向速度導(dǎo)納曲線可以發(fā)現(xiàn),考慮梁彎扭耦合,其具有扭轉(zhuǎn)一階和二階pinned-pinned 模態(tài)引起的峰值(圖中記為ftp1和ftp2).考慮梁彎扭耦合對(duì)鋼軌橫向二階pinnedpinned 共振頻率(圖中記為flp2)影響較小.

    下文通過結(jié)合實(shí)測與計(jì)算的不同垂向預(yù)壓下鋼軌橫向一階彎曲共振頻率,說明考慮梁彎扭耦合的必要性.

    3.3 膠墊垂向預(yù)壓特性對(duì)鋼軌橫向彎曲振動(dòng)特性的影響分析

    不同垂向預(yù)壓下,在橫向單位簡諧荷載激勵(lì)時(shí)鋼軌跨中軌頭原點(diǎn)橫向速度導(dǎo)納如圖8 所示.由圖可知,考慮扣件膠墊預(yù)壓特性主要影響鋼軌橫向中低頻振動(dòng),即橫向一階pinned-pinned 共振頻率以下,而對(duì)其pinned-pinned 共振頻率幾乎沒有影響.隨著垂向預(yù)壓力的增大,鋼軌橫向彎曲共振頻率增大,而該頻率以下的振動(dòng)幅值略微降低.

    圖8 扣件膠墊垂向預(yù)壓特性對(duì)鋼軌橫向速度導(dǎo)納的影響Fig.8 Influence of the vertical preload dependence of rail pads on the lateral mobility of rail

    與考慮梁彎扭耦合相比,不考慮梁彎扭耦合時(shí)鋼軌具有更高的橫向彎曲共振頻率,為詳細(xì)分析不同垂向預(yù)壓下考慮梁彎扭耦合與否仿真的和實(shí)測的鋼軌橫向一階彎曲共振頻率的關(guān)系,圖9 給出了三者隨預(yù)壓變化的取值.

    圖9 不同預(yù)壓下鋼軌橫向彎曲共振頻率的實(shí)測與仿真值Fig.9 Measured and simulated results of lateral BRF of rail under different preloads

    由圖9 可知,隨著預(yù)壓的增大,考慮梁彎扭耦合得到的鋼軌橫向共振頻率變化規(guī)律與實(shí)測的結(jié)果更為接近.由于鋼軌扭轉(zhuǎn)共振頻率高于橫向彎曲共振頻率,考慮梁彎扭耦合特性使得鋼軌橫向彎曲共振頻率高于未考慮的.受扭轉(zhuǎn)振動(dòng)的約束,隨著預(yù)壓的增大,考慮梁彎扭耦合使得鋼軌橫向彎曲共振頻率變化范圍較未考慮的減小.當(dāng)預(yù)壓從30 kN 增加到50 kN時(shí),實(shí)測的橫向彎曲共振頻率增加了13.7 Hz,考慮梁彎扭耦合時(shí)其增加了12.5 Hz,而未考慮時(shí)其增加了21.7 Hz.因此,考慮梁彎扭耦合能更準(zhǔn)確地預(yù)測鋼軌的振動(dòng)特性.

    4 結(jié) 論

    本文主要研究結(jié)論如下:

    1) 本文提出的波譜-辛混合法具有較好的計(jì)算精度和效率.SEM 采用1 個(gè)單元計(jì)算得到的固有頻率值與傳統(tǒng)有限元采用單元數(shù)為1 000 時(shí)更相近.

    2) 通過錘擊試驗(yàn)發(fā)現(xiàn),隨著垂向預(yù)壓力從30 kN增加到50 kN 時(shí),測得的鋼軌橫向一階彎曲共振頻率從132.8 Hz 增加到146.5 Hz.

    3) 考慮梁彎扭耦合使得鋼軌橫向彎曲共振頻率增加了大約29.6 Hz,而其一階pinned-pinned 共振頻率明顯減小了大約47.0 Hz.同時(shí),鋼軌橫向彎曲和扭轉(zhuǎn)pinned-pinned 模態(tài)相互參與彼此的振動(dòng).因此,鋼軌橫向和扭轉(zhuǎn)振動(dòng)頻譜更豐富.

    4) 扣件膠墊垂向預(yù)壓特性主要影響鋼軌橫向中低頻振動(dòng),隨著預(yù)壓的增大鋼軌橫向彎曲共振頻率增大.當(dāng)預(yù)壓從30 kN 增加到50 kN 時(shí),實(shí)測的橫向彎曲共振頻率增加了13.7 Hz,考慮梁彎扭耦合時(shí)其增加了12.5 Hz,而未考慮時(shí)其增加了21.7 Hz.因此,為準(zhǔn)確預(yù)測鋼軌的振動(dòng)響應(yīng),對(duì)軌道結(jié)構(gòu)進(jìn)行數(shù)值仿真時(shí)需考慮梁的彎扭耦合作用.

    猜你喜歡
    膠墊共振頻率扣件
    多激勵(lì)下結(jié)構(gòu)反共振頻率的計(jì)算和配置
    科隆蛋扣件力學(xué)特性分析
    一種改進(jìn)的LBP特征實(shí)現(xiàn)鐵路扣件識(shí)別
    Multi-functional roles of TaSSI2 involved in Fusarium head blight and powdery mildew resistance and drought tolerance
    一種基于工業(yè)機(jī)器人的膠墊自動(dòng)粘貼系統(tǒng)的研究與開發(fā)
    換膠墊250
    非均布荷載下扣件式腳手架的隨機(jī)缺陷穩(wěn)定分析
    三層軟磁鎳超晶格薄膜的共振頻率*
    WJ-8型小阻力扣件軌下橡膠墊板滑出動(dòng)力學(xué)研究
    WJ -8型小阻力扣件軌下膠墊滑出縱向阻力試驗(yàn)研究
    2021天堂中文幕一二区在线观| 成年免费大片在线观看| 一本色道久久久久久精品综合| 午夜福利在线在线| 久久热精品热| av.在线天堂| 在线天堂最新版资源| 一个人看视频在线观看www免费| 精品视频人人做人人爽| 老女人水多毛片| 乱码一卡2卡4卡精品| 国产亚洲av片在线观看秒播厂| 新久久久久国产一级毛片| 在线亚洲精品国产二区图片欧美 | 亚洲精品第二区| 亚洲国产色片| 亚洲国产色片| 狂野欧美激情性bbbbbb| 秋霞在线观看毛片| 久久久久国产精品人妻一区二区| 五月天丁香电影| 下体分泌物呈黄色| 亚洲va在线va天堂va国产| 插阴视频在线观看视频| 亚洲精品日本国产第一区| 乱系列少妇在线播放| 久久久色成人| 中文乱码字字幕精品一区二区三区| 亚洲精品国产av成人精品| 婷婷色综合www| 老司机影院成人| 国产精品99久久99久久久不卡 | 噜噜噜噜噜久久久久久91| 青春草亚洲视频在线观看| 日本av手机在线免费观看| 六月丁香七月| 18禁裸乳无遮挡免费网站照片| 99久国产av精品国产电影| 国产精品熟女久久久久浪| 在线免费十八禁| 亚洲欧美中文字幕日韩二区| 能在线免费看毛片的网站| 禁无遮挡网站| 狂野欧美激情性bbbbbb| 天堂俺去俺来也www色官网| 国产一级毛片在线| 国产一级毛片在线| 欧美高清性xxxxhd video| 极品少妇高潮喷水抽搐| 女人被狂操c到高潮| 新久久久久国产一级毛片| 老女人水多毛片| 99热这里只有是精品50| 精品久久久久久久人妻蜜臀av| 国产精品人妻久久久久久| 亚洲性久久影院| 成年版毛片免费区| 国产男女超爽视频在线观看| 精品国产一区二区三区久久久樱花 | 一区二区三区精品91| 哪个播放器可以免费观看大片| 国产午夜精品久久久久久一区二区三区| 日日啪夜夜撸| 亚洲美女搞黄在线观看| 亚洲国产欧美在线一区| 国产av不卡久久| 在线观看国产h片| 啦啦啦中文免费视频观看日本| 久久午夜福利片| 一本一本综合久久| 麻豆乱淫一区二区| 极品教师在线视频| 国产大屁股一区二区在线视频| 亚洲精品乱码久久久v下载方式| 亚洲天堂国产精品一区在线| 99久久精品一区二区三区| 大香蕉久久网| 一级毛片黄色毛片免费观看视频| 亚洲精品aⅴ在线观看| 国产黄频视频在线观看| 国产精品av视频在线免费观看| 高清av免费在线| 中国国产av一级| 久久久久九九精品影院| 下体分泌物呈黄色| 又爽又黄无遮挡网站| 直男gayav资源| 又黄又爽又刺激的免费视频.| 欧美高清成人免费视频www| 国产成人freesex在线| 国产熟女欧美一区二区| 交换朋友夫妻互换小说| 国产黄频视频在线观看| 国产精品一及| 听说在线观看完整版免费高清| 超碰av人人做人人爽久久| 91久久精品国产一区二区成人| 国产免费视频播放在线视频| 成年人午夜在线观看视频| 国产精品一二三区在线看| 午夜精品国产一区二区电影 | 婷婷色av中文字幕| 中文字幕免费在线视频6| 伊人久久国产一区二区| 久久精品熟女亚洲av麻豆精品| 亚洲精华国产精华液的使用体验| h日本视频在线播放| 纵有疾风起免费观看全集完整版| 人人妻人人澡人人爽人人夜夜| 少妇的逼水好多| 久久精品国产鲁丝片午夜精品| 99热全是精品| 国语对白做爰xxxⅹ性视频网站| 欧美少妇被猛烈插入视频| 午夜日本视频在线| 国产成人freesex在线| 真实男女啪啪啪动态图| 99热全是精品| 亚洲精品国产av成人精品| 色视频在线一区二区三区| 精品久久久久久电影网| 神马国产精品三级电影在线观看| 岛国毛片在线播放| 亚洲人成网站在线观看播放| 久久精品国产亚洲av天美| 水蜜桃什么品种好| 亚洲高清免费不卡视频| 有码 亚洲区| 视频区图区小说| 22中文网久久字幕| 亚洲经典国产精华液单| 亚洲经典国产精华液单| 精品人妻偷拍中文字幕| 97人妻精品一区二区三区麻豆| 国产成人91sexporn| 久久久久久国产a免费观看| 欧美激情在线99| 亚洲欧美中文字幕日韩二区| 在现免费观看毛片| 国产精品久久久久久精品古装| 极品少妇高潮喷水抽搐| 日韩,欧美,国产一区二区三区| 久久久久久久久久人人人人人人| 男的添女的下面高潮视频| 97人妻精品一区二区三区麻豆| 国产欧美另类精品又又久久亚洲欧美| 偷拍熟女少妇极品色| 婷婷色综合www| 啦啦啦啦在线视频资源| 观看免费一级毛片| 日日摸夜夜添夜夜添av毛片| 国产黄片美女视频| 午夜日本视频在线| 日韩制服骚丝袜av| 国产亚洲午夜精品一区二区久久 | 国产一区有黄有色的免费视频| 国产伦理片在线播放av一区| 一区二区三区免费毛片| a级毛色黄片| 最近最新中文字幕大全电影3| 中文资源天堂在线| 国产老妇女一区| 亚洲真实伦在线观看| 国产成人91sexporn| 综合色丁香网| 毛片女人毛片| 国产一级毛片在线| 久热这里只有精品99| 国产有黄有色有爽视频| av播播在线观看一区| 中国三级夫妇交换| 美女xxoo啪啪120秒动态图| 国产成人a∨麻豆精品| 国产成人aa在线观看| 国产亚洲av片在线观看秒播厂| 人人妻人人爽人人添夜夜欢视频 | 欧美精品国产亚洲| 亚洲色图综合在线观看| 在线免费观看不下载黄p国产| 成年av动漫网址| 80岁老熟妇乱子伦牲交| 三级男女做爰猛烈吃奶摸视频| 神马国产精品三级电影在线观看| 日韩大片免费观看网站| 少妇裸体淫交视频免费看高清| 插阴视频在线观看视频| 免费在线观看成人毛片| 久久精品国产亚洲av天美| 免费看不卡的av| 日本黄大片高清| 99热国产这里只有精品6| 一级毛片电影观看| 一个人观看的视频www高清免费观看| 青春草国产在线视频| 美女被艹到高潮喷水动态| 国产av码专区亚洲av| 国产午夜精品久久久久久一区二区三区| xxx大片免费视频| 91精品伊人久久大香线蕉| 五月天丁香电影| 国产乱人偷精品视频| 午夜视频国产福利| 国产免费福利视频在线观看| 亚洲美女视频黄频| 校园人妻丝袜中文字幕| 亚洲人与动物交配视频| 在线观看一区二区三区激情| 成人亚洲精品一区在线观看 | 一级爰片在线观看| 亚洲成色77777| 亚洲精品,欧美精品| 又粗又硬又长又爽又黄的视频| 一区二区三区四区激情视频| 狂野欧美白嫩少妇大欣赏| 在线观看人妻少妇| 亚洲经典国产精华液单| 亚洲精品久久久久久婷婷小说| 久久久久久久久大av| 91在线精品国自产拍蜜月| 麻豆成人午夜福利视频| 国产成人a区在线观看| 老司机影院毛片| 欧美日韩精品成人综合77777| 搡老乐熟女国产| 精品人妻熟女av久视频| 亚洲人成网站高清观看| .国产精品久久| 少妇被粗大猛烈的视频| 一级二级三级毛片免费看| 狂野欧美白嫩少妇大欣赏| 国产精品麻豆人妻色哟哟久久| 九草在线视频观看| 国产乱人视频| 最近最新中文字幕免费大全7| 欧美日韩国产mv在线观看视频 | 久久人人爽人人片av| 成年女人在线观看亚洲视频 | 一区二区三区乱码不卡18| 欧美 日韩 精品 国产| 亚洲无线观看免费| 天天躁夜夜躁狠狠久久av| 亚洲精品色激情综合| 六月丁香七月| 亚洲人成网站在线播| 十八禁网站网址无遮挡 | 亚洲精品成人av观看孕妇| 成人美女网站在线观看视频| 国产午夜福利久久久久久| 久久久成人免费电影| 综合色丁香网| 美女高潮的动态| 赤兔流量卡办理| 汤姆久久久久久久影院中文字幕| a级毛片免费高清观看在线播放| 在线免费十八禁| 中文字幕久久专区| 寂寞人妻少妇视频99o| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 一本久久精品| 直男gayav资源| 在线a可以看的网站| 一区二区av电影网| 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| 搞女人的毛片| av播播在线观看一区| 日韩伦理黄色片| 大码成人一级视频| 青春草国产在线视频| 建设人人有责人人尽责人人享有的 | 国产人妻一区二区三区在| 乱系列少妇在线播放| 成年免费大片在线观看| 欧美日韩亚洲高清精品| 国产久久久一区二区三区| 国产亚洲av片在线观看秒播厂| 精品久久久噜噜| 另类亚洲欧美激情| 亚洲在久久综合| 免费观看av网站的网址| 寂寞人妻少妇视频99o| 22中文网久久字幕| av国产免费在线观看| 久久久精品94久久精品| 简卡轻食公司| 日韩一区二区视频免费看| 我的老师免费观看完整版| 久久人人爽人人爽人人片va| 国产精品一区www在线观看| 一级爰片在线观看| 两个人的视频大全免费| 五月玫瑰六月丁香| av女优亚洲男人天堂| 国产精品麻豆人妻色哟哟久久| 成年女人看的毛片在线观看| 亚洲精品色激情综合| 热99国产精品久久久久久7| 丝袜喷水一区| 亚洲欧美精品自产自拍| 精品午夜福利在线看| 18禁在线无遮挡免费观看视频| 热re99久久精品国产66热6| 欧美激情国产日韩精品一区| 街头女战士在线观看网站| 青青草视频在线视频观看| 热re99久久精品国产66热6| 在线精品无人区一区二区三 | 久久精品人妻少妇| 狂野欧美激情性xxxx在线观看| 亚洲欧美日韩另类电影网站 | 久久精品国产亚洲网站| 另类亚洲欧美激情| 国产黄色免费在线视频| 成人漫画全彩无遮挡| 又爽又黄无遮挡网站| 十八禁网站网址无遮挡 | 成年av动漫网址| 99视频精品全部免费 在线| 久久久久久伊人网av| 国语对白做爰xxxⅹ性视频网站| 乱码一卡2卡4卡精品| 欧美区成人在线视频| 各种免费的搞黄视频| 成人无遮挡网站| 精品亚洲乱码少妇综合久久| 国产白丝娇喘喷水9色精品| 国产成人免费观看mmmm| 高清欧美精品videossex| 亚洲av在线观看美女高潮| 小蜜桃在线观看免费完整版高清| 国产男女超爽视频在线观看| 国产高潮美女av| 中文字幕av成人在线电影| 欧美日韩精品成人综合77777| 亚洲av免费在线观看| 18禁在线播放成人免费| 毛片一级片免费看久久久久| 亚洲人与动物交配视频| 欧美3d第一页| 老师上课跳d突然被开到最大视频| 成人午夜精彩视频在线观看| 黄色一级大片看看| 大话2 男鬼变身卡| 黄片wwwwww| 九色成人免费人妻av| 91久久精品国产一区二区成人| 亚洲人与动物交配视频| 日本三级黄在线观看| 免费看不卡的av| 欧美精品一区二区大全| 久久精品久久久久久噜噜老黄| 国产av不卡久久| 嫩草影院入口| 久久久久性生活片| av线在线观看网站| 少妇的逼水好多| 精品国产一区二区三区久久久樱花 | 亚洲一级一片aⅴ在线观看| 亚洲成人久久爱视频| 高清视频免费观看一区二区| 26uuu在线亚洲综合色| 欧美成人a在线观看| 中文字幕久久专区| freevideosex欧美| 国产黄片视频在线免费观看| 内地一区二区视频在线| 大陆偷拍与自拍| 亚洲国产欧美人成| 欧美丝袜亚洲另类| 国产精品不卡视频一区二区| 亚洲熟女精品中文字幕| 一区二区三区免费毛片| 人妻 亚洲 视频| 99热国产这里只有精品6| 一本一本综合久久| 中国美白少妇内射xxxbb| 视频区图区小说| 久久久久网色| 高清视频免费观看一区二区| 成人国产麻豆网| 亚洲av免费高清在线观看| 日韩欧美一区视频在线观看 | 日日啪夜夜爽| 亚洲成人精品中文字幕电影| 亚洲美女搞黄在线观看| 日本熟妇午夜| 免费不卡的大黄色大毛片视频在线观看| 乱系列少妇在线播放| 特级一级黄色大片| 日日撸夜夜添| 久久久久国产精品人妻一区二区| 日韩国内少妇激情av| 国产精品人妻久久久影院| 秋霞在线观看毛片| 女人久久www免费人成看片| 啦啦啦在线观看免费高清www| 国产精品久久久久久精品古装| 国产av码专区亚洲av| 嫩草影院精品99| 亚洲不卡免费看| 午夜福利高清视频| 毛片一级片免费看久久久久| 日韩电影二区| 精品国产乱码久久久久久小说| 国产老妇伦熟女老妇高清| 亚洲在久久综合| 亚洲伊人久久精品综合| 欧美xxxx性猛交bbbb| 免费大片黄手机在线观看| 亚洲av福利一区| 欧美3d第一页| 欧美三级亚洲精品| 美女xxoo啪啪120秒动态图| 国产 一区 欧美 日韩| 亚洲av国产av综合av卡| 69av精品久久久久久| 另类亚洲欧美激情| 视频区图区小说| 免费黄色在线免费观看| 国产白丝娇喘喷水9色精品| 亚洲国产欧美人成| 亚洲美女搞黄在线观看| 国产永久视频网站| av卡一久久| 国产毛片在线视频| 国产男人的电影天堂91| 嘟嘟电影网在线观看| 久久热精品热| 亚洲精品影视一区二区三区av| 国产又色又爽无遮挡免| 免费av毛片视频| 美女视频免费永久观看网站| 日本黄大片高清| 啦啦啦中文免费视频观看日本| 联通29元200g的流量卡| 99热国产这里只有精品6| 午夜福利视频1000在线观看| 日日摸夜夜添夜夜添av毛片| 成人国产av品久久久| 99久久精品一区二区三区| 久久精品国产自在天天线| 久久精品夜色国产| 欧美激情久久久久久爽电影| 国产成人福利小说| 成人亚洲精品av一区二区| 日本黄大片高清| 美女内射精品一级片tv| 麻豆成人av视频| 各种免费的搞黄视频| 91精品国产九色| freevideosex欧美| 99久久精品国产国产毛片| 伊人久久精品亚洲午夜| 国产亚洲午夜精品一区二区久久 | 欧美bdsm另类| 亚洲精品第二区| 亚洲av日韩在线播放| 欧美成人a在线观看| 亚洲性久久影院| 一边亲一边摸免费视频| 亚洲av在线观看美女高潮| 中文欧美无线码| 在现免费观看毛片| av福利片在线观看| av在线播放精品| 久久亚洲国产成人精品v| 爱豆传媒免费全集在线观看| 欧美日韩亚洲高清精品| 男人舔奶头视频| 毛片女人毛片| 网址你懂的国产日韩在线| 成人黄色视频免费在线看| 深夜a级毛片| 国产精品99久久99久久久不卡 | 欧美3d第一页| 亚洲欧美精品自产自拍| 性插视频无遮挡在线免费观看| 亚洲无线观看免费| 亚洲国产欧美人成| 最近最新中文字幕大全电影3| 91aial.com中文字幕在线观看| 热99国产精品久久久久久7| 波野结衣二区三区在线| 搡老乐熟女国产| 日韩视频在线欧美| 菩萨蛮人人尽说江南好唐韦庄| 干丝袜人妻中文字幕| 久久久久久九九精品二区国产| 亚洲精品国产av成人精品| 99久久九九国产精品国产免费| 综合色av麻豆| 久久精品熟女亚洲av麻豆精品| 高清日韩中文字幕在线| 在线亚洲精品国产二区图片欧美 | 日本wwww免费看| 国产一级毛片在线| 蜜臀久久99精品久久宅男| 水蜜桃什么品种好| 国产精品一区二区性色av| 国语对白做爰xxxⅹ性视频网站| 国产毛片a区久久久久| 青春草亚洲视频在线观看| 国产欧美日韩精品一区二区| 99久久中文字幕三级久久日本| 欧美成人午夜免费资源| 精品久久久久久久久亚洲| 免费在线观看成人毛片| 国内精品美女久久久久久| 一个人观看的视频www高清免费观看| 男女那种视频在线观看| 国产成人一区二区在线| 久久99精品国语久久久| 一个人观看的视频www高清免费观看| 国产大屁股一区二区在线视频| 91久久精品电影网| 欧美国产精品一级二级三级 | 久久精品国产a三级三级三级| 日韩欧美 国产精品| 在线观看国产h片| 国产乱人偷精品视频| 男女边摸边吃奶| 舔av片在线| 别揉我奶头 嗯啊视频| 交换朋友夫妻互换小说| 日本av手机在线免费观看| 国产综合懂色| 久久99精品国语久久久| videossex国产| 另类亚洲欧美激情| 亚洲怡红院男人天堂| 久久久久久久大尺度免费视频| 成人黄色视频免费在线看| 亚洲成人一二三区av| 国产精品成人在线| 国产69精品久久久久777片| 亚洲成人精品中文字幕电影| 国产乱来视频区| 日产精品乱码卡一卡2卡三| 天堂俺去俺来也www色官网| 日产精品乱码卡一卡2卡三| 国产精品一二三区在线看| 91aial.com中文字幕在线观看| 亚洲成人一二三区av| 少妇人妻一区二区三区视频| 久久久久国产精品人妻一区二区| 91午夜精品亚洲一区二区三区| 亚洲国产精品成人久久小说| 永久网站在线| 大香蕉97超碰在线| 国产av码专区亚洲av| 亚洲国产精品成人综合色| 99re6热这里在线精品视频| 高清日韩中文字幕在线| 国产精品国产三级国产专区5o| 欧美区成人在线视频| 有码 亚洲区| 欧美高清性xxxxhd video| 国产精品一区二区性色av| 国产视频首页在线观看| 内地一区二区视频在线| 日韩av不卡免费在线播放| 国产一区亚洲一区在线观看| av在线老鸭窝| 国产欧美另类精品又又久久亚洲欧美| 日日啪夜夜爽| 亚州av有码| 美女被艹到高潮喷水动态| 欧美少妇被猛烈插入视频| 免费av观看视频| 啦啦啦中文免费视频观看日本| 三级经典国产精品| 一区二区三区乱码不卡18| 你懂的网址亚洲精品在线观看| 热99国产精品久久久久久7| 两个人的视频大全免费| 乱系列少妇在线播放| 中文字幕av成人在线电影| 国产精品国产三级国产专区5o| 成人黄色视频免费在线看| 国产91av在线免费观看| 日日撸夜夜添| 免费少妇av软件| 欧美少妇被猛烈插入视频| 欧美老熟妇乱子伦牲交| 久久精品久久精品一区二区三区| 看十八女毛片水多多多| 亚洲综合精品二区| 亚洲精品视频女| 欧美成人精品欧美一级黄| 99热网站在线观看| 久久这里有精品视频免费| 久久精品久久精品一区二区三区| 成人亚洲欧美一区二区av| 亚洲欧美中文字幕日韩二区| 久久99蜜桃精品久久| 欧美日韩亚洲高清精品| 美女国产视频在线观看| 国产色婷婷99| 嫩草影院精品99| 黄色一级大片看看| 欧美日韩国产mv在线观看视频 | 熟女电影av网| 男人和女人高潮做爰伦理| 日韩av在线免费看完整版不卡| 三级国产精品片| 久久精品人妻少妇| 老司机影院成人| h日本视频在线播放| 久久久久久久精品精品| 赤兔流量卡办理| 精品久久久久久久久av| 免费少妇av软件| 黑人高潮一二区| 在线观看人妻少妇| 日韩制服骚丝袜av|