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

    基于多體動力學(xué)的柴油機(jī)曲軸疲勞壽命分析

    2017-05-11 11:53:10常麗田俊龍
    車用發(fā)動機(jī) 2017年2期
    關(guān)鍵詞:軸頸曲軸邊界條件

    常麗, 田俊龍

    (1. 太原科技大學(xué)機(jī)械學(xué)院, 山西 太原 030024; 2. 太原理工大學(xué)機(jī)械學(xué)院, 山西 太原 030024)

    基于多體動力學(xué)的柴油機(jī)曲軸疲勞壽命分析

    常麗1, 田俊龍2

    (1. 太原科技大學(xué)機(jī)械學(xué)院, 山西 太原 030024; 2. 太原理工大學(xué)機(jī)械學(xué)院, 山西 太原 030024)

    為分析4100QBZL柴油機(jī)曲軸的疲勞壽命,建立該曲柄連桿機(jī)構(gòu)的剛?cè)狁詈隙囿w動力學(xué)模型,將多組試驗測量的缸內(nèi)壓力作為驅(qū)動力,進(jìn)行耦合仿真得到曲軸在柔性體模型下的主軸頸、連桿軸頸負(fù)荷仿真結(jié)果,并根據(jù)載荷結(jié)果對曲軸進(jìn)行靜強(qiáng)度校核。最后結(jié)合由多體動力學(xué)軟件得到的載荷譜與有限元分析所得的曲軸在各個工況下的應(yīng)力應(yīng)變分析結(jié)果,以及通過材料的各項屬性擬合出的S-N曲線,對曲軸進(jìn)行了疲勞壽命預(yù)測。結(jié)果表明:曲軸的靜強(qiáng)度及疲勞壽命均達(dá)到了工程設(shè)計要求,曲軸最危險部位的壽命次數(shù)也達(dá)到了1013以上,認(rèn)為曲軸不會發(fā)生疲勞破壞。

    柴油機(jī); 曲軸; 多體動力學(xué); 疲勞壽命

    在國家經(jīng)濟(jì)從快速增長轉(zhuǎn)變?yōu)榉€(wěn)步增長的背景下,2015年度,我國汽車產(chǎn)量仍然以2 450萬輛的產(chǎn)量排在世界首位。但相比于汽車制造的繁榮,我們在汽車的設(shè)計以及零部件的關(guān)鍵技術(shù)上還存在差距,與此同時,在所生產(chǎn)產(chǎn)品的性能、質(zhì)量穩(wěn)定性以及可靠性上與國外相比也有較大的差距[1]。曲軸作為發(fā)動機(jī)的關(guān)鍵零部件,是發(fā)動機(jī)的主要運動件,具有結(jié)構(gòu)復(fù)雜、承受載荷大、載荷頻率高的特點。由于曲軸在工作時長時間承受交變載荷,所以曲軸的失效形式主要是疲勞破壞[2]。柴油機(jī)的可靠與否很大程度上取決于曲軸是否可靠,對曲軸疲勞產(chǎn)生原因的分析以及疲勞壽命的研究就顯得尤為重要。

    1 曲軸系統(tǒng)的多體動力學(xué)分析

    1.1 建立曲柄連桿機(jī)構(gòu)的多體動力學(xué)模型

    通過軟件對曲柄連桿機(jī)構(gòu)的相關(guān)零件,包括飛輪、曲軸、活塞、活塞銷、連桿、軸承組、軸瓦等,進(jìn)行了三維模型的建立,并將所有零件裝配組成曲軸系統(tǒng)的幾何模型[3](見圖1)。發(fā)動機(jī)相關(guān)參數(shù)見表1。

    活塞材料為鋁,密度為2 700 kg/m3;曲軸的材料為鋼,密度為7 850 kg/m3,泊松比為0.29,彈性模量為2.07×1011Pa。

    圖1 曲軸連桿系剛體動力學(xué)模型

    缸徑/mm100行程/mm105氣缸數(shù)4燃燒室型式直噴ω型燃燒室排量/L3.298壓縮比17.5∶1標(biāo)定功率/kW70標(biāo)定轉(zhuǎn)速/r·min-13200最大扭矩/N·m245最大扭矩轉(zhuǎn)速/r·min-12400供油提前角/(°)14±2

    1.2 邊界條件的建立

    1.2.1 運動邊界條件的建立

    添加約束:取缸體為基礎(chǔ),主軸軸承與缸體為固定副,軸承與曲軸為轉(zhuǎn)動副,連桿軸瓦與曲軸以球鉸鏈代替,軸瓦與連桿為轉(zhuǎn)動副,連桿軸瓦與連桿為轉(zhuǎn)動副,活塞銷與連桿為轉(zhuǎn)動副,活塞與活塞銷為轉(zhuǎn)動副,活塞與氣缸為平動副。

    添加驅(qū)動:為了最真實地模擬發(fā)動機(jī)的運行情況,在靠近飛輪端的主軸頸與曲軸之間的旋轉(zhuǎn)副上添加驅(qū)動,驅(qū)動方程為曲軸轉(zhuǎn)角與時間的函數(shù)[4]。

    1.2.2 工況邊界條件的確定

    試驗用發(fā)動機(jī)為4100QBZL增壓中冷柴油發(fā)動機(jī),發(fā)動機(jī)通過ET2000測控系統(tǒng)調(diào)控,缸內(nèi)壓力傳感器為Kistler 6125B型,壓力傳感器通過在氣缸蓋頂部打孔的方式安裝,數(shù)據(jù)采集設(shè)備為YE6232B動態(tài)數(shù)據(jù)采集器。

    試驗時以0號柴油為燃料,在1 600 r/min,2 000 r/min,2 400 r/min的轉(zhuǎn)速下通過電渦流測功機(jī)進(jìn)行加載,分別以10 N·m,100 N·m,180 N·m負(fù)載代表發(fā)動機(jī)的輕載、中載和重載,測定發(fā)動機(jī)各缸的缸內(nèi)壓力曲線。

    采用試驗測得的缸內(nèi)壓力作為主動力驅(qū)動模型,取0.1°曲軸轉(zhuǎn)角為步長,測得曲軸旋轉(zhuǎn)720°,即發(fā)動機(jī)一個工作循環(huán)的7 200個數(shù)據(jù),并在ADAMS/View中利用Akima樣條擬合函數(shù)將壓力數(shù)據(jù)樣條擬合為一條逼近數(shù)據(jù)點的光滑曲線,即缸內(nèi)壓力樣條曲線Spline。圖2示出1 600 r/min轉(zhuǎn)速100 N·m負(fù)載時的壓力曲線。此發(fā)動機(jī)點火次序為1—3—4—2,根據(jù)工作順序?qū)訔l曲線依次加載到各缸活塞表面。氣體爆發(fā)壓力以體積力的形式加載在每個活塞表面上計算得出集中力載荷[5]。

    函數(shù)方程:

    AKISPL(Measure angle,0,Spline,0)。

    其中,Measure angle為曲軸旋轉(zhuǎn)角度的測量函數(shù),變化范圍0°~720°。

    圖2 缸內(nèi)壓力樣條曲線

    1.3 曲柄連桿機(jī)構(gòu)的剛?cè)狁詈戏抡?/p>

    將原有的剛性體曲軸進(jìn)行柔性化并替換[6],得到曲柄連桿機(jī)構(gòu)的剛?cè)狁詈夏P?,并測試靜平衡,測試通過后開始動力學(xué)仿真。取發(fā)動機(jī)常用轉(zhuǎn)速1 600 r/min,仿真步數(shù)為1 000步,時間0.07 s,仿真結(jié)束后進(jìn)入后處理模塊,分析仿真結(jié)果。

    圖3示出連桿軸頸在剛?cè)狁詈夏P拖碌氖芰?,圖4示出主軸頸在剛?cè)狁詈夏P拖碌氖芰ΑGS作為長柔性體桿件,在外載荷作用下,會發(fā)生彎曲以及扭轉(zhuǎn)形變。從圖3可以看出,連桿軸頸的載荷幅值大約在25 kN。由圖4可以看出,第一至第五主軸頸在發(fā)動機(jī)一個工作循環(huán)內(nèi)載荷的變化基本相似,載荷譜幅值在50 kN以內(nèi),突變的激振力載荷介于100 kN到150 kN之間。可以看出,曲軸在實際工作中情況非常復(fù)雜。由剛?cè)狁詈夏P头抡娴贸鲋鬏S頸與連桿軸頸的最大激振力,這是曲軸動態(tài)特性研究的重要邊界條件。

    圖3 剛?cè)狁詈夏P拖逻B桿軸頸在X方向上的載荷

    圖4 剛?cè)狁詈夏P拖轮鬏S頸在X方向上的載荷

    2 曲軸的靜力學(xué)分析

    發(fā)動機(jī)在運行過程中,曲軸承受著周期性變化的燃燒壓力、往復(fù)質(zhì)量慣性力、旋轉(zhuǎn)運動質(zhì)量慣性力及其力矩(轉(zhuǎn)矩以及彎矩),4100QBZL柴油機(jī)的工作順序為1—3—4—2,通過試驗測出發(fā)動機(jī)實際運轉(zhuǎn)時的燃燒壓力數(shù)據(jù),在進(jìn)行多體動力學(xué)仿真時,氣體力作為主動力以體積力的形式加載在氣缸內(nèi)。將多體動力學(xué)仿真得到的曲軸主軸頸、連桿軸頸的集中力作為靜力學(xué)分析的主要載荷[7]。

    2.1 邊界條件的建立

    2.1.1 位移邊界條件的建立

    以Solid185對曲軸進(jìn)行自由網(wǎng)格劃分,在飛輪端主軸頸圓周面約束其所有節(jié)點X,Y,Z3個方向的移動自由度,自由端軸頸圓周面所有節(jié)點約束其X,Y方向的移動自由度。

    2.1.2 載荷邊界條件的建立

    根據(jù)有限寬度油膜壓力的分布特點,可假設(shè)邊界條件為載荷沿著軸線方向按照二次拋物線分布,在沿著軸的徑向方向按照余弦規(guī)律分布[8](見圖5)。

    圖5 軸頸受力邊界條件

    曲軸沿著軸頸方向載荷方程:

    式中:Qc為作用在軸頸上的集中力;R為軸半徑;x為軸頸承受載荷的長度。

    2.1.3 重力以及旋轉(zhuǎn)慣性力的加載

    定義曲軸密度為7 850kg/m3,重力加速度g=9.8m/s2,取發(fā)動機(jī)轉(zhuǎn)速為1 600r/min,方向同曲軸實際運轉(zhuǎn)方向相同,旋轉(zhuǎn)慣性力會自動加載到每個節(jié)點上。

    曲軸進(jìn)行限定邊界條件之后,根據(jù)有限寬度油膜分布規(guī)律,分別以發(fā)動機(jī)每缸最高燃燒壓力下的載荷值對連桿軸頸進(jìn)行加載,加約束和載荷后的曲軸見圖6。

    圖6 加約束和載荷后的曲軸

    2.2 靜力學(xué)計算結(jié)果分析

    對力和位移邊界條件約束后進(jìn)入仿真環(huán)境,運行結(jié)束后通過后處理查看應(yīng)力云圖。

    圖7至圖10示出了柴油機(jī)在一個工作循環(huán)內(nèi)的應(yīng)力分布云圖,圖11示出最大應(yīng)力局部圖。從圖11可以看出,曲軸的最大應(yīng)力均位于主軸頸與曲柄的過渡圓角處或連桿軸頸與曲柄的過渡圓角處。當(dāng)?shù)谝桓装l(fā)火和第二缸發(fā)火時,最大應(yīng)力位于連桿軸頸與曲柄臂的連接處,應(yīng)力值分別為23.75 MPa和33.67 MPa。第三缸發(fā)火時,最大應(yīng)力位于連桿與軸頸接觸點的正下方,應(yīng)力為50.63 MPa。第四缸發(fā)火時最大應(yīng)力位于相鄰自由端軸頸與曲柄的連接處,最大應(yīng)力為66.71 MPa。曲軸所受應(yīng)力最大和應(yīng)力最為集中的區(qū)域,也是曲軸最容易發(fā)生破壞的地方,在疲勞壽命預(yù)測時需要進(jìn)行重點分析。

    圖7 第一缸發(fā)火時的整體應(yīng)力云圖

    圖8 第二缸發(fā)火時的整體應(yīng)力云圖

    圖9 第三缸發(fā)火時的整體應(yīng)力云圖

    圖10 第四缸發(fā)火時的整體應(yīng)力云圖

    圖11 最大應(yīng)力局部圖

    通過曲軸的靜力學(xué)分析,利用多體動力學(xué)仿真的載荷譜作為力的邊界約束條件,曲軸在發(fā)動機(jī)中的實際約束條件作為位移邊界條件,進(jìn)行加載,計算出曲軸的應(yīng)力最大處以及應(yīng)力值,為疲勞壽命分析提供數(shù)據(jù)支持。

    3 曲軸的疲勞壽命分析

    在靜態(tài)載荷分析的基礎(chǔ)上,對柴油機(jī)曲軸進(jìn)行疲勞強(qiáng)度分析,并根據(jù)疲勞壽命云圖找出曲軸的最危險部位,從而推算出曲軸的整體使用壽命[9]。

    曲軸疲勞是典型的高周疲勞[10],所以采用工程上最常用的基于線性疲勞累積損傷理論的名義應(yīng)力法,對曲軸進(jìn)行疲勞壽命分析。

    3.1 疲勞壽命的分析流程

    疲勞壽命的分析流程(見圖12)一般分為以下幾個步驟:

    1) 載荷譜的獲取,采用試驗測試或者仿真模擬分析的方法獲得載荷數(shù)據(jù);

    2) 應(yīng)力和應(yīng)變載荷的獲??;

    3) 疲勞分析模型的確定,這是疲勞分析的核心步驟,選擇適合的疲勞壽命分析模型,并根據(jù)疲勞損傷累積方法,結(jié)合材料參數(shù)以及載荷數(shù)據(jù),進(jìn)行疲勞壽命預(yù)測;

    4) 輸入材料參數(shù);

    5) 進(jìn)行疲勞分析并查看結(jié)果。

    圖12 疲勞壽命分析模塊

    3.2 曲軸的疲勞壽命計算

    1) 有限元模型的導(dǎo)入

    將有限元分析的應(yīng)力應(yīng)變結(jié)果文件載入到DesignLife中,顯示當(dāng)前分析結(jié)果。

    2) 載荷數(shù)據(jù)導(dǎo)入

    進(jìn)入ADAMS后處理,將載荷數(shù)據(jù)導(dǎo)出至Designlife,載荷譜見圖13。

    圖13 疲勞載荷曲線

    3) 設(shè)置材料的S-N曲線

    可以通過試驗或公式擬合得到材料的S-N曲線,由于曲軸的疲勞破壞均為高周疲勞,通過疲勞試驗測定材料的疲勞極限會耗費大量的人力物力,故采用公式擬合法確定其S-N曲線[11-12]。

    4100QBZL柴油機(jī)曲軸材料為42CrMo,其主要成分見表2。

    表2 曲軸材料成分

    將表3中材料參數(shù)導(dǎo)入后,建立材料的S-N曲線,設(shè)置完成之后,得到該材料經(jīng)UTS修正的S-N曲線(見圖14)。

    表3 42CrMo材料特性

    圖14 曲軸材料的S-N曲線

    3.3 仿真結(jié)果分析

    在完成所有參數(shù)設(shè)置后開始仿真,結(jié)束進(jìn)入后處理階段。查看后處理模塊窗口,疲勞壽命分析的結(jié)果見圖15。

    圖15 曲軸最低壽命區(qū)域的局部應(yīng)力云圖

    由圖15可知,曲軸的疲勞壽命最低的部位在6 966節(jié)點,節(jié)點壽命為8.835×1013,該節(jié)點位于靠近飛輪端第一連桿軸頸與連桿的接觸位置,與ANSYS中分析的應(yīng)力最大點位置基本吻合。其次為自由端連桿軸頸與主軸頸的過渡圓角處,節(jié)點壽命為4.453×1014。根據(jù)壽命損傷云圖,該曲軸所有節(jié)點壽命均高于1013。載荷時間歷程為30 s,城市平均時速60 km/h,由此估算出該曲軸設(shè)計滿足行駛里程百萬公里以上的要求,認(rèn)為該零件在整個運行生命周期內(nèi)不會發(fā)生破壞。

    4 結(jié)束語

    建立曲柄連桿機(jī)構(gòu)的剛?cè)狁詈夏P?,得到曲軸在穩(wěn)態(tài)過程中的各個主軸頸,連桿軸頸的壓力載荷譜,為進(jìn)一步的曲軸靜力學(xué)研究提供載荷邊界條件。對曲軸進(jìn)行靜力學(xué)分析,利用多體動力學(xué)仿真的載荷譜作為力的邊界約束條件,曲軸在發(fā)動機(jī)中的實際約束條件作為位移邊界條件,計算出曲軸的應(yīng)力最大處以及應(yīng)力值,為疲勞壽命分析提供數(shù)據(jù)支持。對4100柴油機(jī)作了強(qiáng)度校核,根據(jù)Mine疲勞累積損傷理論,結(jié)合多體動力學(xué)軟件得到的載荷譜與有限元分析所得的曲軸在各個工況下的應(yīng)力應(yīng)變分析結(jié)果,以及通過材料的各項屬性擬合出的S-N曲線,對曲軸作了疲勞壽命預(yù)測,分析結(jié)果顯示,曲軸所有節(jié)點壽命均高于1013,滿足行駛里程百萬公里以上的要求,認(rèn)為該零件在整個運行生命周期內(nèi)不會發(fā)生破壞。

    [1] 門峰,王今.中國汽車產(chǎn)業(yè)發(fā)展趨勢預(yù)測[J].汽車工業(yè)研究,2011(2):2-5.

    [2] Aleksandar Vencl.Aleksandar Rac Diesel engine crankshaft journal bearings failures: Case study[J].Engineering Failure Analysis,2014(44):217-228.

    [3] 左鵬,武志斐,王鐵,等.F-T柴油摻燒甲醇對柴油機(jī)燃燒及振動特性的影響[J].內(nèi)燃機(jī)工程,2016,37(5):27-31.

    [4] 段秀兵,郝志勇,岳東鵬,等.汽車發(fā)動機(jī)曲軸扭振的多體動力學(xué)分析[J].汽車工程,2005(2):233-235,250.

    [5] 戴曦.8L250型柴油機(jī)曲軸動力學(xué)仿真及強(qiáng)度校核[D].武漢:武漢理工大學(xué),2012.

    [6] Deshmukh S V, Wayzode R R, Alvi N G. Dynamic Load and Stress Analysis of a Crankshaft[J]. Indian Streams Research journal,2012,Ⅱ(Ⅳ):765.

    [7] 王勝.基于UG/ANSYS的曲軸參數(shù)化建模系統(tǒng)設(shè)計及有限元分析[D].成都:西南交通大學(xué),2012.

    [8] Kharabe S,Ukhande M,Mane R,et al. An Evaluation of Crankshaft Bending Rigidity Using Finite Element Analysis[C].SAE Paper 2013-01-0126.

    [9] Wang Z M, Lv P M. Crankshaft Fatigue Test and Statistical Distribution of Fatigue Experimental Data[J]. Applied Mechanics & Materials,2013,401-403:93-96.

    [10] 陳淵博,郝志勇,張煥宇.基于彎曲疲勞試驗的柴油機(jī)曲軸疲勞壽命分析及改進(jìn)[J].內(nèi)燃機(jī)工程,2011,32(1):75-78,84.

    [11] Utkin V S.Calculation of crankshaft reliability in terms of fatigue strength, with limited statistical data[J].Russian Engineering Research,2010,30(8):763-767.

    [12] 孫楠楠,李國祥,白書戰(zhàn),等.基于應(yīng)變-壽命理論的曲軸疲勞分析研究[J]. 內(nèi)燃機(jī)工程,2014,35(6):60-64,83.

    [編輯: 潘麗麗]

    Fatigue Life Analysis of Diesel Engine Crankshaft Based on Multibody Dynamics

    CHANG Li1, TIAN Junlong2

    (1. School of Mechanical Engineering, Taiyuan University of Science and Technology, Taiyuan 030024, China; 2. School of Mechanical Engineering, Taiyuan University of Technology, Taiyuan 030024, China)

    In order to analyze the fatigue life of 4100QBZL diesel engine crankshaft, the coupled multibody dynamics model of crankshaft system was established. The load simulation results of main journal and connecting rod journal under the flexible body model of crankshaft were acquired and the static strength was then checked. Finally, the crankshaft fatigue life prediction was conducted by combining the load spectrum collected from multibody dynamics software, the stress and strain finite element analysis results of the crankshaft under various conditions and the fitted SN curve according to the material properties. The results show that both the crankshaft static strength and fatigue life can meet the design requirements of the project. The life of critical position for crankshaft is beyond 1013times, so it is thought that crankshaft fatigue failure will not happen.

    diesel engine; crankshaft; multibody dynamics; fatigue life

    2016-05-13;

    2017-02-17

    常麗(1989—),男,碩士,研究方向為發(fā)動機(jī)關(guān)鍵零部件研究;changli0504@qq.com。

    10.3969/j.issn.1001-2222.2017.02.013

    TK423.3

    B

    1001-2222(2017)02-0072-06

    猜你喜歡
    軸頸曲軸邊界條件
    一類帶有Stieltjes積分邊界條件的分?jǐn)?shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動邊值問題的漸近解
    淺析曲軸磨床的可靠性改進(jìn)
    一種曲軸豎置熱處理裝置設(shè)計
    曲軸軸頸磨削變形的可疊加機(jī)理分析
    基于支持向量機(jī)的曲軸動平衡控制方法
    曲軸連桿軸頸表面鍍覆層的改性效果比較
    曲軸軸頸車-車梳刀具工作原理及結(jié)構(gòu)設(shè)計
    應(yīng)用數(shù)控加工技術(shù)提高曲軸軸頸的粗加工精度
    帶Robin邊界條件的2維隨機(jī)Ginzburg-Landau方程的吸引子
    亚洲av成人精品一二三区| 在现免费观看毛片| 女的被弄到高潮叫床怎么办| 免费黄频网站在线观看国产| 一个人看视频在线观看www免费| 国产高潮美女av| 亚洲欧洲国产日韩| 伦理电影免费视频| 久久午夜福利片| 狂野欧美白嫩少妇大欣赏| 国产一区有黄有色的免费视频| 国产在线视频一区二区| 国产精品偷伦视频观看了| 国产亚洲欧美精品永久| 久久99热这里只有精品18| 亚洲怡红院男人天堂| 寂寞人妻少妇视频99o| 国产毛片在线视频| 不卡视频在线观看欧美| 亚洲欧美日韩东京热| 久久久久久久久久久免费av| 欧美xxⅹ黑人| 国产成人精品久久久久久| 国产成人精品久久久久久| 这个男人来自地球电影免费观看 | 欧美国产精品一级二级三级 | 日本爱情动作片www.在线观看| 中文乱码字字幕精品一区二区三区| 久久久久久久国产电影| 人人妻人人爽人人添夜夜欢视频 | 久久精品国产鲁丝片午夜精品| 亚洲国产精品999| 夜夜看夜夜爽夜夜摸| 又大又黄又爽视频免费| 亚洲成人手机| 伊人久久国产一区二区| 人妻制服诱惑在线中文字幕| 少妇人妻一区二区三区视频| av国产精品久久久久影院| 国产精品一二三区在线看| 精品一品国产午夜福利视频| 99视频精品全部免费 在线| av国产免费在线观看| 一级a做视频免费观看| 夜夜骑夜夜射夜夜干| 九九爱精品视频在线观看| 99re6热这里在线精品视频| 这个男人来自地球电影免费观看 | 黄片无遮挡物在线观看| 嫩草影院新地址| 99热这里只有是精品在线观看| a级一级毛片免费在线观看| 肉色欧美久久久久久久蜜桃| 国产精品人妻久久久影院| 精品久久久精品久久久| 中文在线观看免费www的网站| 免费大片黄手机在线观看| 国产精品偷伦视频观看了| 噜噜噜噜噜久久久久久91| 精品人妻熟女av久视频| 亚洲av中文字字幕乱码综合| 国产极品天堂在线| 视频中文字幕在线观看| 国产一级毛片在线| www.av在线官网国产| av播播在线观看一区| 久久ye,这里只有精品| 日韩 亚洲 欧美在线| 美女中出高潮动态图| 国产av精品麻豆| 国产精品偷伦视频观看了| 美女xxoo啪啪120秒动态图| 亚洲欧美一区二区三区国产| 91狼人影院| 国产亚洲5aaaaa淫片| 国产亚洲91精品色在线| 亚洲国产色片| 大又大粗又爽又黄少妇毛片口| 高清午夜精品一区二区三区| 久久久国产一区二区| 一本一本综合久久| 精品一品国产午夜福利视频| 在线观看免费日韩欧美大片 | 18禁裸乳无遮挡动漫免费视频| av国产久精品久网站免费入址| 九九爱精品视频在线观看| 三级经典国产精品| 99久久中文字幕三级久久日本| 99久久综合免费| 久久精品国产亚洲av天美| 成人免费观看视频高清| 亚洲精品色激情综合| 在现免费观看毛片| 亚洲欧美成人精品一区二区| 国产极品天堂在线| 美女xxoo啪啪120秒动态图| 99热这里只有是精品在线观看| 人人妻人人看人人澡| 日本-黄色视频高清免费观看| 亚洲av.av天堂| 一级二级三级毛片免费看| 麻豆乱淫一区二区| 日韩欧美一区视频在线观看 | 国产精品欧美亚洲77777| 亚洲国产色片| 欧美人与善性xxx| 日本午夜av视频| 观看美女的网站| 一级av片app| 免费看光身美女| 成人黄色视频免费在线看| 久久99热这里只频精品6学生| 青春草国产在线视频| 尾随美女入室| 自拍偷自拍亚洲精品老妇| 综合色丁香网| 啦啦啦啦在线视频资源| 国产白丝娇喘喷水9色精品| 国产精品av视频在线免费观看| 91精品国产国语对白视频| 免费黄网站久久成人精品| www.av在线官网国产| 日本wwww免费看| 99久久综合免费| 97在线人人人人妻| 久久热精品热| 亚洲色图综合在线观看| 国产日韩欧美亚洲二区| 国产高潮美女av| 精品午夜福利在线看| 高清欧美精品videossex| 久久精品国产亚洲av天美| 久久毛片免费看一区二区三区| 99热这里只有精品一区| 久久久久性生活片| 国产精品伦人一区二区| 欧美97在线视频| 少妇高潮的动态图| 亚洲精品亚洲一区二区| 日韩伦理黄色片| 免费av中文字幕在线| 丰满少妇做爰视频| 热re99久久精品国产66热6| 成人一区二区视频在线观看| 成人亚洲精品一区在线观看 | 观看美女的网站| 中文乱码字字幕精品一区二区三区| 99久久人妻综合| 人妻制服诱惑在线中文字幕| 欧美亚洲 丝袜 人妻 在线| 丰满人妻一区二区三区视频av| 99久久综合免费| tube8黄色片| 寂寞人妻少妇视频99o| 多毛熟女@视频| 三级国产精品片| 国产在线视频一区二区| 精品午夜福利在线看| 亚洲欧美一区二区三区国产| 伦理电影大哥的女人| 国产成人免费观看mmmm| 国产人妻一区二区三区在| 久久精品熟女亚洲av麻豆精品| 日韩av不卡免费在线播放| 久久久久人妻精品一区果冻| 国产精品一二三区在线看| 国产成人精品福利久久| 国产亚洲91精品色在线| 纯流量卡能插随身wifi吗| 国产高清不卡午夜福利| av一本久久久久| 特大巨黑吊av在线直播| 永久网站在线| 久久久久久久精品精品| 夜夜爽夜夜爽视频| 亚洲熟女精品中文字幕| 热99国产精品久久久久久7| 精品视频人人做人人爽| 五月伊人婷婷丁香| 精品久久久精品久久久| 国产伦理片在线播放av一区| 久久久成人免费电影| 女性被躁到高潮视频| 国产成人午夜福利电影在线观看| 日本wwww免费看| 老司机影院成人| 成人午夜精彩视频在线观看| 亚洲成人中文字幕在线播放| 99久久精品国产国产毛片| 久久这里有精品视频免费| 日日啪夜夜撸| 精品久久久精品久久久| h视频一区二区三区| av播播在线观看一区| 六月丁香七月| 亚洲三级黄色毛片| 国产精品无大码| 欧美成人a在线观看| tube8黄色片| 久久久午夜欧美精品| 最近中文字幕高清免费大全6| av黄色大香蕉| 只有这里有精品99| 亚洲欧美清纯卡通| 人妻一区二区av| 一区二区av电影网| 国产乱来视频区| 九草在线视频观看| 色哟哟·www| 国产无遮挡羞羞视频在线观看| 黄色日韩在线| 中国国产av一级| 精品99又大又爽又粗少妇毛片| 亚洲电影在线观看av| av一本久久久久| 九九爱精品视频在线观看| 边亲边吃奶的免费视频| 国产av国产精品国产| 国产一区有黄有色的免费视频| 日产精品乱码卡一卡2卡三| 91精品国产国语对白视频| 国产淫语在线视频| 日本vs欧美在线观看视频 | 在线观看免费视频网站a站| av播播在线观看一区| 丰满少妇做爰视频| 一个人看视频在线观看www免费| 亚洲精品乱久久久久久| 国产91av在线免费观看| 2022亚洲国产成人精品| 热99国产精品久久久久久7| 精品少妇黑人巨大在线播放| 精品久久久久久久久av| 精品人妻一区二区三区麻豆| 国产 一区精品| 美女cb高潮喷水在线观看| 中文资源天堂在线| 中文精品一卡2卡3卡4更新| 国产69精品久久久久777片| 五月开心婷婷网| 国产国拍精品亚洲av在线观看| 欧美3d第一页| 欧美精品一区二区大全| 美女内射精品一级片tv| 国产伦在线观看视频一区| 亚洲av.av天堂| 日韩一区二区三区影片| 亚洲第一区二区三区不卡| 美女xxoo啪啪120秒动态图| 老师上课跳d突然被开到最大视频| 久久99精品国语久久久| 国产片特级美女逼逼视频| 国产成人freesex在线| 一本色道久久久久久精品综合| 色视频www国产| 黄色日韩在线| 视频中文字幕在线观看| 日韩视频在线欧美| 日本wwww免费看| 国产中年淑女户外野战色| 伦理电影大哥的女人| 国国产精品蜜臀av免费| av在线老鸭窝| 视频区图区小说| 精品人妻熟女av久视频| 你懂的网址亚洲精品在线观看| 亚洲国产精品一区三区| 免费人成在线观看视频色| 亚洲成人中文字幕在线播放| 亚洲精品第二区| 国产精品一区二区在线观看99| 91久久精品国产一区二区三区| 国产深夜福利视频在线观看| 日本午夜av视频| 国产美女午夜福利| 肉色欧美久久久久久久蜜桃| 亚洲aⅴ乱码一区二区在线播放| 大又大粗又爽又黄少妇毛片口| 美女主播在线视频| 一级黄片播放器| 精品国产一区二区三区久久久樱花 | 午夜视频国产福利| .国产精品久久| 婷婷色综合大香蕉| 男女边吃奶边做爰视频| 亚洲精品日本国产第一区| 亚洲人成网站高清观看| 成人毛片60女人毛片免费| 丰满乱子伦码专区| 三级国产精品片| 国产精品99久久久久久久久| 最黄视频免费看| 91在线精品国自产拍蜜月| 国产深夜福利视频在线观看| 肉色欧美久久久久久久蜜桃| 国产国拍精品亚洲av在线观看| 欧美xxxx性猛交bbbb| 亚洲精华国产精华液的使用体验| 一区二区三区精品91| 精品久久久精品久久久| 女人久久www免费人成看片| 成人影院久久| 久久精品熟女亚洲av麻豆精品| 日本wwww免费看| 亚洲精品第二区| 人体艺术视频欧美日本| 国产成人免费无遮挡视频| 80岁老熟妇乱子伦牲交| 一边亲一边摸免费视频| 精品久久久精品久久久| 久久国产精品大桥未久av | 麻豆乱淫一区二区| 婷婷色麻豆天堂久久| 久久久久网色| 成人毛片60女人毛片免费| 亚洲va在线va天堂va国产| 女性生殖器流出的白浆| 午夜福利视频精品| 国产在线免费精品| 久久国产亚洲av麻豆专区| 91狼人影院| 亚洲精品国产色婷婷电影| 日韩视频在线欧美| 一本—道久久a久久精品蜜桃钙片| 亚洲精品乱码久久久v下载方式| 欧美日韩在线观看h| 亚洲精品日本国产第一区| 全区人妻精品视频| 日韩免费高清中文字幕av| 一个人看的www免费观看视频| 国产精品女同一区二区软件| 国产伦理片在线播放av一区| 一级毛片久久久久久久久女| 99re6热这里在线精品视频| 青春草亚洲视频在线观看| 精品一品国产午夜福利视频| 免费久久久久久久精品成人欧美视频 | 不卡视频在线观看欧美| 最近中文字幕高清免费大全6| 亚洲无线观看免费| 超碰av人人做人人爽久久| 人人妻人人看人人澡| 欧美国产精品一级二级三级 | 午夜福利在线观看免费完整高清在| 精品久久久久久久末码| 中文字幕制服av| 午夜视频国产福利| 国产女主播在线喷水免费视频网站| 美女福利国产在线 | 亚洲国产高清在线一区二区三| 国产亚洲最大av| 99久久精品一区二区三区| 男女下面进入的视频免费午夜| 久久久久网色| 亚洲av欧美aⅴ国产| 亚洲色图综合在线观看| 免费看光身美女| 亚洲欧洲日产国产| 中文在线观看免费www的网站| www.色视频.com| 美女高潮的动态| 国产av码专区亚洲av| 美女内射精品一级片tv| 伦精品一区二区三区| 丰满人妻一区二区三区视频av| 涩涩av久久男人的天堂| 久久精品久久久久久噜噜老黄| 91精品伊人久久大香线蕉| 99国产精品免费福利视频| 精品午夜福利在线看| 狠狠精品人妻久久久久久综合| 深夜a级毛片| 2021少妇久久久久久久久久久| 日韩国内少妇激情av| 超碰av人人做人人爽久久| av在线app专区| 国内少妇人妻偷人精品xxx网站| 免费久久久久久久精品成人欧美视频 | 久久久久久人妻| 成年美女黄网站色视频大全免费 | 国产一区二区三区av在线| 麻豆精品久久久久久蜜桃| 一个人看的www免费观看视频| 六月丁香七月| 小蜜桃在线观看免费完整版高清| 成人特级av手机在线观看| 久久久久久人妻| 国产精品av视频在线免费观看| 妹子高潮喷水视频| 亚洲精品日韩av片在线观看| 99热全是精品| 亚洲美女黄色视频免费看| 91精品伊人久久大香线蕉| 伦理电影大哥的女人| 欧美性感艳星| 国产欧美亚洲国产| 国产视频首页在线观看| 欧美人与善性xxx| 色哟哟·www| 另类亚洲欧美激情| 777米奇影视久久| 日韩电影二区| 国产精品人妻久久久影院| 美女中出高潮动态图| 亚洲国产精品专区欧美| 乱码一卡2卡4卡精品| 免费观看a级毛片全部| 男人狂女人下面高潮的视频| 久久精品人妻少妇| 干丝袜人妻中文字幕| 22中文网久久字幕| 插阴视频在线观看视频| 国产在线免费精品| 精品久久久久久久末码| 免费观看av网站的网址| 少妇猛男粗大的猛烈进出视频| av福利片在线观看| 卡戴珊不雅视频在线播放| 久久久久久久久久人人人人人人| 亚洲色图综合在线观看| 九九在线视频观看精品| 久久人人爽av亚洲精品天堂 | 舔av片在线| 尾随美女入室| 一级毛片电影观看| 伊人久久精品亚洲午夜| 成人漫画全彩无遮挡| av黄色大香蕉| 日日啪夜夜撸| 国产91av在线免费观看| 欧美精品一区二区免费开放| 中文字幕免费在线视频6| 午夜福利在线观看免费完整高清在| 一级毛片我不卡| av在线蜜桃| 日韩免费高清中文字幕av| 永久免费av网站大全| 在现免费观看毛片| 你懂的网址亚洲精品在线观看| 男人爽女人下面视频在线观看| 国产成人freesex在线| 青青草视频在线视频观看| 亚洲不卡免费看| 日韩人妻高清精品专区| 亚洲伊人久久精品综合| 插阴视频在线观看视频| 久久精品国产亚洲av涩爱| 日韩中文字幕视频在线看片 | av专区在线播放| 亚洲伊人久久精品综合| 麻豆乱淫一区二区| 免费大片黄手机在线观看| 亚洲精品乱久久久久久| 在线观看免费视频网站a站| 亚洲熟女精品中文字幕| 欧美日韩国产mv在线观看视频 | 男女免费视频国产| 亚洲婷婷狠狠爱综合网| 国产色婷婷99| 亚洲熟女精品中文字幕| 一级黄片播放器| 日本-黄色视频高清免费观看| 久久精品久久精品一区二区三区| 欧美区成人在线视频| 免费大片18禁| 免费av不卡在线播放| 国产成人a∨麻豆精品| 国产有黄有色有爽视频| 日本一二三区视频观看| 不卡视频在线观看欧美| 日日摸夜夜添夜夜爱| 午夜免费男女啪啪视频观看| 久久久久人妻精品一区果冻| 国产爽快片一区二区三区| 搡女人真爽免费视频火全软件| 国产精品蜜桃在线观看| 99久久中文字幕三级久久日本| 极品少妇高潮喷水抽搐| 精品久久国产蜜桃| 一级毛片aaaaaa免费看小| 天堂8中文在线网| 在线播放无遮挡| 七月丁香在线播放| 久久久精品免费免费高清| 色哟哟·www| 不卡视频在线观看欧美| 啦啦啦中文免费视频观看日本| 国产精品久久久久久久电影| 一级毛片我不卡| av天堂中文字幕网| 乱系列少妇在线播放| 精品人妻偷拍中文字幕| 国产欧美日韩一区二区三区在线 | 精品酒店卫生间| 一区二区三区四区激情视频| 久久精品国产亚洲av涩爱| 性色avwww在线观看| 在线天堂最新版资源| 一本一本综合久久| 永久免费av网站大全| 美女福利国产在线 | 色网站视频免费| 男人爽女人下面视频在线观看| 高清黄色对白视频在线免费看 | 日韩视频在线欧美| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品aⅴ在线观看| 国产中年淑女户外野战色| 又大又黄又爽视频免费| 成人一区二区视频在线观看| 久久精品久久精品一区二区三区| 亚洲精华国产精华液的使用体验| 99久久综合免费| 国产乱人视频| 亚洲精品乱码久久久v下载方式| 亚洲内射少妇av| 99久久精品一区二区三区| 亚洲精品色激情综合| 欧美 日韩 精品 国产| 国产精品秋霞免费鲁丝片| 精品亚洲乱码少妇综合久久| 一级毛片电影观看| 亚洲人成网站高清观看| 视频中文字幕在线观看| h视频一区二区三区| av国产精品久久久久影院| 看十八女毛片水多多多| 亚洲人与动物交配视频| 中文字幕av成人在线电影| 国产欧美日韩一区二区三区在线 | 人妻一区二区av| 日韩成人av中文字幕在线观看| 国产伦在线观看视频一区| 2018国产大陆天天弄谢| 久热这里只有精品99| 男男h啪啪无遮挡| av福利片在线观看| 国产中年淑女户外野战色| 只有这里有精品99| 免费黄频网站在线观看国产| 九色成人免费人妻av| 精品酒店卫生间| 成人无遮挡网站| 亚洲精品乱码久久久久久按摩| 免费在线观看成人毛片| 久久久久精品久久久久真实原创| 精品人妻一区二区三区麻豆| 99九九线精品视频在线观看视频| 国产精品99久久久久久久久| 日韩亚洲欧美综合| 久久精品国产自在天天线| 尤物成人国产欧美一区二区三区| 国产高清国产精品国产三级 | 国精品久久久久久国模美| 欧美激情国产日韩精品一区| 老司机影院成人| 少妇被粗大猛烈的视频| av黄色大香蕉| 欧美一区二区亚洲| 国产精品.久久久| 免费黄色在线免费观看| 日韩制服骚丝袜av| 国产成人精品一,二区| 免费观看性生交大片5| 网址你懂的国产日韩在线| 亚洲精品成人av观看孕妇| 久久人妻熟女aⅴ| 欧美 日韩 精品 国产| 又粗又硬又长又爽又黄的视频| 久久久欧美国产精品| 插阴视频在线观看视频| 亚洲精品aⅴ在线观看| 欧美成人午夜免费资源| 舔av片在线| 国产欧美亚洲国产| 最近的中文字幕免费完整| 久久99精品国语久久久| 成人美女网站在线观看视频| xxx大片免费视频| 人妻 亚洲 视频| 能在线免费看毛片的网站| 亚洲精品aⅴ在线观看| 久久精品国产亚洲av天美| 久久人人爽人人爽人人片va| 亚洲精华国产精华液的使用体验| 国产亚洲最大av| 免费观看av网站的网址| 在线亚洲精品国产二区图片欧美 | 国产成人freesex在线| av国产久精品久网站免费入址| 男女免费视频国产| 高清毛片免费看| 亚洲av.av天堂| 另类亚洲欧美激情| 高清欧美精品videossex| 日产精品乱码卡一卡2卡三| 国产久久久一区二区三区| 亚洲av电影在线观看一区二区三区| 在线观看免费视频网站a站| 欧美成人一区二区免费高清观看| 国产精品福利在线免费观看| 亚洲图色成人| 久久综合国产亚洲精品| 国产白丝娇喘喷水9色精品| 建设人人有责人人尽责人人享有的 | 欧美精品人与动牲交sv欧美| 免费看av在线观看网站| 七月丁香在线播放| 国产成人午夜福利电影在线观看| av卡一久久| 日韩人妻高清精品专区| av在线老鸭窝| 草草在线视频免费看| 日日摸夜夜添夜夜添av毛片| 大香蕉97超碰在线|