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

    有砟軌道結構彈性波傳播特性研究

    2019-08-02 03:20:34趙才友
    鐵道學報 2019年6期
    關鍵詞:周期性扣件鋼軌

    易 強,王 平,趙才友,盛 曦,盧 俊

    (1.西南交通大學高速鐵路線路工程教育部重點實驗室,四川成都 610031;2.西南交通大學土木工程學院,四川成都 610031)

    鐵路軌道結構直線段通常設計為相同的基本單元沿線路縱向規(guī)則排列,對于有砟軌道,基本單元由一個扣件間距內鋼軌、扣件系統(tǒng)、軌枕、道床以及下部基礎組成,因此鐵路軌道結構可視為周期結構。近代固體物理學研究表明周期結構具有重要的物理特性,即彈性波帶隙特性。彈性波在周期結構中傳播時會產生一系列特殊的物理特性,周期結構具有能夠抑制彈性波傳播的頻率范圍,這些頻率范圍稱為帶隙[1]。另一方面,結構振動噪聲產生的本質原因通??蓺w結為結構中彈性波傳播效應以及結構彈性波與周圍介質的相互耦合作用[2],因此分析研究彈性波在軌道結構中的傳播特性對于軌道結構振動噪聲控制具有重要的意義。固體物理領域提出的聲子晶體概念為周期結構中彈性波的研究工作提供了新的研究思路,聲子晶體是一類由特殊設計的人工結構單元周期排列構成的新材料或新結構[3],屬于人工周期結構范疇。因此可以從軌道結構周期特征出發(fā),基于聲子晶體基本理論和基本方法研究彈性波在軌道結構中的傳播特性。

    目前針對人工周期結構的研究已取得較多成果,最初的研究主要關注實際工程中廣泛應用的結構,如周期性梁桿結構和板殼結構。文獻[4]研究周期性梁、板結構的振動特性,提出傳遞矩陣法、空間諧波法等計算方法,并首次給出無限周期梁板結構的彈性波帶隙特性。文獻[5]研究周期性Timoshenko梁中壓縮波、扭轉波、彎曲波的耦合振動特性。文獻[6]提出聲子晶體概念后,聲子晶體具有的帶隙特性得到廣泛關注。國內外學者針對帶隙產生機制、帶隙計算方法以及聲子晶體的應用等方面開展大量研究工作[2]。國防科技大學溫熙森教授團隊將聲子晶體帶隙原理引入典型機械結構減振降噪設計中,開展彈性波在周期性結構中傳播特性及控制的研究。文獻[7]采用平面波展開法研究彎曲波在局域共振梁結構中的傳播特性,實現寬頻帶隙設計。文獻[8]針對薄板結構分別通過理論和試驗測試驗證聲子晶體帶隙在控制低頻振動方面的能力。文獻[9-10]將聲子晶體理論引入充液管路系統(tǒng)設計中,利用帶隙特性實現充液管路的減振設計。

    對于鐵路軌道結構,也有大量學者從周期特征出發(fā)開展相應的研究工作。文獻[11]針對周期彈性點支撐軌道結構提出廣義狀態(tài)矩陣以求解其導納函數,并進行試驗驗證。文獻[12]將軌道結構劃分為周期尺寸為10 mm的單元,并將鋼軌截面通過有限單元描述,以得到50~5 000 Hz范圍內的軌道振動特性。文獻[13]針對連續(xù)彈性支撐梁模型、周期離散點支撐梁模型以及考慮鋼軌截面變形的連續(xù)支撐模型分別進行計算分析,結果表明橫向振動特性與鋼軌截面變形密切相關,采用考慮截面變形的軌道模型能夠有效預測橫向振動,而垂向振動則需考慮軌道結構周期性。文獻[14]分析周期性軌道結構的振動響應特性,并指出在較寬且明顯的衰減域內鋼軌噪聲輻射較小,且pinned-pinned振動可以在周期性支撐的鋼軌中自由傳播。文獻[15]采用多梁模型模擬鋼軌截面,運用格林公式和疊加原理計算有限長離散支撐鋼軌橫向振動,得到較準確的結果。文獻[16-17]利用2.5維有限元法計算移動荷載下周期性軌道結構的動力響應以及振動傳播常數隨頻率的變化曲線。文獻[18]在計算周期性軌道結構垂向導納以及振動衰減率特性時,在pinned-pinned共振頻率附近發(fā)現了振動傳播的阻帶(帶隙),但并未對這一特征進行深入分析。

    有砟鐵路軌道是典型的無限周期結構,目前,周期性軌道結構動力響應方面的研究已取得了大量成果,但有關彈性波在軌道結構中傳播特性的研究尚待深入。本文結合聲子晶體基本理論研究彈性波在普通有砟軌道結構中的傳播特性,并利用功率流方法分析周期性軌道結構中的能量傳播問題。

    1 有砟軌道結構彈性波帶隙特性

    以有砟軌道結構為例,軌道結構的彈性主要由扣件系統(tǒng)以及道床提供,如圖1(a)所示。此時將軌道結構簡化為由鋼軌、扣件、軌枕、道床組成的無限周期結構,并建立雙層彈性點支承軌道結構動力學模型,如圖1(b)所示。

    圖1 有砟軌道主要組成及其計算模型

    當彈性介質受到外力作用時,并非彈性介質中所有部分都立刻產生位移、應力和應變,而是外力作用開始之后,在作用處產生變形,使該點處質點產生振動并通過質點之間的相互作用在介質內由近及遠向外傳播,這種振動狀態(tài)在彈性介質中的傳播過程,稱為彈性波。對于均勻、連續(xù)、線性介質中彈性波的傳播,其波動方程可表示為[19]

    (1)

    式中:u為位移分量;ρ為介質密度;C11和C44為彈性模量元;p=x,y,z。

    研究彈性波在周期結構中傳播時,基于Bloch定理,對于任意一個給定的Bloch波數k,可以求出一系列對應的本征值和相應的本征矢,每一個本征值和本征矢都對應一個頻率和相應的運動狀態(tài)。波數k取遍整個不可約Brillouin區(qū),得到本征頻率隨波數k的變化曲線,稱為頻散曲線,也叫作能帶結構圖,其包含了彈性波傳播的相位關系及幅值衰減關系。

    根據Bloch定理,周期結構中的彈性波在每個周期的響應與相鄰周期相關[2]。

    um+1(x,t)=e-iklum(x,t)

    (2)

    軌道結構中鋼軌采用Timoshenko梁進行描述,因此鋼軌自由波動方程可表示為[19]

    (3)

    式中:G為鋼軌剪切模量;A為鋼軌橫截面積;κ為剪切因子;E為鋼軌彈性模量;I為鋼軌截面慣性矩;ρ為鋼軌密度;v為鋼軌垂向位移;φ為鋼軌彎曲轉角;u為鋼軌縱向位移。其自由波解為[19]

    (4)

    式中:V為垂向位移振幅;Vψ為轉角振幅;U為縱向位移振幅;kn為自由梁中波數。轉角位移比值關系為

    (5)

    此時,彎曲波中存在2對互為相反數的波數,縱波中存在1對互為相反數的波數,因此位移響應可寫為

    (6)

    鋼軌狀態(tài)向量可表示為

    ψ=[qf]=[vφuQMN]

    (7)

    式中:q為位移向量,包括垂向位移v、彎曲轉角φ、縱向位移u;f為力向量,包括剪力Q、彎矩M、軸力N。

    即可得到相鄰周期單元之間的傳遞矩陣為

    T=Tl·Ts

    (8)

    式中:Tl為彈性波傳播距離為l時的狀態(tài)向量傳遞矩陣;Ts為扣件兩端鋼軌狀態(tài)向量傳遞矩陣

    其中:kv、kl表示鋼軌下部支撐動態(tài)剛度。

    En為狀態(tài)向量因子

    結合Bloch定理可得周期性軌道結構彈性波傳播特征方程為

    |T-e-iklI|=0

    (9)

    式中:k為沿周期性軌道結構傳播的特征波(包括彎曲波和縱波)的波數,也稱Bloch波數,實部表示彈性波相位改變,虛部表示彈性波阻尼,即波的傳播衰減系數。求解式(9)可得波數k與頻率f之間的關系,即頻散特性。由于周期單元之間通過垂向位移、彎曲轉角、縱向位移3個自由度耦合,因此通過求解上述特征方程,可得周期性軌道中三種特征波的頻散特性曲線。

    1.1 軌道結構頻散特性

    軌道結構參數取值見表1。

    表1 有砟軌道結構參數

    圖2 周期軌道結構中彎曲波頻散曲線

    圖3 周期軌道結構中縱波頻散曲線

    周期性軌道結構中彎曲波頻散曲線如圖2所示,縱波頻散曲線如圖3所示,灰色區(qū)域為帶隙范圍。

    由圖2和圖3可知,在不考慮阻尼情況下,彎曲波和縱波在周期性軌道結構中的傳播存在明顯的通帶、禁帶特征。鋼軌中彎曲波由兩類特征波組成[20]:第一類彎曲波在全頻段內波數虛部均不為0,即全頻段均為帶隙范圍,因此第一類彎曲波為近場波;第二類彎曲波則在0~1 500 Hz范圍內存在3個帶隙,分別為0~129 Hz、182~262 Hz、1 080~1 125 Hz,因此鋼軌中彎曲波帶隙特性由第二類彎曲波決定??v波在0~300 Hz范圍內存在兩個帶隙,分別為0~24 Hz、44~80 Hz。在這些頻率范圍內,軌道結構中彈性波波數虛部均不為0,因此在傳播過程中快速衰減。

    1.2 軌道結構參數對帶隙特性的影響

    彈性波帶隙特性是周期性軌道結構的固有屬性,軌道結構參數直接影響其帶隙特征。扣件、道床垂向剛度以及扣件間距對軌道結構彎曲波帶隙的影響如圖4所示,不同顏色代表帶隙對彈性波傳播的衰減能力。由圖4可知,隨著扣件剛度增加,帶隙寬度逐漸增大;隨著道床剛度的增加,第一階帶隙寬度逐漸增大,第二階帶隙寬度逐漸減小,而第三階帶隙位置與寬度均保持不變;隨著扣件間距增加,第一階和第二階帶隙寬度略有減小,第三階帶隙逐漸向低頻移動但帶隙寬度保持不變。因此可以通過改變扣件剛度、道床剛度以及扣件間距等措施,控制軌道結構中振動的傳播。

    圖4 軌道結構參數對彎曲波帶隙的影響

    周期性軌道結構頻散關系表明彈性波在軌道結構中的傳播具有帶隙特性,帶隙范圍內彈性波的傳播受到抑制,彈性波在通帶范圍內可自由傳播。彈性波在軌道結構的傳播過程中必然會向外界傳播振動或輻射噪聲,對于鋼軌振動,可以設計軌道結構帶隙,將其限制在局部位置,抑制振動的傳播,進一步通過阻尼耗能等措施實現振動的控制;由于鋼軌振動輻射的噪聲與振動衰減特性密切相關[18],利用帶隙特性可實現振動的衰減,從而降低鋼軌噪聲。另一方面,由于彈性波在軌道結構中的傳播過程實質上是振動能量的傳遞過程,因此有必要從能量的角度出發(fā),進一步深入分析周期軌道結構中彈性波的傳播。

    2 周期性軌道結構動力響應

    單一的力或位移傳遞率不能有效地評價結構中的能量傳遞,功率流方法綜合了力和速度響應的大小及相位關系,給出了振動傳輸的一個絕對度量[21]。通過對功率流分析可以更加深入地研究彈性波能量在軌道結構中的傳播,揭示振動傳輸機制。

    功率流分析前需要得到軌道結構的動力響應結果,因此首先計算周期性軌道結構動力響應。由于軌道結構具有明顯的周期性,彈性波在周期性軌道結構中的傳播具有帶隙特征,且不同成分的彈性波具有不同的傳播性質。因此根據周期結構中彈性波的傳播特性,利用彈性波疊加原理,即可求得周期性軌道結構在荷載作用下的動力響應。無限長自由梁在集中力荷載作用下的系統(tǒng)動力響應可表示為[22]

    (10)

    式中:F為外力荷載;xr為響應位置;x0為荷載激勵位置。

    圖5 無限長周期支撐梁

    在鋼軌上方作用外力F0時,如圖5所示,將周期支撐視為支反力作用于無限長鋼軌上,根據Bloch定理,支反力滿足周期條件[5]

    (11)

    式中:Fsr為激勵右側第s跨處支反力;Fsl為激勵左側第s跨處支反力;s為正整數;u為周期結構中的傳播常數。傳播常數與Bloch波數之間的關系為[5]

    u=-ikl

    (12)

    將式(11)代入式(10),通過計算在所有支反力和外力作用下無限長梁的位移響應,利用彈性波疊加原理即可得到軌道結構某一位置處的位移。因此,在受載跨內xr處位移響應可表示為[23]

    (13)

    式中:0≤|xr|≤l且將坐標原點設置為0l扣件所在位置。

    式(13)中的兩個無限級數由軌道結構中自由跨支反力組成,由于滿足周期條件,可簡化為

    (14)

    對于周期性軌道結構,將鋼軌簡化為Timoshenko梁,鋼軌彎曲波沿正負方向存在兩種不同的傳播常數±u1、±u2,因此每個支反力需要分為兩個部分,分別對應不同成分的彈性波。在扣件處支反力和鋼軌位移滿足以下平衡條件

    -Fs/kv=v(xs)

    (15)

    式中:Fs為第s個支撐點作用力;kv為軌下等效動剛度;v為該支撐點處鋼軌位移。兩類彎曲波成分引起的扣件支反力可表示為

    (16)

    第s個支撐點處總反力為

    (17)

    將式(17)代入式(14)可得受載跨內任意位置的鋼軌位移響應

    (18)

    此時響應函數中存在6個未知變量:F0l、F0r、F11l、F12l、F11r、F12r,需要6個平衡條件才能求解。選擇支撐位置為2r、1r、0r、0l、1l、2l處力和位移平衡條件,即可求得6個未知數,從而得到軌道結構動力響應。以0r處平衡條件為例:將xr=l帶入式(18)并結合力和位移平衡條件(式(15)),可得到0r處平衡方程

    (19)

    同理可寫出其他位置處平衡方程,求解得到F0l、F0r、F11l、F12l、F11r、F12r并代入式(18)即可求得無限長周期軌道結構在簡諧荷載作用下的動力響應。

    以表1中有砟軌道結構參數為例,得到無限長周期軌道結構在簡諧荷載(幅值為1 N)作用下鋼軌原點的導納函數,如圖6所示。由于振動的本質為彈性波在結構中的傳播,因此周期性軌道結構固有頻率與其彈性波帶隙頻率一一對應。同時建立有限長軌道結構模型(由160個元胞組成),采用有限元方法計算得到扣件跨中位置處原點的導納函數,如圖7所示??梢钥闯觯瑑煞N方法基本能夠得到一致的結果。但對于有限長軌道結構,當頻率處于通帶范圍內時,由于此時彈性波傳播距離較遠,邊界處的反射將使計算結果不準確;帶隙范圍內彈性波在較短距離內衰減,采用有限長軌道結構即可準確描述。

    圖6 無限長周期性軌道結構鋼軌導納函數

    圖7 有限、無限長軌道結構鋼軌導納函數對比

    3 周期性軌道結構功率流分析

    若某點所受簡諧力為Feiωt,該點對應的速度響應為Veiωt,功率流計算公式為[21]

    (20)

    式中:下標i為節(jié)點編號。

    將鋼軌考慮為Timoshenko梁時,鋼軌中彎曲波功率流可分為剪力和彎矩分別攜帶的功率流兩部分。當在鋼軌跨中位置以及1/4跨位置施加單位簡諧垂向激勵時,鋼軌剪力和彎矩響應如圖8、圖9所示。在跨中激勵時,鋼軌剪力保持不變,始終為輸入力幅值的一半,而彎矩則在軌道結構固有頻率處出現極值。在1/4跨位置激勵時,鋼軌剪力和彎矩均在固有頻率處出現極值。

    圖8 跨中位置施加激勵時產生的剪力和彎矩

    圖9 1/4跨位置施加激勵時產生的剪力和彎矩

    在鋼軌跨中位置施加單位簡諧激勵時,輸入軌道結構功率流以及鋼軌傳播功率流如圖10所示。由計算結果可知,輸入軌道結構中的功率流也存在帶隙特性,在禁帶范圍內不能向軌道中輸入功率流,且在通帶范圍內,輸入系統(tǒng)的功率流一半沿正向傳播,一半沿負方向傳播。當考慮軌道結構阻尼后,通帶內輸入的功率流變化不明顯,而帶隙范圍內可向軌道結構輸入功率流。

    圖10 輸入軌道結構功率流(灰色區(qū)域為帶隙范圍)

    在不同位置處激勵時,軌道結構輸入功率流隨頻率的變化曲線如圖11所示。由圖11可以看出,在第一帶隙、第一通帶和第二帶隙范圍內,不同位置的激勵輸入軌道結構的功率流大小基本相同;當頻率處于第二階通帶頻率范圍內時(258~1 013 Hz),靠近跨中位置輸入功率流較大;但在第三階通帶頻率范圍內則靠近扣件位置輸入功率流較大。這是由不同頻率下軌道結構的振動模態(tài)決定的,當作用在跨中位置時,可激勵產生pinned-pinned共振,此時輸入鋼軌和傳播的功率流達到最大。

    圖11 不同位置激勵時輸入軌道結構功率流(考慮阻尼情況)

    當在扣件上方激勵時,輸入功率流、傳播功率流以及激勵位置處軌枕功率流如圖12所示。由圖12可知軌枕振動能量主要集中在300 Hz以內。而在第二階帶隙范圍內,軌枕功率流接近鋼軌功率流,此時軌枕發(fā)生共振,軌枕作為局域振子能夠吸收大量振動能量。彎曲波功率流在周期性軌道結構中的傳播系數隨頻率的變化如圖13所示,其中傳輸系數的參考值為激勵位置處的鋼軌功率流。從圖13可以看出,彎曲波功率流在軌道結構中的傳播呈明顯的帶隙特征,在帶隙頻率范圍內彎曲波能量的傳播得到明顯的抑制。

    圖12 鋼軌、軌枕功率流(扣件上方激勵)

    圖13 彎曲波功率流沿鋼軌縱向傳播情況

    在不同的軌道結構參數下,在鋼軌跨中激勵時輸入軌道結構功率流如圖14所示。與軌道結構參數對帶隙的調控規(guī)律類似,當增加扣件剛度時,第一帶隙范圍內功率流降低,且峰值向高頻移動,在第二帶隙范圍內整體向高頻移動,增加扣件剛度導致第三帶隙拓寬;增加道床剛度則對第一帶隙范圍改變明顯,對第三帶隙沒有影響;改變扣件間距對第一和第二帶隙峰值略有影響,但主要影響第三帶隙的位置。因此,通過設計合理的軌道結構參數可對彈性波傳播進行調控,達到振動噪聲控制的目的。

    圖14 軌道結構參數對功率流影響

    4 結論

    本文從軌道結構周期特性出發(fā),以有砟軌道結構為例,結合聲子晶體基本理論,開展彈性波在周期性軌道結構中傳播特性的研究。基于Bloch定理并結合傳遞矩陣法,得到周期性軌道結構頻散曲線及其帶隙特征;基于彈性波疊加原理并結合周期特征求解無限長周期軌道結構在簡諧荷載作用下的動力響應,利用功率流方法分析振動能量在軌道結構中的傳播特性。得到以下結論:

    (1)周期性軌道結構具有明顯的帶隙特性,在0~1 500 Hz范圍內,鋼軌彎曲波存在3個帶隙,縱波存在2個帶隙,帶隙范圍以內彈性波的傳播受到抑制。

    (2) 對于有限長軌道結構,當頻率處于通帶范圍內時,由于彈性波傳播較遠,在邊界處反射將造成計算結果不準確;帶隙范圍內彈性波在較短距離內衰減,采用有限長軌道結構即可準確描述。

    (3) 周期性軌道結構中功率流傳播同樣存在通帶、禁帶特性。在禁帶頻率范圍內,振源不能向軌道結構輸入功率流;在通帶范圍內能夠輸入功率流且輸入的功率流中一半正向傳播,一半負向傳播。

    (4) 軌道結構參數對彈性波傳播特性產生明顯的影響,因此可以通過設計合理的軌道結構來調控彈性波在軌道結構中的傳播,為軌道結構減振降噪設計提供參考。

    猜你喜歡
    周期性扣件鋼軌
    科隆蛋扣件力學特性分析
    數列中的周期性和模周期性
    中等數學(2019年5期)2019-08-30 03:52:20
    一種改進的LBP特征實現鐵路扣件識別
    一類整數遞推數列的周期性
    中等數學(2018年12期)2018-02-16 07:48:40
    鋼軌焊后雙中頻感應加熱工藝研究
    非均布荷載下扣件式腳手架的隨機缺陷穩(wěn)定分析
    基于擴頻碼周期性的單通道直擴通信半盲分離抗干擾算法
    高速鐵路鋼軌疲勞過程的超聲非線性系數表征
    國內外高速鐵路鋼軌性能對比研究
    諧振式浮軌扣件的安裝工藝
    国产亚洲av高清不卡| 人人澡人人妻人| 国产在视频线精品| 精品亚洲成a人片在线观看| 国产精品二区激情视频| 国产色视频综合| 日韩中文字幕视频在线看片| 不卡av一区二区三区| 99久久人妻综合| 蜜桃在线观看..| 肉色欧美久久久久久久蜜桃| 亚洲av男天堂| 高清欧美精品videossex| 狠狠精品人妻久久久久久综合| 日本欧美视频一区| 香蕉丝袜av| 国产国语露脸激情在线看| 午夜激情久久久久久久| 国产精品免费大片| 两性午夜刺激爽爽歪歪视频在线观看 | 成人三级做爰电影| 亚洲va日本ⅴa欧美va伊人久久 | 日本一区二区免费在线视频| 黄色毛片三级朝国网站| 国产高清视频在线播放一区 | 国产成人精品久久二区二区免费| 18禁观看日本| 女性生殖器流出的白浆| 一进一出抽搐动态| 亚洲国产欧美网| 在线亚洲精品国产二区图片欧美| 日韩,欧美,国产一区二区三区| 又大又爽又粗| 国产人伦9x9x在线观看| 亚洲欧美日韩另类电影网站| 日韩欧美一区二区三区在线观看 | 欧美精品高潮呻吟av久久| videosex国产| 精品国产一区二区久久| 王馨瑶露胸无遮挡在线观看| 精品亚洲乱码少妇综合久久| 精品国产一区二区三区久久久樱花| 美女主播在线视频| 亚洲精品久久午夜乱码| 亚洲国产成人一精品久久久| 亚洲中文av在线| 制服诱惑二区| 国精品久久久久久国模美| 另类精品久久| 日本黄色日本黄色录像| 中文字幕人妻熟女乱码| 午夜福利视频精品| 97在线人人人人妻| 日本a在线网址| 欧美黄色淫秽网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品久久久久久婷婷小说| 午夜久久久在线观看| 亚洲精品国产av成人精品| 国产深夜福利视频在线观看| 国产成人精品在线电影| 国产精品亚洲av一区麻豆| 亚洲av片天天在线观看| 精品亚洲成国产av| 亚洲精品中文字幕一二三四区 | 欧美 亚洲 国产 日韩一| 午夜免费成人在线视频| 黄频高清免费视频| 熟女少妇亚洲综合色aaa.| 又黄又粗又硬又大视频| 自线自在国产av| 久久国产精品男人的天堂亚洲| 女人久久www免费人成看片| 亚洲性夜色夜夜综合| 日本wwww免费看| 中文字幕av电影在线播放| 每晚都被弄得嗷嗷叫到高潮| 午夜福利免费观看在线| av又黄又爽大尺度在线免费看| 亚洲va日本ⅴa欧美va伊人久久 | 99久久精品国产亚洲精品| 一区二区三区乱码不卡18| 久久久国产欧美日韩av| 欧美精品一区二区大全| 国产精品一区二区免费欧美 | 国产成人免费观看mmmm| 老熟妇乱子伦视频在线观看 | 爱豆传媒免费全集在线观看| 一级a爱视频在线免费观看| av视频免费观看在线观看| 免费看十八禁软件| 亚洲 欧美一区二区三区| 久久久久久久大尺度免费视频| 亚洲成av片中文字幕在线观看| 亚洲天堂av无毛| 久久国产亚洲av麻豆专区| 成人国产一区最新在线观看| 亚洲国产精品一区二区三区在线| 两个人看的免费小视频| 久久人人97超碰香蕉20202| 色94色欧美一区二区| e午夜精品久久久久久久| 我要看黄色一级片免费的| 国产xxxxx性猛交| www.自偷自拍.com| 国产有黄有色有爽视频| 一本色道久久久久久精品综合| 色婷婷久久久亚洲欧美| 欧美黑人欧美精品刺激| 老司机午夜十八禁免费视频| 在线 av 中文字幕| 捣出白浆h1v1| 一本色道久久久久久精品综合| 在线 av 中文字幕| 欧美午夜高清在线| 国产国语露脸激情在线看| 搡老熟女国产l中国老女人| 无遮挡黄片免费观看| 亚洲欧美一区二区三区黑人| 两个人看的免费小视频| 亚洲专区中文字幕在线| 纵有疾风起免费观看全集完整版| 亚洲av电影在线观看一区二区三区| 999久久久精品免费观看国产| 亚洲精品国产精品久久久不卡| 日本wwww免费看| 我要看黄色一级片免费的| 久久亚洲精品不卡| 欧美日韩福利视频一区二区| 国产老妇伦熟女老妇高清| 少妇被粗大的猛进出69影院| 99热国产这里只有精品6| 国产免费一区二区三区四区乱码| 两个人免费观看高清视频| 母亲3免费完整高清在线观看| 国产人伦9x9x在线观看| 精品亚洲乱码少妇综合久久| 1024视频免费在线观看| 午夜福利视频精品| 他把我摸到了高潮在线观看 | 久久 成人 亚洲| 亚洲少妇的诱惑av| 久久精品成人免费网站| 日本av手机在线免费观看| 精品亚洲乱码少妇综合久久| 最新的欧美精品一区二区| 精品少妇久久久久久888优播| 国产深夜福利视频在线观看| 亚洲七黄色美女视频| 久久久久精品国产欧美久久久 | 久久久久国内视频| 天天躁日日躁夜夜躁夜夜| 男女无遮挡免费网站观看| 嫁个100分男人电影在线观看| 人人妻人人添人人爽欧美一区卜| 91成人精品电影| 亚洲av欧美aⅴ国产| 国产成人av激情在线播放| 国产深夜福利视频在线观看| 亚洲精华国产精华精| 两个人看的免费小视频| 啦啦啦免费观看视频1| 亚洲av国产av综合av卡| 亚洲av片天天在线观看| 久久中文看片网| 一区二区av电影网| 国产精品99久久99久久久不卡| 中文字幕人妻丝袜一区二区| 日韩中文字幕视频在线看片| 久久久精品区二区三区| 性色av乱码一区二区三区2| 欧美国产精品一级二级三级| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品久久二区二区91| 两个人看的免费小视频| 日韩视频在线欧美| 黄色片一级片一级黄色片| 老司机亚洲免费影院| 亚洲成人国产一区在线观看| 一本—道久久a久久精品蜜桃钙片| 久久国产精品影院| 国内毛片毛片毛片毛片毛片| 黄色a级毛片大全视频| 日本五十路高清| 电影成人av| 亚洲精品久久午夜乱码| 午夜免费观看性视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲一区中文字幕在线| 国产一区二区三区综合在线观看| 亚洲av男天堂| 女人被躁到高潮嗷嗷叫费观| 欧美xxⅹ黑人| av在线播放精品| 国产日韩欧美在线精品| 国产欧美日韩综合在线一区二区| 久久久久国产精品人妻一区二区| 俄罗斯特黄特色一大片| 色视频在线一区二区三区| 日韩欧美国产一区二区入口| 香蕉国产在线看| 一区二区三区乱码不卡18| 亚洲精品一二三| 精品国产一区二区久久| 亚洲精品国产区一区二| 老鸭窝网址在线观看| 精品国内亚洲2022精品成人 | 国产精品免费视频内射| 黄频高清免费视频| 一级毛片女人18水好多| 久久久精品免费免费高清| 18禁黄网站禁片午夜丰满| 老司机深夜福利视频在线观看 | 丰满少妇做爰视频| 亚洲成av片中文字幕在线观看| 爱豆传媒免费全集在线观看| 欧美+亚洲+日韩+国产| 一级片免费观看大全| 久久 成人 亚洲| 人妻 亚洲 视频| 国产一区二区在线观看av| 久久人人爽人人片av| 一个人免费在线观看的高清视频 | 在线观看一区二区三区激情| 亚洲色图 男人天堂 中文字幕| 久久精品亚洲熟妇少妇任你| 女性生殖器流出的白浆| 我的亚洲天堂| 在线永久观看黄色视频| av又黄又爽大尺度在线免费看| 黄色 视频免费看| 男女免费视频国产| 黑人操中国人逼视频| 蜜桃在线观看..| 日韩大片免费观看网站| 国产深夜福利视频在线观看| 波多野结衣一区麻豆| 成年人午夜在线观看视频| 啦啦啦 在线观看视频| 国产成人啪精品午夜网站| 99久久人妻综合| 汤姆久久久久久久影院中文字幕| 成年人午夜在线观看视频| 国产三级黄色录像| 国产精品一区二区免费欧美 | 日本av手机在线免费观看| 日韩 欧美 亚洲 中文字幕| 国产极品粉嫩免费观看在线| 久久久久久免费高清国产稀缺| 日本a在线网址| 亚洲精品av麻豆狂野| 久热这里只有精品99| 99久久精品国产亚洲精品| 精品国产超薄肉色丝袜足j| 51午夜福利影视在线观看| 啦啦啦视频在线资源免费观看| 日韩欧美一区二区三区在线观看 | 女人爽到高潮嗷嗷叫在线视频| 黄片小视频在线播放| 可以免费在线观看a视频的电影网站| 亚洲国产欧美日韩在线播放| 欧美亚洲 丝袜 人妻 在线| 99九九在线精品视频| 99re6热这里在线精品视频| 午夜久久久在线观看| 97在线人人人人妻| 免费观看a级毛片全部| 日韩 亚洲 欧美在线| 叶爱在线成人免费视频播放| 久久国产亚洲av麻豆专区| 国产91精品成人一区二区三区 | 中文字幕精品免费在线观看视频| 亚洲国产日韩一区二区| 人人妻人人添人人爽欧美一区卜| a级片在线免费高清观看视频| 爱豆传媒免费全集在线观看| 男女午夜视频在线观看| 一区二区av电影网| 狠狠精品人妻久久久久久综合| 丝袜人妻中文字幕| 国产精品1区2区在线观看. | 日韩欧美一区二区三区在线观看 | 精品福利永久在线观看| 老司机深夜福利视频在线观看 | 一边摸一边抽搐一进一出视频| 欧美日本中文国产一区发布| 两个人看的免费小视频| 亚洲成人国产一区在线观看| 久久ye,这里只有精品| 在线观看一区二区三区激情| 国产成人欧美在线观看 | 亚洲精品在线美女| 麻豆av在线久日| 欧美人与性动交α欧美软件| 丝袜美腿诱惑在线| 亚洲av片天天在线观看| 捣出白浆h1v1| 免费在线观看日本一区| 在线精品无人区一区二区三| 国产成人av激情在线播放| 少妇被粗大的猛进出69影院| www.熟女人妻精品国产| 国产又爽黄色视频| 欧美一级毛片孕妇| 狠狠精品人妻久久久久久综合| 青草久久国产| 另类精品久久| 国产成人影院久久av| 亚洲国产精品成人久久小说| cao死你这个sao货| 丰满少妇做爰视频| 叶爱在线成人免费视频播放| av电影中文网址| 日韩熟女老妇一区二区性免费视频| 9色porny在线观看| 人妻一区二区av| 一级毛片电影观看| 国产欧美日韩一区二区三 | 国产日韩欧美亚洲二区| 极品人妻少妇av视频| 在线亚洲精品国产二区图片欧美| av福利片在线| 91字幕亚洲| 欧美人与性动交α欧美精品济南到| 麻豆乱淫一区二区| 嫁个100分男人电影在线观看| 国产精品1区2区在线观看. | 亚洲第一青青草原| 精品第一国产精品| 国产精品二区激情视频| 精品少妇久久久久久888优播| 免费在线观看完整版高清| 丰满人妻熟妇乱又伦精品不卡| 欧美精品高潮呻吟av久久| 两人在一起打扑克的视频| 天天躁夜夜躁狠狠躁躁| 亚洲欧美成人综合另类久久久| 美女视频免费永久观看网站| 69精品国产乱码久久久| 欧美 日韩 精品 国产| 亚洲美女黄色视频免费看| 97精品久久久久久久久久精品| cao死你这个sao货| 日韩 亚洲 欧美在线| 天天操日日干夜夜撸| 免费在线观看视频国产中文字幕亚洲 | 搡老乐熟女国产| 一区二区日韩欧美中文字幕| 亚洲精品国产色婷婷电影| 在线观看免费视频网站a站| 777米奇影视久久| 无遮挡黄片免费观看| 国产区一区二久久| 天天影视国产精品| 高清av免费在线| 青春草视频在线免费观看| 国产男人的电影天堂91| 欧美97在线视频| 色综合欧美亚洲国产小说| 国产又爽黄色视频| 国产亚洲午夜精品一区二区久久| 精品人妻在线不人妻| 大片电影免费在线观看免费| av欧美777| 国产精品免费大片| 国产免费av片在线观看野外av| 熟女少妇亚洲综合色aaa.| 国产淫语在线视频| 一级毛片电影观看| 黄片小视频在线播放| 国产日韩一区二区三区精品不卡| 欧美黄色片欧美黄色片| 啦啦啦中文免费视频观看日本| 日本黄色日本黄色录像| 91精品伊人久久大香线蕉| 久久精品人人爽人人爽视色| 丰满少妇做爰视频| 久久精品国产亚洲av高清一级| 在线十欧美十亚洲十日本专区| 亚洲av日韩在线播放| 成年人午夜在线观看视频| 永久免费av网站大全| 久久久久网色| 1024视频免费在线观看| 另类精品久久| 日本一区二区免费在线视频| 男女午夜视频在线观看| 婷婷丁香在线五月| 两个人免费观看高清视频| 高清在线国产一区| 91精品三级在线观看| av在线播放精品| 国产精品.久久久| 精品高清国产在线一区| 国产有黄有色有爽视频| 亚洲一区中文字幕在线| 国产黄色免费在线视频| 在线观看www视频免费| 中文欧美无线码| 亚洲精品在线美女| 国产国语露脸激情在线看| 久久精品国产a三级三级三级| 亚洲伊人久久精品综合| 国产xxxxx性猛交| 国产激情久久老熟女| 久久精品国产亚洲av香蕉五月 | 91麻豆av在线| 在线永久观看黄色视频| 日本av手机在线免费观看| 欧美激情极品国产一区二区三区| 咕卡用的链子| 高清黄色对白视频在线免费看| a 毛片基地| 丝袜人妻中文字幕| 欧美日韩亚洲高清精品| 菩萨蛮人人尽说江南好唐韦庄| 日本五十路高清| 爱豆传媒免费全集在线观看| 人妻 亚洲 视频| 国产日韩一区二区三区精品不卡| 岛国毛片在线播放| 欧美成狂野欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 99久久综合免费| 最近中文字幕2019免费版| 亚洲专区字幕在线| 人成视频在线观看免费观看| 99热全是精品| 免费在线观看黄色视频的| 精品亚洲乱码少妇综合久久| av网站在线播放免费| 日韩欧美一区视频在线观看| 色视频在线一区二区三区| 五月开心婷婷网| 国产亚洲精品一区二区www | 免费女性裸体啪啪无遮挡网站| 男人爽女人下面视频在线观看| 久久免费观看电影| 国产精品久久久久久精品古装| 精品国产乱码久久久久久小说| 亚洲伊人色综图| av又黄又爽大尺度在线免费看| 人妻 亚洲 视频| 国产精品一区二区免费欧美 | 国产老妇伦熟女老妇高清| 亚洲精品自拍成人| 老汉色∧v一级毛片| 精品国产一区二区久久| 精品一区二区三卡| 黑人欧美特级aaaaaa片| 女性生殖器流出的白浆| 色婷婷久久久亚洲欧美| 亚洲五月婷婷丁香| 日韩制服丝袜自拍偷拍| 老司机福利观看| 亚洲成人国产一区在线观看| 美女福利国产在线| 亚洲一码二码三码区别大吗| 黄网站色视频无遮挡免费观看| 欧美日韩一级在线毛片| 久久久国产欧美日韩av| 亚洲精品国产色婷婷电影| 欧美日本中文国产一区发布| 无限看片的www在线观看| 中亚洲国语对白在线视频| 老汉色∧v一级毛片| 日本vs欧美在线观看视频| 免费少妇av软件| 成人影院久久| 成人免费观看视频高清| 不卡av一区二区三区| 这个男人来自地球电影免费观看| 91字幕亚洲| 国产男人的电影天堂91| 欧美黑人欧美精品刺激| 一区二区三区乱码不卡18| 一级片免费观看大全| 999精品在线视频| 性少妇av在线| 久久亚洲精品不卡| 亚洲色图综合在线观看| 成在线人永久免费视频| 久热爱精品视频在线9| 久久香蕉激情| 巨乳人妻的诱惑在线观看| 国产日韩一区二区三区精品不卡| 国产一区二区三区在线臀色熟女 | 免费观看a级毛片全部| 涩涩av久久男人的天堂| www.熟女人妻精品国产| 欧美日韩中文字幕国产精品一区二区三区 | 日日夜夜操网爽| 免费女性裸体啪啪无遮挡网站| 亚洲精品第二区| 日本一区二区免费在线视频| 国产淫语在线视频| 精品国内亚洲2022精品成人 | 久久精品熟女亚洲av麻豆精品| 91老司机精品| 肉色欧美久久久久久久蜜桃| 丰满饥渴人妻一区二区三| 亚洲国产av新网站| 在线观看免费视频网站a站| 日本av手机在线免费观看| 欧美激情高清一区二区三区| 国产精品一二三区在线看| 少妇人妻久久综合中文| 国产又色又爽无遮挡免| 王馨瑶露胸无遮挡在线观看| 亚洲精品自拍成人| 乱人伦中国视频| 夫妻午夜视频| videosex国产| 久久久久国产精品人妻一区二区| 人人妻人人澡人人看| 亚洲国产看品久久| 一区二区三区四区激情视频| 男女免费视频国产| 免费在线观看影片大全网站| 中文字幕色久视频| 男女床上黄色一级片免费看| 老司机深夜福利视频在线观看 | 国产亚洲av片在线观看秒播厂| 一级毛片电影观看| 亚洲av电影在线进入| 国产精品99久久99久久久不卡| 国产极品粉嫩免费观看在线| 欧美黑人精品巨大| 在线观看免费高清a一片| 亚洲一码二码三码区别大吗| 久久久久久久大尺度免费视频| 一级片'在线观看视频| 日韩,欧美,国产一区二区三区| 老司机影院成人| 后天国语完整版免费观看| 在线观看人妻少妇| 91老司机精品| 亚洲国产精品一区二区三区在线| 成人av一区二区三区在线看 | 捣出白浆h1v1| a级毛片黄视频| 精品国产乱码久久久久久小说| 一区福利在线观看| 国产免费视频播放在线视频| 中亚洲国语对白在线视频| 丝瓜视频免费看黄片| 精品一区在线观看国产| 国产av精品麻豆| 国产亚洲av高清不卡| 夜夜夜夜夜久久久久| 亚洲第一青青草原| 我要看黄色一级片免费的| 欧美日韩视频精品一区| 久久狼人影院| 男男h啪啪无遮挡| 视频在线观看一区二区三区| 精品亚洲成a人片在线观看| 黑人欧美特级aaaaaa片| 亚洲精品中文字幕一二三四区 | 日韩熟女老妇一区二区性免费视频| 日韩 亚洲 欧美在线| 菩萨蛮人人尽说江南好唐韦庄| 91麻豆精品激情在线观看国产 | 十八禁人妻一区二区| 午夜福利视频在线观看免费| 久久久久久久久免费视频了| 亚洲欧美一区二区三区黑人| 国产欧美日韩一区二区三区在线| 国产精品免费大片| 亚洲av电影在线进入| 建设人人有责人人尽责人人享有的| 日韩大片免费观看网站| 捣出白浆h1v1| 亚洲欧美精品综合一区二区三区| 99久久精品国产亚洲精品| 久久天堂一区二区三区四区| 高清欧美精品videossex| 国产熟女午夜一区二区三区| 国产真人三级小视频在线观看| 国产精品 欧美亚洲| 99国产精品免费福利视频| 一本大道久久a久久精品| 免费观看a级毛片全部| 欧美精品高潮呻吟av久久| 男女高潮啪啪啪动态图| 日韩三级视频一区二区三区| 欧美激情久久久久久爽电影 | 美女福利国产在线| a在线观看视频网站| 亚洲国产毛片av蜜桃av| 窝窝影院91人妻| 国产黄色免费在线视频| 国产一区有黄有色的免费视频| 一级毛片精品| 免费看十八禁软件| 国产精品麻豆人妻色哟哟久久| 丁香六月欧美| 亚洲精品久久久久久婷婷小说| 天天影视国产精品| 亚洲中文日韩欧美视频| 亚洲精品久久久久久婷婷小说| 午夜影院在线不卡| 中文字幕av电影在线播放| 亚洲精品久久久久久婷婷小说| 巨乳人妻的诱惑在线观看| 啦啦啦视频在线资源免费观看| 丝袜人妻中文字幕| 久久国产精品男人的天堂亚洲| 丝袜美足系列| 电影成人av| 悠悠久久av| 欧美黄色淫秽网站| 1024香蕉在线观看| 搡老熟女国产l中国老女人| 久久久国产成人免费|