亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        考慮轉(zhuǎn)向架影響的重載貨車車廂模態(tài)特性*

        2015-06-13 09:38:06張琪昌于躍斌呂大力
        振動、測試與診斷 2015年6期
        關(guān)鍵詞:剛體側(cè)墻轉(zhuǎn)向架

        王 婧,張琪昌,于躍斌,呂大力,靳 剛,趙 峰

        (1.天津大學天津市非線性動力學與混沌控制重點實驗室 天津,300072)(2.齊齊哈爾軌道交通裝備有限責任公司 齊齊哈爾,161002)

        ?

        考慮轉(zhuǎn)向架影響的重載貨車車廂模態(tài)特性*

        王 婧1,張琪昌1,于躍斌2,呂大力2,靳 剛1,趙 峰1

        (1.天津大學天津市非線性動力學與混沌控制重點實驗室 天津,300072)(2.齊齊哈爾軌道交通裝備有限責任公司 齊齊哈爾,161002)

        為深入研究轉(zhuǎn)向架對重載貨車動力學特性的影響,以某空車工況的40t軸重礦石車為研究對象,在實車裝配模型的基礎(chǔ)上建立了整車及車廂的力學仿真模型,并應用有限元法對二者進行了模態(tài)分析。對實際車輛進行了模態(tài)實驗研究,通過計算結(jié)果及實驗結(jié)果的對比,從剛體振型及彈性體振型兩個方面考察了轉(zhuǎn)向架對重載貨車車廂模態(tài)特性的影響。研究表明,轉(zhuǎn)向架的質(zhì)量、剛度及邊界條件均會對模態(tài)計算結(jié)果產(chǎn)生影響,整車模型能更精確地描述實車的模態(tài)特性。

        結(jié)構(gòu)振動; 模態(tài)分析; 模態(tài)實驗; 有限元分析; ANSYS軟件;鐵路重載貨車

        引 言

        重載快捷運輸是世界鐵路發(fā)展的一大趨勢,鐵路的客、貨運能力很大程度上決定著國民經(jīng)濟和社會的發(fā)展水平。為實現(xiàn)鐵路重載快捷運輸?shù)目缭绞桨l(fā)展,需要不斷提高貨車車輛的設(shè)計研發(fā)和制造技術(shù)水平。其中,結(jié)合車輛的動力學特性進行車體結(jié)構(gòu)的設(shè)計及優(yōu)化是亟待解決的核心問題[1],而開展模態(tài)計算是進行動力學優(yōu)化設(shè)計的重要前提。

        模態(tài)分析是動力學設(shè)計及優(yōu)化的基礎(chǔ),國內(nèi)外學者在相關(guān)領(lǐng)域進行了較為深入的研究。張大鈞等[2]采用隨機沖擊激振方法和分區(qū)激振測試技術(shù)對北京型內(nèi)燃機車車體進行了實驗模態(tài)分析。有限元模擬方面,湯禮鵬等[3]用非線性接觸單元模擬心盤、旁承與車體的連接形式,以梁單元和彈簧單元模擬承載。郭志全等[4]通過彈簧阻尼單元及節(jié)點耦合的方式將車體和轉(zhuǎn)向架連結(jié)為一體。王卉子[5]用集中力表示轉(zhuǎn)向架對車體的作用,但并沒有考慮旁承對車體的支撐影響。文獻[6-8]將轉(zhuǎn)向架構(gòu)架強度與疲勞計算相結(jié)合進行分析,然而涉及車體與轉(zhuǎn)向架的結(jié)合研究較少。文獻[9]的車輛模型中省略了轉(zhuǎn)向架,用剛性處理替代基本連接區(qū),計算結(jié)果與實際情況誤差較大。常寧等[10-12]用梁單元和彈簧單元模擬承載方式。以上方法在計算時均為針對車廂模型的計算,并未建立完整的整車(包括車廂、轉(zhuǎn)向架及輪對)模型,因此無法真實地反應出車體與轉(zhuǎn)向架之間的動態(tài)特性[13],使得計算精度難以保證[14]。

        隨著鐵路運輸?shù)目旖莼l(fā)展,重載列車的結(jié)構(gòu)更趨于復雜。在轉(zhuǎn)向架和車體組成的多自由度耦合動力系統(tǒng)中,轉(zhuǎn)向架對車體結(jié)構(gòu)的影響無法忽略[15-18]。在鐵路工程計算中,大多采用分別計算車體和轉(zhuǎn)向架的模態(tài),并據(jù)此進行動力學設(shè)計的方法。筆者以某型40t軸重礦石車為研究對象,將整車、車廂模型的有限元計算結(jié)果與實驗結(jié)果進行對比分析,討論轉(zhuǎn)向架對車輛模態(tài)的影響,對提高計算精度及復雜結(jié)構(gòu)重載貨車的動力學設(shè)計有著重要的參考及借鑒意義。

        1 車輛模型的動力學模態(tài)計算分析

        1.1 有限元模型的建立

        筆者以某型40t軸重礦石車實車裝配圖紙為依據(jù),根據(jù)車體結(jié)構(gòu)、裝配和約束特點,將車輛的裝配模型轉(zhuǎn)化為力學模型,精簡掉對整車的應力分布和力學特性影響較小的構(gòu)件,建立了車輛結(jié)構(gòu)的動力學有限元模型,并以此為基礎(chǔ)進行了車輛系統(tǒng)的模態(tài)分析。

        在Hypermesh中對車廂和轉(zhuǎn)向架兩個模型進行組裝,在磨耗板與斜楔間建立可滑動接觸以模擬側(cè)架對搖枕的作用,二系減振通過軸向8根彈簧(總軸向剛度為9 612N/mm)和橫向2根彈簧(總橫向剛度為6 489N/mm)來模擬。一系減振的橡膠塊以實體單元劃分,采用Link8單元模擬底架與底板的連接,并對交界面添加接觸邊界條件。將旁承部位的約束等效為垂向彈簧,釋放心盤部位繞垂向軸的旋轉(zhuǎn)自由度。車體采用Shell63單元模擬,賦予不同的實常數(shù)表征不同厚度的車體構(gòu)造板材。圖1為整車的有限元模型,約束后輪與鐵軌接觸部分的所有自由度,前車輪約束除縱向位移外的其余自由度,單元數(shù)為159 154。圖2為車廂的有限元模型,該模型心盤部分除繞垂向軸的旋轉(zhuǎn)自由度外,固定約束其余自由度,并考慮旁承對車廂的垂向支撐作用[19](總支撐剛度為8 800N/mm),單元數(shù)為72 915。各部分材料參數(shù)如表1所示。

        圖1 整車有限元模型

        圖2 車廂有限元模型

        表1 各部分材料參數(shù)

        Tab.1 The material parameter

        部件彈性模量/kPa泊松比密度/(kg·mm-3)單元類型搖枕、側(cè)架1.75×1080.37.85×10-6 Solid45輪對2.06×1080.37.85×10-6 Solid45一系減振2.1×1080.37.85×10-6 Solid45二系減振2.1×1080.37.85×10-6 Combin14車體1.75×1080.37.85×10-6 Shell63底架2.06×1080.37.85×10-6 Shell63焊接點2.06×1080.37.85×10-6 Link8沖擊座、上心盤1.75×1080.37.85×10-6 Solid45

        1.2 模態(tài)計算結(jié)果

        應用Block Lanczos模態(tài)提取法對整車模型及車廂模型進行模態(tài)分析,得到了前10階固有頻率及振型,如表2所示。由于重點關(guān)注車廂的模態(tài),因此在計算及實驗結(jié)果中略去了轉(zhuǎn)向架的局部模態(tài)。

        表2 模態(tài)計算結(jié)果

        Tab.2 The results of modal calculation

        階數(shù)振型整車模型頻率/Hz車廂模型頻率/Hz1剛體-車體繞縱向軸轉(zhuǎn)動4.27213.8192剛體-車體繞垂向軸轉(zhuǎn)動5.987—3沿縱向中心軸扭轉(zhuǎn)11.42715.0954側(cè)墻1階反向彎曲22.40019.9235側(cè)墻2階反向彎曲29.59525.3046側(cè)墻1階同向彎曲28.27132.2997側(cè)墻2階同向彎曲36.47040.0638側(cè)墻3階反向彎曲,端墻1階反向彎曲39.77037.7779側(cè)墻3階同向彎曲,端墻1階同向彎曲46.53049.26510側(cè)墻3階反向彎曲,端墻1階反向彎曲49.76050.394

        同時,以車體繞縱向軸轉(zhuǎn)動、側(cè)墻1階反向彎曲、側(cè)墻2階同向彎曲振型為例,給出整車模型的振型圖,如圖3所示。

        圖3 整車模型的振型圖

        2 實車動力學實驗

        為驗證有限元計算結(jié)果,對轉(zhuǎn)向架支撐的實車進行了振動實驗研究。實驗采用單點激勵多點拾振的方法,模態(tài)實驗系統(tǒng)如圖4所示。將車輛沿其縱軸方向分成7個截面,每個截面設(shè)置7個測點,在前后兩轉(zhuǎn)向架的側(cè)梁、搖枕分別設(shè)置4個測點,并在車輛兩端面的上梁的中點分別設(shè)置4個測點,測點數(shù)共77個,分布情況如圖5所示。采用沖擊激勵的方式激發(fā)實驗車進入振動狀態(tài),每個測點測試x,y,z三個方向的響應信號,實驗激勵點如圖5所示,分析帶寬為40Hz,測得實驗結(jié)果如表3所示。

        圖4 模態(tài)實驗系統(tǒng)

        圖5 實驗測點布置示意圖

        表3 模態(tài)實驗結(jié)果

        Tab.3 The results of modal testing

        階數(shù)振型實驗結(jié)果頻率/Hz1剛體-車體繞縱向軸轉(zhuǎn)動4.8432剛體-車體繞垂向軸轉(zhuǎn)動6.8623沿縱向中心軸扭轉(zhuǎn)11.7814側(cè)墻1階反向彎曲22.5605側(cè)墻2階反向彎曲31.7096側(cè)墻1階同向彎曲32.6767側(cè)墻2階同向彎曲38.7168側(cè)墻3階反向彎曲,端墻1階反向彎曲41.8809側(cè)墻3階同向彎曲,端墻1階同向彎曲47.67510側(cè)墻3階反向彎曲,端墻1階反向彎曲49.407

        3 轉(zhuǎn)向架對模態(tài)計算的影響

        在結(jié)構(gòu)的動力計算中,帶有轉(zhuǎn)向架的模型結(jié)構(gòu)更為復雜,相比于無轉(zhuǎn)向架模型,計算量大為增加。因此,筆者研究了轉(zhuǎn)向架對結(jié)構(gòu)模態(tài)計算的影響,并探討不同模型的適用性。

        3.1 轉(zhuǎn)向架對剛體振型的影響

        以實驗結(jié)果中的兩階剛體振型為基準,將兩種模型的計算結(jié)果和實驗結(jié)果進行比較,可以發(fā)現(xiàn):

        1) 實驗及整車模型計算中均出現(xiàn)2階剛體振型,而由于未考慮轉(zhuǎn)向架的影響,車廂模型計算結(jié)果中僅出現(xiàn)1階;

        2) 整車模型計算所得剛體振型的出現(xiàn)順序與實驗結(jié)果一致,說明此模型與實際工況相符,而車廂模型過于簡化,未出現(xiàn)第1階剛體振型;

        3) 車廂模型的計算誤差大于整車模型的誤差,由表2中的計算結(jié)果及表4中的誤差分析可知,在車廂模型的計算結(jié)果中未出現(xiàn)第2階剛體振型(車體繞垂向軸轉(zhuǎn)動),而整車模型則可以得到全部的剛體振型,且模型可修改性較好。

        表4 各模型剛體振型頻率計算結(jié)果誤差

        Tab.4 Errors of rigid modal calculation

        階數(shù)振型整車模型誤差/%車廂模型誤差/%1剛體-車體繞縱向軸轉(zhuǎn)動-13.37185.342剛體-車體繞垂向軸轉(zhuǎn)動-12.75—

        根據(jù)以上數(shù)據(jù)分析可以看出,整車模型計算所得振型跟實驗一致性較好,而且顯著提高了剛體振型的計算精度。其原因如下:

        1) 雖然兩個模型中車廂部分的約束相同,但整車模型在模態(tài)計算時考慮了轉(zhuǎn)向架的高度,而車廂模型未加以考慮,造成車廂模型重心偏低,因此頻率偏大;

        2) 整車模型中車廂下方是與轉(zhuǎn)向架相連接的,轉(zhuǎn)向架本身的減振系統(tǒng)增大了結(jié)構(gòu)的阻尼,且轉(zhuǎn)向架的減振系統(tǒng)與車廂構(gòu)成了耦合的彈性體系,而車廂模型下方?jīng)]有減振結(jié)構(gòu),相比于整車模型而言其剛度更大,因此在模態(tài)計算中車廂模型的振動頻率更大。

        3.2 轉(zhuǎn)向架對彈性體振型的影響

        以實驗結(jié)果為基準,從以下幾個方面對數(shù)據(jù)進行了分析。

        1) 模態(tài)的完整性:整車及車廂模型均能得到完整的彈性體模態(tài)振型。

        2) 振型順序:實驗與模型計算所得彈性體振型出現(xiàn)順序基本一致。

        3) 頻率計算精度:整車模型的頻率誤差在0.71%~13.48%之間,如表5所示,說明對實車模型的簡化是合理的;車廂模型中彈性體振型頻率誤差在1%~28%之間,大于整車模型的誤差,說明車廂模型的精確性需要進一步探討。

        表5 各模型彈性體振型頻率計算結(jié)果誤差表

        Tab.5 Errors table of elastic modal calculation

        階數(shù)振型整車模型誤差/%車廂模型誤差/%3沿縱向中心軸扭轉(zhuǎn)-3.0028.134側(cè)墻1階反向彎曲-0.71-11.695側(cè)墻2階反向彎曲-6.67-20.206側(cè)墻1階同向彎曲-13.48-1.157側(cè)墻2階同向彎曲-5.803.488側(cè)墻3階反向彎曲,端墻1階反向彎曲-5.04-9.809側(cè)墻3階同向彎曲,端墻1階同向彎曲-2.403.3410側(cè)墻3階反向彎曲,端墻1階反向彎曲0.712.00方差0.00170.0182

        由表5中的方差值可以看出,就整體計算精度而言,整車模型與實際運行中的車輛更為接近。綜上所述,整車模型的計算結(jié)果與實驗結(jié)果的一致性較好,車廂模型的誤差較大。由此可見,轉(zhuǎn)向架對模態(tài)的計算精度有較大的影響,考慮轉(zhuǎn)向架的影響因素對探索車輛動力學特性具有重要的意義。因此,若不關(guān)心剛體振型,傳統(tǒng)的車廂模型單元數(shù)最少,具有較高的計算效率;若想提高計算精度,則需考慮轉(zhuǎn)向架的影響,整車模型不僅可以獲得完整、可靠的剛體振型,同時也提高了彈性體模態(tài)振型的計算精度。

        4 結(jié) 論

        1) 轉(zhuǎn)向架對貨車車廂的剛體模態(tài)振型及彈性體模態(tài)振型均有影響。剛體振型方面,轉(zhuǎn)向架與車廂存在振動耦合現(xiàn)象,對轉(zhuǎn)向架的簡化使車廂模型重心偏低,約束剛度偏大,因此車廂模型的剛體模態(tài)振型計算結(jié)果誤差較大;彈性體模態(tài)方面,轉(zhuǎn)向架的剛度及質(zhì)量均對模態(tài)頻率產(chǎn)生影響,整車模型的精確度更高。因此,建立整車模型是提高計算精度的必要條件。

        2) 整車有限元模型的模態(tài)計算結(jié)果與實驗結(jié)果吻合度較好,可為工程計算中類似模型的處理提供參考,也為深入了解車輛結(jié)構(gòu)系統(tǒng)的動力學特性提供了幫助。根據(jù)筆者的模態(tài)分析結(jié)果可進行結(jié)構(gòu)動力學修改,以進一步提高車輛的運動特性。

        [1] Dickens J M, Nakagawa J M, Wittbrodt M J. A critique of mode acceleration and modal truncation augmentation methods for modal response analysis[J]. Computers & Structures, 1997, 62(6): 985-998.

        [2] 張大鈞,蕭龍翔,霍拳忠. 北京型內(nèi)燃機車車體的試驗模態(tài)分析[J]. 鐵道學報, 1990, 12(2): 1-9.

        Zhang Dajun, Xiao Longxiang, Huo Quanzhong. Experimental modal analysis of Beijing diesel locomotive carbody[J]. Journal of the China Railway Society, 1990, 12(2): 1-9. (in Chinese)

        [3] 湯禮鵬,劉國偉. 30t軸重浴盆式敞車車體模態(tài)分析[J]. 機械, 2010, 37(11): 3-5.

        Tang Lipeng, Liu Guowei. Modal analysis of 30 t-axleload bathtub type gondola carbody[J]. Machinery, 2010, 37(11): 3-5. (in Chinese)

        [4] 郭志全,徐燕申,楊江天. 基于FEM的新型運煤敞車的結(jié)構(gòu)模態(tài)分析[J]. 機械強度, 2006, 28(6): 919-922.

        Guo Zhiquan, Xu Yanshen, Yang Jiangtian. Modal analyses based on FEM of open boxcar[J]. Journal of Mechanical Strength, 2006, 28(6): 919-922. (in Chinese)

        [5] 王卉子. C80型重載貨車車體結(jié)構(gòu)的CAE分析[D]. 大連: 大連交通大學, 2006.

        [6] 任萬勇,羅冠煒. 轉(zhuǎn)8A貨車轉(zhuǎn)向架重要承載部件的模態(tài)試驗與分析[J]. 工程力學, 1999, 3(3): 848-852.

        Ren Wanyong, Luo Guanwei. Modal test and analysis on important bearing part of zhuan 8A freight car bogies[J]. Engineering Mechanics, 1999, 3(3): 848-852. (in Chinese)

        [7] 魏培虎,張濟民,王開文. 貨車轉(zhuǎn)向架動力學參數(shù)測試方法研究[J]. 鐵道車輛, 2002, 40(6): 16-18.

        Wei Peihu, Zhang Jimin, Wang Kaiwen. Research on measurement method for dynamics parameters of freight car bogies[J]. Rolling Stock, 2002, 40(6): 16-18. (in Chinese)

        [8] 彭永明,封全保,李強. HX_D2B機車轉(zhuǎn)向架構(gòu)架有限元分析[J]. 鐵道機車車輛, 2010, 30(3): 42-44.

        Peng Yongming, Feng Quanbao, Li Qiang. Finite-element analyses on the bogie frame for HX_D2B locomotive[J]. Railway Locomotive & Car, 2010, 30(3): 42-44. (in Chinese)

        [9] 楊愛國,張志強,楊江天. 基于有限元建模的敞車輕量化設(shè)計[J]. 中國鐵道科學, 2007, 28(3): 79-83.

        Yang Aiguo, Zhang Zhiqiang, Yang Jiangtian. Lightweight design of open top wagon using updated finite element models[J]. China Railway Science, 2007, 28(3): 79-83. (in Chinese)

        [10]常寧,姚曙光,劉國偉. 基于有限元的重載敞車車體結(jié)構(gòu)靜力分析[C]∥2008年(第十屆)中國科協(xié)年會論文集. 北京: 國防工業(yè)出版社,2008:588-595.

        [11]趙俊杰,魯寨軍,徐程. 新型108 t通用敞車車體結(jié)構(gòu)強度分析[J]. 機車車輛工藝, 2006, 2(1): 6-8.

        Zhao Junjie, Lu Zhaijun, Xu Cheng. Structural strength analysis of the new 108t general-purpose open-top car body[J]. Locomotive and Rolling Stock Technology, 2006, 2(1): 6-8. (in Chinese)

        [12]常寧,姚曙光,劉國偉. 軸重40t重載敞車車體有限元分析[J]. 鐵道機車車輛, 2009, 29(1): 8-12.

        Chang Ning, Yao Shuguang, Liu Guowei. FEM analysis of 40 t axleload heavy haul gondola car-body[J]. Railway Locomotive & Car, 2009, 29(1): 8-12. (in Chinese)

        [13]曲天威,王惠玉. 機車車體和轉(zhuǎn)向架模態(tài)分析研究[J]. 鐵道機車車輛, 2012, 32(3): 5-8.

        Qu Tianwei, Wang Huiyu. Modal analysis of locomotive body and bogies[J].Railway Locomotive & Car, 2012, 32(3): 5-8. (in Chinese)

        [14]顧楓. SS1型電力機車轉(zhuǎn)向架模態(tài)分析研究[J]. 機車電傳動, 1996(6): 14-18.

        Gu Feng. Modal study of bogie frame of SS1 electric locomotive[J]. Electric Drive for Locomotive, 1996(6): 14-18. (in Chinese)

        [15]宮島,周勁松,孫文靜. 高速列車彈性車體與轉(zhuǎn)向架耦合振動分析[J]. 交通運輸工程學報, 2011, 11(4): 41-47.

        Gong Dao, Zhou Jinsong, Sun Wenjing. Coupled vibration analysis of flexible car body and bogie for high-speed train[J]. Journal Traffic and Transportation Engineering, 2011, 11(4): 41-47. (in Chinese)

        [16]Dietz S, Netter H, Sachau D. Fatigue life prediction of a railway bogie under dynamic loads through simulation[J]. Vehicle System Dynamics, 1998, 29(6): 385-402.

        [17]鄔平波,薛世海,楊晨輝.基于彈性車體模型的高速客車動態(tài)響應[J].交通運輸工程學報, 2005, 5(2): 5-8.

        Wu Pingbo, Xue Shihai, Yang Chenhui. Dynamic response of high-speed passenger car based on flexible car body model[J]. Journal Traffic and Transportation Engineering, 2005, 5(2): 5-8. (in Chinese)

        [18]孫玉華,董大偉,閆兵. 雙層隔振系統(tǒng)模態(tài)匹配分析及其振動特性[J].振動、測試與診斷, 2014,34(4):727-731.

        Sun Yuhua, Dong Dawei, Yan Bing. Modal matching analysis and vibration characteristics of two-stage vibration isolation system[J]. Journal of Vibration, Measurement & Diagnosis, 2014,34(4):727-731.(in Chinese)

        [19]孫保衛(wèi),王紅,商躍進. 車輛彈性旁承垂向剛度分析方法研究[J].蘭州交通大學學報,2009, 28(1): 130-133.

        Sun Baowei, Wang Hong, Shang Yuejin. A study of the vertical stiffness of vehicle elastic side bearing[J]. Journal of Lanzhou Jiaotong University, 2009, 28(1): 130-133. (in Chinese)

        10.16450/j.cnki.issn.1004-6801.2015.06.028

        *天津市自然科學基金重點資助項目(12JCZDJC28000)

        2014-06-09;

        2014-10-31

        U272.2

        王婧,女,1989年1月生,博士研究生。主要研究方向為非線性振動及模態(tài)研究。 E-mail:guu120@163.com

        猜你喜歡
        剛體側(cè)墻轉(zhuǎn)向架
        基于ABAQUS的動車組側(cè)墻滑槽連接結(jié)構(gòu)設(shè)計優(yōu)化
        差值法巧求剛體轉(zhuǎn)動慣量
        20t軸重米軌轉(zhuǎn)向架裝用搖枕的研制
        基于SPS模式的轉(zhuǎn)向架軸箱裝配線仿真研究
        基于結(jié)構(gòu)應力的轉(zhuǎn)向架構(gòu)架焊縫疲勞強度研究
        地鐵站側(cè)墻鋼筋混凝土結(jié)構(gòu)保護層厚度控制
        建筑科技(2018年6期)2018-08-30 03:40:58
        城軌不銹鋼車輛鼓型側(cè)墻弧焊工藝研究
        車載冷發(fā)射系統(tǒng)多剛體動力學快速仿真研究
        209 P型轉(zhuǎn)向架軸箱彈簧斷裂的原因及預防
        MIG—V工作站在高速動車鋁合金車體側(cè)墻焊接中的應用
        焊接(2015年1期)2015-07-18 11:07:33
        天天干夜夜躁| 久久久亚洲av成人网站| av国产传媒精品免费| 欧美 国产 综合 欧美 视频| 蜜臀久久99精品久久久久久小说| 亚洲国产AV无码男人的天堂| 国产成人九九精品二区三区| 无码日韩AⅤ一区二区三区| 国产免费无码9191精品| 97人妻蜜臀中文字幕| 国产高清在线精品一区二区三区| 亚洲国产综合久久天堂| 无码av中文一区二区三区桃花岛 | 国产成本人片无码免费2020| 国产精品无套内射迪丽热巴| 国产精品99久久久久久98AV| 亚洲日韩精品AⅤ片无码富二代| 天堂影院久久精品国产午夜18禁 | 45岁妇女草逼视频播放| 后入到高潮免费观看| 亚洲а∨天堂久久精品2021| 无码国产精品一区二区免费16| 青青草视频华人绿色在线 | 一区二区三区中文字幕脱狱者 | 少妇无码av无码专区线| 亚洲AV日韩AV永久无码电影| 91超碰在线观看免费| 冲田杏梨av天堂一区二区三区| 国产美腿丝袜一区二区| 欧美人与善在线com| 无码人妻精品丰满熟妇区| 免费人成黄页在线观看视频国产| 久久久久成人精品免费播放| 国产成人亚洲合色婷婷| 日韩高清不卡一区二区三区| 在线观看老湿视频福利| 日韩视频中文字幕精品偷拍 | 青青草伊人视频在线观看| 色和尚色视频在线看网站| 国产v片在线播放免费无码| 亚洲av无码一区二区三区四区 |