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

    螺栓連接鼓筒轉(zhuǎn)子結(jié)構(gòu)動力學(xué)特性分析

    2019-08-10 03:48:22孟春曉張文勝馬輝秦朝燁樊富友
    振動工程學(xué)報 2019年3期

    孟春曉 張文勝 馬輝 秦朝燁 樊富友

    摘要: 以螺栓連接鼓筒轉(zhuǎn)子為研究對象,基于有限元軟件,首先建立了減縮梁殼彈簧混合單元模型;其次通過仿真與實驗的固有特性對比驗證了減縮模型的準確性;最后分析了螺栓個數(shù)、螺栓松動、系統(tǒng)轉(zhuǎn)速對螺栓連接鼓筒轉(zhuǎn)子結(jié)構(gòu)結(jié)合面的時變剛度和系統(tǒng)的響應(yīng)特性的影響規(guī)律。研究結(jié)果表明:螺栓個數(shù)、螺栓松動與系統(tǒng)轉(zhuǎn)速對連接結(jié)構(gòu)的時變剛度和響應(yīng)特性均有一定的影響,隨著螺栓個數(shù)的增多,結(jié)合面連接剛度不斷增大,系統(tǒng)的非線性特性不斷減弱;螺栓松動導(dǎo)致連接剛度發(fā)生較大波動,并且隨著松動個數(shù)的增多,系統(tǒng)非線性不斷增強;隨著系統(tǒng)轉(zhuǎn)速的增大,螺栓連接時變剛度增大,但轉(zhuǎn)速越高其剛度波動越大,導(dǎo)致系統(tǒng)非線性增強。

    關(guān)鍵詞: 轉(zhuǎn)子動力學(xué); 螺栓連接; 鼓筒轉(zhuǎn)子; 時變剛度; 螺栓松動

    中圖分類號:O347.6; TH113.1??文獻標志碼: A??文章編號: 1004-4523(2019)03.0517.09

    DOI:10.16385/j.cnki.issn.10044523.2019.03.017

    引?言

    螺栓連接的盤鼓式轉(zhuǎn)子結(jié)構(gòu)廣泛應(yīng)用于大型航空發(fā)動機[1],其結(jié)構(gòu)可靠性直接影響了航空發(fā)動機轉(zhuǎn)軸的轉(zhuǎn)動特性與穩(wěn)定性[2]。并且由于連接結(jié)構(gòu)的存在,轉(zhuǎn)子系統(tǒng)的動力學(xué)特性變得越來越復(fù)雜,而螺栓的松動又會導(dǎo)致結(jié)構(gòu)的響應(yīng)特性發(fā)生較大變化。

    針對螺栓連接結(jié)構(gòu),文獻[34]提出了多種建模方法,其中基于接觸的預(yù)應(yīng)力模態(tài)法,可準確模擬螺栓的接觸非線性,具有較高的精度。趙丹等[5]基于該方法,以盤盤螺栓連接轉(zhuǎn)子結(jié)構(gòu)為研究對象,采用預(yù)緊力單元法模擬螺栓預(yù)緊力,研究了其對系統(tǒng)固有特性的影響。針對渦扇發(fā)動機中的螺栓連接,翟學(xué)等[6]分別采用完全耦合建模方法與彈簧參數(shù)化建模方法模擬螺栓連接,并通過與實驗?zāi)B(tài)結(jié)構(gòu)對比,指出參數(shù)化建模方法在保證建模精度的同時,提高了求解效率。馬雙超[7]基于實驗結(jié)果,對連接界面的材料參數(shù)進行了辨識,通過改變材料參數(shù)來調(diào)整螺栓連接剛度,進而對其進行模態(tài)分析,指出隨著預(yù)緊力的增大,系統(tǒng)固有頻率增大。

    為解決螺栓連接接觸非線性所帶來的求解效率較低問題,可采用模態(tài)綜合法對其進行降維處理。目前,模態(tài)綜合法的基本理論從最初Hurty提出的固定界面法,已經(jīng)發(fā)展為包含固定界面法、自由界面法、混合界面法的多種模型減縮方法[8],并且被應(yīng)用到了轉(zhuǎn)子結(jié)構(gòu)中。文獻[9]將Guyan減縮的基本原理與轉(zhuǎn)子動力學(xué)相結(jié)合,對某渦扇發(fā)動機轉(zhuǎn)子進行了降維處理,并驗證了其可行性。文獻[10]提出了適用于具有轉(zhuǎn)速效應(yīng)的單轉(zhuǎn)子或多轉(zhuǎn)子的旋轉(zhuǎn)子結(jié)構(gòu)減縮方法。孫傳宗等[11]應(yīng)用CraigBampton模態(tài)綜合法,針對復(fù)雜雙轉(zhuǎn)子結(jié)構(gòu)進行了模型減縮,并通過對比臨界轉(zhuǎn)速和固有特性驗證了其模型的正確性。

    在螺栓連接轉(zhuǎn)子結(jié)構(gòu)的響應(yīng)特性分析方面,Qin等[12]推導(dǎo)了盤鼓連接結(jié)構(gòu)的彎曲剛度解析模型,將其引入轉(zhuǎn)子系統(tǒng)有限元模型中,進而采用諧波平衡法對轉(zhuǎn)子系統(tǒng)動力學(xué)方程進行求解。Klompas[13]采用彎矩來模擬螺栓連接結(jié)構(gòu),研究了連接結(jié)構(gòu)對渦輪轉(zhuǎn)子力學(xué)特性的影響。針對含螺栓連接Jeffcott轉(zhuǎn)子,Isa等[14]建立其集中參數(shù)模型,進而探究了系統(tǒng)動力學(xué)特性,基于該模型計算了轉(zhuǎn)子的不平衡響應(yīng),并開展了試驗研究??紤]航空發(fā)動機轉(zhuǎn)子盤鼓組合界面螺栓松動引起的時變剛度特性,Qin等[15]分析了螺栓松動對轉(zhuǎn)子系統(tǒng)動力學(xué)特性的影響規(guī)律。

    由上述文獻可知,目前,國內(nèi)外學(xué)者們針對轉(zhuǎn)子系統(tǒng)中螺栓連接結(jié)構(gòu)的研究依然很少,并且多數(shù)學(xué)者局限于對其固有特性的分析,對結(jié)構(gòu)的連接非線性以及響應(yīng)特性的研究較少[2,16],且并未考慮轉(zhuǎn)速效應(yīng)對其的影響??紤]到螺栓的接觸非線性是導(dǎo)致系統(tǒng)響應(yīng)特性求解效率低的關(guān)鍵因素,本文首先通過靜力學(xué)求得螺栓連接剛度,將彈簧單元引入到轉(zhuǎn)子系統(tǒng)中,進而將模型進行減縮處理,并通過對比實驗與仿真的固有特性結(jié)果,驗證所建模型的可靠性;其次,通過引入時變剛度考慮螺栓連接的非線性特性,探究螺栓擰緊個數(shù)、螺栓松動與系統(tǒng)轉(zhuǎn)速對結(jié)構(gòu)響應(yīng)特性的影響。

    2.2?螺栓個數(shù)的影響

    選擇預(yù)緊力F=800 N,螺栓松動個數(shù)m=0,轉(zhuǎn)速Ω=0,基于圖7 (a)螺栓連接結(jié)構(gòu)的時變剛度Kt求解流程,分析螺栓個數(shù)n=8,16,24時,螺栓連接結(jié)構(gòu)的時變剛度,提取0°360°螺栓連接的時變剛度,如圖8所示,為清晰看出時變剛度曲線的變化趨勢,選取n=16時的部分曲線進行放大處理。

    由圖8可知,當螺栓數(shù)目相同時,螺栓連接時變剛度存在波動現(xiàn)象,但變化趨勢不明顯;隨著螺栓數(shù)目的增加,時變剛度明顯增大,可以看出螺栓數(shù)目對系統(tǒng)的整體剛度影響較大。

    由圖9可以看出:激勵頻率fe相同時,隨著螺栓數(shù)目的增加,三維譜圖中倍頻成分不斷減少,非線性現(xiàn)象也越來越弱,這是因為螺栓數(shù)目越多,系統(tǒng)的連接界面的連接剛度越大,剛度波動減小,從而導(dǎo)致系統(tǒng)出現(xiàn)弱非線性;螺栓數(shù)目相同時,隨激勵頻率fe的增加,一倍頻幅值不斷增大;激勵頻率fe相同時,隨著螺栓數(shù)目的增加,一倍頻幅值不斷減小,這是因為螺栓數(shù)目越多,系統(tǒng)的整體剛度也越大,從而導(dǎo)致一倍頻幅值減小。

    2.3?螺栓松動的影響

    螺栓松動失效經(jīng)常發(fā)生,當螺栓發(fā)生松動時,螺栓連接界面的接觸特性與螺栓連接剛度均會發(fā)生變化,進而影響其動力學(xué)響應(yīng)?;谶@一現(xiàn)狀,本節(jié)分析不同螺栓松動個數(shù)對結(jié)構(gòu)時變剛度以及動力學(xué)響應(yīng)特性的影響。圖10為固定預(yù)緊力F=800 N,螺栓個數(shù)n=8,轉(zhuǎn)速Ω=0,螺栓松動個數(shù)m=0,1,2時,螺栓連接結(jié)構(gòu)的時變剛度變化曲線。

    從圖10可以看出:螺栓松動導(dǎo)致kx,ky,kθx和kθy四個方向的時變剛度曲線波動相對較大,kz和kθz兩個方向的剛度曲線波動很小,這主要是由于當m=1或2時,受z向力作用其松動位置與未松動位置的彈性變形差異較小,而結(jié)構(gòu)在剪切力作用下,當力與松動位置接近共線時會導(dǎo)致結(jié)構(gòu)發(fā)生較大變形;Mx與Fy,My與Fx對結(jié)構(gòu)的作用效果相同,因此剛度變化趨勢相近;而x與y方向上相差90°,導(dǎo)致kx和ky(kθx和kθy)兩個方向的時變剛度曲線波動趨勢相差90°的相位差。

    圖11給出了不同的螺栓松動個數(shù)(m=1, 2)和不同的激振頻率fe對鼓筒轉(zhuǎn)子系統(tǒng)振動響應(yīng)影響規(guī)律的三維譜圖。

    由圖11可以看出,激勵頻率fe相同時,隨著螺栓松動數(shù)目的增加,一倍頻幅值呈不斷增大趨勢,三維譜圖中倍頻成分不斷增多,幅值放大現(xiàn)象增多,非線性現(xiàn)象變強,尤其是當m=2時,出現(xiàn)了明顯的二倍頻(2fe),這是因為螺栓數(shù)目松動越多,系統(tǒng)的整體剛度越小,振動越劇烈,從而導(dǎo)致了結(jié)構(gòu)的強非線性。

    2.4?轉(zhuǎn)速的影響

    螺栓連接鼓筒轉(zhuǎn)子結(jié)構(gòu)在實際工作中多處于旋轉(zhuǎn)狀態(tài),其在旋轉(zhuǎn)過程中由于旋轉(zhuǎn)軟化、離心剛化和陀螺效應(yīng),螺栓連接的接觸狀態(tài)與連接剛度會受到一定影響。本小節(jié)主要討論預(yù)緊力F=800 N,螺栓個數(shù)n=8,螺栓松動個數(shù)m=2,轉(zhuǎn)速分別為Ω=0,3000,6000 r/min時,結(jié)合面時變剛度的變化。?圖12給出了不同轉(zhuǎn)速下系統(tǒng)時變剛度的變化曲線。

    從圖12可以看出:隨著系統(tǒng)轉(zhuǎn)速的增加,螺栓連接結(jié)構(gòu)的時變剛度不斷增加,剛度波動更加明顯,這是因為轉(zhuǎn)速越高,離心剛化的作用越明顯,其對結(jié)合面的接觸剛度和整體結(jié)構(gòu)剛度影響越大,并且轉(zhuǎn)速越高,對松動螺栓區(qū)域和非松動螺栓區(qū)域的剛度影響也越大,導(dǎo)致結(jié)合面的剛度波動較大。

    圖13給出了不同轉(zhuǎn)速和不同激振頻率fe對鼓筒轉(zhuǎn)子系統(tǒng)振動響應(yīng)影響規(guī)律的三維譜圖。

    由圖13可知,隨著轉(zhuǎn)速的增大,系統(tǒng)一倍頻幅值不斷減小,這是由于轉(zhuǎn)速加快導(dǎo)致連接剛度增大,響應(yīng)的幅值減小;隨著轉(zhuǎn)速的升高,二倍頻(2fe)變得更加明顯,這主要是由于轉(zhuǎn)速的升高導(dǎo)致連接處剛度波動較大,結(jié)構(gòu)非線性增強。

    3?結(jié)?論

    本文建立了螺栓連接鼓筒轉(zhuǎn)子結(jié)構(gòu)的減縮梁殼彈簧混合有限元模型,并基于該模型分析了螺栓個數(shù)、螺栓松動以及系統(tǒng)轉(zhuǎn)速對其時變剛度和響應(yīng)特性的影響,得出以下結(jié)論:

    (1) 本文提出的基于界面接觸求得螺栓連接剛度,并將其引入到減縮梁殼彈簧混合模型中的建模方法具有較高的精度和計算效率。

    (2) 隨著螺栓個數(shù)的增多,螺栓連接時變剛度增大;螺栓松動使螺栓連接的時變剛度波動較大,并且隨著螺栓松動個數(shù)的增多,結(jié)合面時變剛度減小;隨著系統(tǒng)轉(zhuǎn)速的增大,螺栓連接時變剛度不斷增大,剛度波動更加明顯。

    (3) 隨著螺栓連接個數(shù)的增多,系統(tǒng)非線性減弱;螺栓松動會導(dǎo)致系統(tǒng)非線性增強,且松動個數(shù)越多,系統(tǒng)非線性越強;轉(zhuǎn)速效應(yīng)在引起剛度波動增大的同時會使得系統(tǒng)非線性增強。

    參考文獻:

    [1]?秦朝燁, 王洪玉, 褚福磊. 螺栓連接對風扇轉(zhuǎn)子動力學(xué)特性影響分析[J]. 振動與沖擊,2012, 31(S):273276.

    Qin Zhaoye, Wang Hongyu, Chu Fulei. Influence of the bolted joint oil dynamics of fan rotor in aeroengine[J]. Journal of Vibration and Shock,2012, 31(S):273276.

    [2]?Liu S, Ma Y, Zhang D, et al. Studies on dynamic characteristics of the joint in the aeroengine rotor system[J]. Mechanical Systems & Signal Processing, 2012, 29(5):120136.

    [3]?Montgomery J. Methods for modeling bolts in the bolted joint[C]. ANSYS Users Conference, Orlando, Florida, 2002.

    [4]?Kim J, Yoon J C, Kang B S. Finite element analysis and modeling of structure with bolted joints [J]. Applied Mathematical Modelling, 2007, 31(5): 895911.

    [5]?趙?丹, 艾延廷, 翟?學(xué), 等. 盤盤螺栓連接結(jié)構(gòu)模態(tài)頻率分析[J]. 航空發(fā)動機, 2012, 38(5): 5557.

    Zhao Dan, Ai Yanting, Zhai Xue, et al. Analysis of modal frequencies for bolted plateplate structure[J]. Aeroengine, 2012, 38(5): 5557.

    [6]?Zhai X, Zhai Q G, Wang J J, et al. Research on parameterized modeling technology for the bolted joints structure[J]. Applied Mechanics and Materials,2014, 565: 211216.

    [7]?馬雙超. 航空發(fā)動機機匣模型確認與動力學(xué)特性研究[D].南京:南京航空航天大學(xué), 2012.

    Ma Shuangchao. Structure dynamic analysis and model validation of aeroengine casings[D]. Nanjing:Nanjing University of Aeronautics and Astronautics,2012.

    [8]?邱吉寶. 計算結(jié)構(gòu)動力學(xué)[M]. 合肥:中國科學(xué)技術(shù)大學(xué)出版社, 2009.

    Qiu Jibao. Computational Structural Dynamics[M]. Hefei:Press of University of Science and Technology of China, 2009.

    [9]?左彥飛, 王建軍, 馬威猛. 3D有限元轉(zhuǎn)子模型減縮的旋轉(zhuǎn)子結(jié)構(gòu)法[J]. 航空動力學(xué)報, 2014, 29(4):894900.

    Zuo Yanfei, Wang Jianjun, Ma Weimeng. Rotationg substructure method for 3D finite element rotor model reduction[J]. Journal of Aerospace Power, 2014, 29(4):894900.

    [10]左彥飛, 王建軍. 3D有限元轉(zhuǎn)子系統(tǒng)動力減縮的部件模態(tài)綜合方法及應(yīng)用[J]. 航空動力學(xué)報, 2016, 31(8):19291934.

    Zuo Yanfei, Wang Jianjun. Component mode synthesis for dynamic reduction of 3D finite element rotor system and its application[J]. Journal of Aerospace Power, 2016, 31(8):19291934.

    [11]孫傳宗, 陳予恕, 侯?磊. 復(fù)雜結(jié)構(gòu)雙轉(zhuǎn)子系統(tǒng)的建模及模型縮減[J]. 航空動力學(xué)報, 2017, 32(7):17471753.

    Sun Chuanzong, Chen Yushu, Hou Lei. Modeling method and reduction of dualrotor system with complicated structures[J]. Journal of Aerospace Power, 2017, 32(7):17471753.

    [12]Qin Z Y, Han Q K, Chu F L. Analytical model of bolted diskdrum joints and its application to dynamic analysis of jointed rotor[J].Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2014, 228(4): 646663.

    [13]Klompas N. Effects of anomalous rotor joints on turbomachine dynamics[J]. Journal of Engineering for Power, 1983, 105(4): 927934.

    [14]Mat Isa A A, Penny J E T, Garvey S D. Dynamics of bolted and laminated rotors[C]. Proceedings of the International Modal Analysis Conference, San Antonio, Texas, 2000, 1:867872.

    [15]Qin Z Y, Han Q K, Chu F L. Bolt loosening at rotating joint interface and its influence on rotor dynamics [J]. Engineering Failure Analysis, 2016, 59: 456466.

    [16]劉?存. 航空發(fā)動機轉(zhuǎn)子連接螺栓預(yù)緊力與疲勞壽命研究[D]. 南京:南京航空航天大學(xué), 2009.

    Liu Cun. Research on preload and fatigue life of connecting bolts in aeroengine rotors[D]. Nanjing:Nanjing University of Aeronautics and Astronautics,2009.

    [17]Herting D N. A general purpose, multistage, component modal synthesis method[J]. Finite Elements in Analysis and Design, 1985, 1(2):153164.

    Abstract: In this paper, a reduced beamshellspring mixed element model is established, and the natural characteristic of the system is analyzed based on the model. The accuracy of the model is verified by comparing the results with the experimental test results. Finally, the influence of the bolt number, the bolt looseness and the system speed on the timevarying stiffness and response characteristics of the bolted drum rotor structure are discussed. The results show that the number of bolt, the loosening of the bolt and the speed of the system have certain influences on the time varying stiffness and response characteristics of rotor system with bolt joints. With the increase of the number of bolts, the stiffness of the bolt connection will increase and the nonlinear characteristics of the system will be weak. The loosening of the bolt leads to the larger fluctuation of the connection stiffness, and as the number of loosened bolt increases, the nonlinearity of the system will increase. The stiffness of the bolt connection increases with the increase of the speed of the system, but the higher the speed is, the greater the fluctuation of the stiffness of the system will be, which will lead to the enhancement of nonlinear characteristics of the system.Key words: rotor dynamics; bolted joint; drum rotor; timevarying stiffness; bolt loosening

    作者簡介:孟春曉 (1993), 女, 碩士。電話: 13840342437; Email: 2494893725@qq.com

    通訊作者:馬?輝 (1978), 男, 教授。電話: (024)83684491; Email: mahui_2007@163.com

    两性夫妻黄色片| 国产男女内射视频| 久久久欧美国产精品| 精品福利观看| 91大片在线观看| 久久这里只有精品19| 精品国产乱码久久久久久男人| 久久天堂一区二区三区四区| 涩涩av久久男人的天堂| 国产又色又爽无遮挡免| 97精品久久久久久久久久精品| 欧美日韩av久久| 午夜福利影视在线免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲黑人精品在线| 亚洲精品粉嫩美女一区| 午夜福利一区二区在线看| 18禁观看日本| 妹子高潮喷水视频| 亚洲av片天天在线观看| 99久久国产精品久久久| 少妇裸体淫交视频免费看高清 | 久久久国产欧美日韩av| 五月天丁香电影| 天堂8中文在线网| 少妇人妻久久综合中文| 国产av精品麻豆| 久热这里只有精品99| 自线自在国产av| 国产一级毛片在线| 日韩制服骚丝袜av| 国产亚洲精品久久久久5区| 他把我摸到了高潮在线观看 | 999久久久精品免费观看国产| 久久精品亚洲熟妇少妇任你| 日本av手机在线免费观看| 日本一区二区免费在线视频| 久久精品亚洲av国产电影网| 性高湖久久久久久久久免费观看| 精品人妻1区二区| 水蜜桃什么品种好| 中文字幕高清在线视频| 色综合欧美亚洲国产小说| 国产色视频综合| 欧美日韩亚洲综合一区二区三区_| 两个人免费观看高清视频| 老司机深夜福利视频在线观看 | 久久久国产欧美日韩av| 亚洲欧美日韩另类电影网站| 久久99热这里只频精品6学生| 国产男人的电影天堂91| 精品熟女少妇八av免费久了| 狠狠狠狠99中文字幕| 欧美激情极品国产一区二区三区| 亚洲国产看品久久| 久久久久久久国产电影| 国产精品一二三区在线看| 在线观看舔阴道视频| 亚洲色图综合在线观看| 中国国产av一级| 夫妻午夜视频| 日本a在线网址| 国产真人三级小视频在线观看| av线在线观看网站| 成人av一区二区三区在线看 | 欧美激情极品国产一区二区三区| 欧美av亚洲av综合av国产av| 国产一级毛片在线| 侵犯人妻中文字幕一二三四区| 国产免费视频播放在线视频| 亚洲精品av麻豆狂野| 精品国产乱码久久久久久小说| 操美女的视频在线观看| 精品久久蜜臀av无| 亚洲avbb在线观看| 日本wwww免费看| 大型av网站在线播放| 91成年电影在线观看| 精品人妻熟女毛片av久久网站| 久久狼人影院| 丝袜喷水一区| 亚洲 国产 在线| 成年av动漫网址| 777久久人妻少妇嫩草av网站| 欧美日韩av久久| 久久人人爽人人片av| 老熟妇乱子伦视频在线观看 | 午夜两性在线视频| 国产精品熟女久久久久浪| 搡老熟女国产l中国老女人| 亚洲自偷自拍图片 自拍| 人妻一区二区av| 丁香六月天网| 91字幕亚洲| 激情视频va一区二区三区| 色婷婷久久久亚洲欧美| 国产免费一区二区三区四区乱码| 黄色a级毛片大全视频| √禁漫天堂资源中文www| 最新在线观看一区二区三区| 亚洲人成电影免费在线| 在线观看免费日韩欧美大片| 天堂中文最新版在线下载| 90打野战视频偷拍视频| 精品人妻熟女毛片av久久网站| 国产精品久久久人人做人人爽| 国产精品麻豆人妻色哟哟久久| 男人舔女人的私密视频| 国产精品久久久久久精品古装| 侵犯人妻中文字幕一二三四区| 亚洲精品中文字幕在线视频| 亚洲av成人不卡在线观看播放网 | 欧美黄色片欧美黄色片| 国产淫语在线视频| 久久精品熟女亚洲av麻豆精品| 女人高潮潮喷娇喘18禁视频| www.自偷自拍.com| 欧美+亚洲+日韩+国产| 91字幕亚洲| 精品人妻一区二区三区麻豆| 亚洲精品国产av蜜桃| 国产高清视频在线播放一区 | 成在线人永久免费视频| 91麻豆精品激情在线观看国产 | 国产精品影院久久| 精品少妇一区二区三区视频日本电影| 高清欧美精品videossex| 国产99久久九九免费精品| 久久影院123| 成人影院久久| 久久久精品94久久精品| 2018国产大陆天天弄谢| 少妇的丰满在线观看| 成人国产av品久久久| 精品久久久久久久毛片微露脸 | 国产精品偷伦视频观看了| 电影成人av| 最近最新免费中文字幕在线| 一级,二级,三级黄色视频| 少妇被粗大的猛进出69影院| 少妇被粗大的猛进出69影院| 男女床上黄色一级片免费看| 女性生殖器流出的白浆| 久久精品aⅴ一区二区三区四区| 午夜精品国产一区二区电影| 蜜桃国产av成人99| 久热爱精品视频在线9| 每晚都被弄得嗷嗷叫到高潮| 蜜桃国产av成人99| 国产精品国产三级国产专区5o| 9热在线视频观看99| 老司机在亚洲福利影院| 99国产精品一区二区三区| 久久久久久久久免费视频了| e午夜精品久久久久久久| bbb黄色大片| cao死你这个sao货| 欧美日本中文国产一区发布| 国产精品久久久av美女十八| 999久久久国产精品视频| av在线app专区| 久久中文看片网| 欧美激情 高清一区二区三区| 亚洲性夜色夜夜综合| 久久久久久久久久久久大奶| 国产成人一区二区三区免费视频网站| tocl精华| 国产色视频综合| 国产精品偷伦视频观看了| 久久久国产成人免费| 极品人妻少妇av视频| 91麻豆av在线| 精品一区在线观看国产| 黄色a级毛片大全视频| 色播在线永久视频| 色播在线永久视频| 水蜜桃什么品种好| 中文字幕制服av| 精品一品国产午夜福利视频| 18禁国产床啪视频网站| www.av在线官网国产| 国内毛片毛片毛片毛片毛片| 亚洲欧美一区二区三区黑人| 欧美日韩成人在线一区二区| 日韩一卡2卡3卡4卡2021年| 久久人妻熟女aⅴ| 精品国产乱码久久久久久男人| www.999成人在线观看| 国产亚洲欧美在线一区二区| 99精国产麻豆久久婷婷| 91av网站免费观看| cao死你这个sao货| 亚洲av国产av综合av卡| 中文欧美无线码| 天堂中文最新版在线下载| 亚洲五月婷婷丁香| 精品久久蜜臀av无| 日本猛色少妇xxxxx猛交久久| 久久影院123| 亚洲色图综合在线观看| 99久久99久久久精品蜜桃| 精品卡一卡二卡四卡免费| 在线十欧美十亚洲十日本专区| 久久人人爽人人片av| 飞空精品影院首页| 亚洲精品中文字幕在线视频| 日韩一区二区三区影片| 女人高潮潮喷娇喘18禁视频| 十八禁高潮呻吟视频| 美女午夜性视频免费| 欧美黑人精品巨大| 国产主播在线观看一区二区| 一本一本久久a久久精品综合妖精| 日韩,欧美,国产一区二区三区| 国产男女超爽视频在线观看| 午夜福利免费观看在线| 免费女性裸体啪啪无遮挡网站| av超薄肉色丝袜交足视频| 香蕉国产在线看| 日韩一卡2卡3卡4卡2021年| 少妇人妻久久综合中文| 欧美黄色淫秽网站| 高清在线国产一区| 久久久久久久精品精品| 精品久久久精品久久久| 9热在线视频观看99| 婷婷成人精品国产| 亚洲成人免费av在线播放| 久久久精品免费免费高清| 十八禁网站网址无遮挡| 日本黄色日本黄色录像| 黄色视频,在线免费观看| 男女边摸边吃奶| 精品国产乱码久久久久久小说| 首页视频小说图片口味搜索| 亚洲色图综合在线观看| 欧美精品亚洲一区二区| 精品亚洲成国产av| 久久狼人影院| 高清视频免费观看一区二区| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕人妻丝袜制服| 日韩有码中文字幕| 久久久国产一区二区| 亚洲国产精品一区三区| 中文字幕高清在线视频| a 毛片基地| 中文字幕制服av| 久久狼人影院| 黄色视频不卡| 满18在线观看网站| 黄片大片在线免费观看| 日本vs欧美在线观看视频| 久久久久精品国产欧美久久久 | 丝袜美足系列| 国产精品免费视频内射| 女人被躁到高潮嗷嗷叫费观| 国产亚洲av片在线观看秒播厂| 国产一卡二卡三卡精品| 一个人免费在线观看的高清视频 | 波多野结衣一区麻豆| 亚洲色图 男人天堂 中文字幕| 一本—道久久a久久精品蜜桃钙片| 淫妇啪啪啪对白视频 | 国产在线一区二区三区精| 日韩免费高清中文字幕av| 亚洲欧美精品综合一区二区三区| 中文字幕制服av| 午夜福利一区二区在线看| 久久久国产精品麻豆| 超色免费av| 大片电影免费在线观看免费| 91精品伊人久久大香线蕉| 秋霞在线观看毛片| 天天躁狠狠躁夜夜躁狠狠躁| 黄色视频不卡| 国产精品偷伦视频观看了| 大片电影免费在线观看免费| 国产在视频线精品| 日韩欧美国产一区二区入口| 欧美激情极品国产一区二区三区| 亚洲激情五月婷婷啪啪| 成人国语在线视频| 真人做人爱边吃奶动态| 亚洲精品成人av观看孕妇| 日韩 欧美 亚洲 中文字幕| 俄罗斯特黄特色一大片| 国产成人一区二区三区免费视频网站| xxxhd国产人妻xxx| av免费在线观看网站| 亚洲伊人色综图| 天天躁日日躁夜夜躁夜夜| 一本综合久久免费| 久久久久久久大尺度免费视频| videos熟女内射| 亚洲精品美女久久久久99蜜臀| 美女大奶头黄色视频| 国产欧美日韩综合在线一区二区| 伦理电影免费视频| 国产精品熟女久久久久浪| 成人国产一区最新在线观看| 欧美老熟妇乱子伦牲交| 国产免费福利视频在线观看| 热99re8久久精品国产| 天堂8中文在线网| 水蜜桃什么品种好| 精品亚洲乱码少妇综合久久| 欧美老熟妇乱子伦牲交| 在线观看免费日韩欧美大片| 可以免费在线观看a视频的电影网站| svipshipincom国产片| 中文字幕av电影在线播放| 18禁观看日本| 久久国产精品大桥未久av| 欧美精品人与动牲交sv欧美| 成年女人毛片免费观看观看9 | 黑人巨大精品欧美一区二区mp4| 男人添女人高潮全过程视频| 中文字幕人妻丝袜制服| 日韩一卡2卡3卡4卡2021年| 少妇猛男粗大的猛烈进出视频| 亚洲熟女毛片儿| 1024香蕉在线观看| 午夜福利在线观看吧| 国产高清国产精品国产三级| a 毛片基地| 久久国产精品男人的天堂亚洲| 精品久久久久久久毛片微露脸 | 黑人猛操日本美女一级片| 在线十欧美十亚洲十日本专区| av天堂在线播放| 性高湖久久久久久久久免费观看| 久久性视频一级片| 国产高清国产精品国产三级| 丝袜脚勾引网站| 日韩视频在线欧美| 精品免费久久久久久久清纯 | 91字幕亚洲| 在线天堂中文资源库| 国产一区二区三区av在线| 亚洲精品国产色婷婷电影| 手机成人av网站| 欧美97在线视频| 亚洲人成77777在线视频| 久久精品久久久久久噜噜老黄| 亚洲欧洲精品一区二区精品久久久| 亚洲国产精品成人久久小说| 国产又爽黄色视频| 狠狠婷婷综合久久久久久88av| 久久久久视频综合| 精品少妇内射三级| 亚洲三区欧美一区| 他把我摸到了高潮在线观看 | 国产免费av片在线观看野外av| 在线观看一区二区三区激情| 国产一区二区激情短视频 | 五月天丁香电影| 午夜影院在线不卡| 亚洲九九香蕉| 最新在线观看一区二区三区| 午夜91福利影院| 极品少妇高潮喷水抽搐| 一个人免费在线观看的高清视频 | 亚洲av电影在线观看一区二区三区| a 毛片基地| 在线观看免费视频网站a站| 男女免费视频国产| 人妻久久中文字幕网| 国产老妇伦熟女老妇高清| 国产在线观看jvid| 亚洲七黄色美女视频| 这个男人来自地球电影免费观看| 爱豆传媒免费全集在线观看| 在线观看免费高清a一片| 欧美 日韩 精品 国产| 12—13女人毛片做爰片一| 国产真人三级小视频在线观看| 两个人免费观看高清视频| 丝袜美腿诱惑在线| 中文字幕色久视频| 黄色视频在线播放观看不卡| 啦啦啦免费观看视频1| 黄色视频,在线免费观看| 免费在线观看黄色视频的| 日日摸夜夜添夜夜添小说| 99国产极品粉嫩在线观看| 女人久久www免费人成看片| 操出白浆在线播放| 99久久99久久久精品蜜桃| 日韩视频在线欧美| 国产成人影院久久av| av在线播放精品| 亚洲专区中文字幕在线| 俄罗斯特黄特色一大片| 久久久久国产一级毛片高清牌| 欧美 亚洲 国产 日韩一| 极品少妇高潮喷水抽搐| 久久人妻福利社区极品人妻图片| 欧美一级毛片孕妇| 国精品久久久久久国模美| 亚洲 国产 在线| 免费高清在线观看日韩| 精品国产一区二区久久| 欧美精品啪啪一区二区三区 | 亚洲欧美一区二区三区黑人| 伦理电影免费视频| 动漫黄色视频在线观看| 三上悠亚av全集在线观看| 中文字幕人妻丝袜一区二区| 在线av久久热| 午夜福利影视在线免费观看| 精品人妻在线不人妻| 欧美人与性动交α欧美精品济南到| 欧美日韩国产mv在线观看视频| 日韩电影二区| 色精品久久人妻99蜜桃| 搡老乐熟女国产| 亚洲精品国产区一区二| 国产亚洲欧美在线一区二区| 精品久久久久久电影网| 欧美性长视频在线观看| 欧美变态另类bdsm刘玥| 成年人黄色毛片网站| 精品高清国产在线一区| 日韩精品免费视频一区二区三区| 两性夫妻黄色片| 中文字幕高清在线视频| 国产一区有黄有色的免费视频| 久久久久国内视频| 男男h啪啪无遮挡| 人妻 亚洲 视频| bbb黄色大片| 一本综合久久免费| 午夜福利一区二区在线看| 两性夫妻黄色片| 中文字幕高清在线视频| 午夜影院在线不卡| 视频在线观看一区二区三区| 黄色怎么调成土黄色| 午夜福利视频在线观看免费| 1024视频免费在线观看| 亚洲精品av麻豆狂野| 人妻 亚洲 视频| 欧美黑人欧美精品刺激| 国产成人av教育| 女性被躁到高潮视频| 国产一卡二卡三卡精品| 国产成人av教育| 亚洲一区中文字幕在线| 日本精品一区二区三区蜜桃| 精品欧美一区二区三区在线| 悠悠久久av| 久久人妻熟女aⅴ| 人人澡人人妻人| 精品一区二区三卡| 欧美日韩福利视频一区二区| 色老头精品视频在线观看| 亚洲国产日韩一区二区| 夫妻午夜视频| 老熟妇仑乱视频hdxx| 久久99一区二区三区| 天堂8中文在线网| 午夜福利视频在线观看免费| 久久精品亚洲av国产电影网| 午夜日韩欧美国产| a级片在线免费高清观看视频| 最新在线观看一区二区三区| 建设人人有责人人尽责人人享有的| 国产一区二区激情短视频 | a级毛片黄视频| 亚洲国产日韩一区二区| 中文精品一卡2卡3卡4更新| 热re99久久精品国产66热6| 日韩视频一区二区在线观看| 多毛熟女@视频| 老熟女久久久| a在线观看视频网站| 18禁黄网站禁片午夜丰满| 亚洲av男天堂| 亚洲三区欧美一区| 黄色视频,在线免费观看| 成人影院久久| videos熟女内射| 后天国语完整版免费观看| 亚洲午夜精品一区,二区,三区| 五月开心婷婷网| 午夜激情久久久久久久| 亚洲精品国产av蜜桃| 午夜福利在线观看吧| 无限看片的www在线观看| 十八禁网站免费在线| 国产精品免费大片| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久亚洲精品国产蜜桃av| 麻豆乱淫一区二区| 精品少妇黑人巨大在线播放| 亚洲精品国产一区二区精华液| 另类精品久久| 久久精品国产亚洲av高清一级| 在线观看免费高清a一片| 久久亚洲国产成人精品v| 国产精品免费大片| 久久久国产欧美日韩av| 午夜老司机福利片| 12—13女人毛片做爰片一| 亚洲av日韩在线播放| 亚洲av成人一区二区三| 国产免费av片在线观看野外av| 亚洲免费av在线视频| 每晚都被弄得嗷嗷叫到高潮| 国产成人精品无人区| 午夜福利视频精品| 亚洲av美国av| 欧美另类一区| 国产真人三级小视频在线观看| 日本vs欧美在线观看视频| 国产一级毛片在线| 一级,二级,三级黄色视频| 无限看片的www在线观看| 19禁男女啪啪无遮挡网站| 伦理电影免费视频| 国产成人系列免费观看| 国产欧美日韩综合在线一区二区| 又紧又爽又黄一区二区| 大香蕉久久网| 伦理电影免费视频| 日本欧美视频一区| 精品一区二区三区四区五区乱码| 韩国高清视频一区二区三区| cao死你这个sao货| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产毛片av蜜桃av| 国产成人av教育| av线在线观看网站| 操美女的视频在线观看| 伊人久久大香线蕉亚洲五| 老司机午夜十八禁免费视频| 国产精品久久久av美女十八| 2018国产大陆天天弄谢| 久久久国产欧美日韩av| 日韩有码中文字幕| 动漫黄色视频在线观看| 老鸭窝网址在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 两性午夜刺激爽爽歪歪视频在线观看 | 一区二区三区激情视频| 女性生殖器流出的白浆| 亚洲精品乱久久久久久| 少妇粗大呻吟视频| 亚洲人成电影观看| 欧美激情久久久久久爽电影 | 大片免费播放器 马上看| 捣出白浆h1v1| 亚洲情色 制服丝袜| 亚洲视频免费观看视频| 久久久精品94久久精品| 精品一区二区三卡| 国产欧美日韩综合在线一区二区| 欧美精品亚洲一区二区| 久久久久久久久久久久大奶| 亚洲国产精品成人久久小说| 精品免费久久久久久久清纯 | 99国产精品免费福利视频| 亚洲精品国产色婷婷电影| h视频一区二区三区| av又黄又爽大尺度在线免费看| 精品亚洲乱码少妇综合久久| 18禁黄网站禁片午夜丰满| 美女视频免费永久观看网站| 一区二区三区精品91| 精品少妇一区二区三区视频日本电影| 大型av网站在线播放| 老司机影院成人| 国产伦人伦偷精品视频| 夜夜骑夜夜射夜夜干| 亚洲性夜色夜夜综合| 丝袜美足系列| 久久久国产精品麻豆| 国产免费视频播放在线视频| 热re99久久精品国产66热6| 免费观看av网站的网址| 一本久久精品| 久久人妻熟女aⅴ| 91九色精品人成在线观看| 国产一区二区三区在线臀色熟女 | 18禁裸乳无遮挡动漫免费视频| av在线app专区| 久久久久国产一级毛片高清牌| 国产又爽黄色视频| 色综合欧美亚洲国产小说| 国产又爽黄色视频| 大码成人一级视频| 美女中出高潮动态图| 久久久国产一区二区| 欧美97在线视频| 亚洲avbb在线观看| 午夜免费观看性视频| av一本久久久久| 深夜精品福利| 久久99热这里只频精品6学生| 十八禁网站免费在线| av国产精品久久久久影院| 久久久久精品人妻al黑| av在线播放精品| 大香蕉久久网| 丝袜美足系列| 黄色视频,在线免费观看| 国产成+人综合+亚洲专区| 免费观看av网站的网址| 热99国产精品久久久久久7| 国产成人精品久久二区二区91| 亚洲国产精品成人久久小说| 久久久久久久久久久久大奶| 亚洲专区国产一区二区| 欧美亚洲 丝袜 人妻 在线|