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

    基于瞬態(tài)分析的柴油機(jī)活塞疲勞壽命預(yù)測

    2017-11-08 11:32:42何盼攀劉建敏王普凱劉艷斌康琦
    車用發(fā)動(dòng)機(jī) 2017年5期
    關(guān)鍵詞:傳熱系數(shù)瞬態(tài)活塞

    何盼攀,劉建敏,王普凱,劉艷斌,康琦

    (1.裝甲兵工程學(xué)院機(jī)械工程系,北京 100072; 2.裝甲兵工程學(xué)院訓(xùn)練部,北京 100072)

    基于瞬態(tài)分析的柴油機(jī)活塞疲勞壽命預(yù)測

    何盼攀1,劉建敏2,王普凱1,劉艷斌1,康琦1

    (1.裝甲兵工程學(xué)院機(jī)械工程系,北京 100072; 2.裝甲兵工程學(xué)院訓(xùn)練部,北京 100072)

    以某增壓柴油機(jī)活塞為研究對(duì)象,建立了由曲柄連桿機(jī)構(gòu)和缸套組成的裝配體有限元模型,計(jì)算了活塞在熱載荷、機(jī)械載荷和熱-機(jī)耦合作用下的應(yīng)力分布,在此基礎(chǔ)上將計(jì)算得到的熱-機(jī)耦合應(yīng)力場作為疲勞載荷,采用名義應(yīng)力法對(duì)活塞進(jìn)行疲勞壽命計(jì)算。結(jié)果表明:活塞的短壽命區(qū)域出現(xiàn)在活塞銷座內(nèi)側(cè)上部,最低循環(huán)次數(shù)為8.823×107次,折合1 470.5 h。從計(jì)算結(jié)果看,活塞的結(jié)構(gòu)較為合理,能滿足柴油機(jī)的使用要求。

    柴油機(jī);活塞;瞬態(tài)分析;熱-機(jī)耦合;疲勞壽命

    隨著柴油機(jī)強(qiáng)化程度的不斷提高,柴油機(jī)各部件承受的負(fù)荷不斷增加,特別是對(duì)于活塞而言,不僅受到瞬變的高溫燃?xì)庾饔?,而且承受著沖擊性的高頻機(jī)械載荷,在由此產(chǎn)生的熱-機(jī)耦合應(yīng)力的反復(fù)作用下,活塞容易出現(xiàn)疲勞失效。因此,對(duì)活塞進(jìn)行疲勞壽命預(yù)測對(duì)于活塞結(jié)構(gòu)的改進(jìn)和故障的預(yù)判具有重要的價(jià)值[1]。

    疲勞壽命計(jì)算往往建立在有限元仿真的基礎(chǔ)之上,以往由于計(jì)算機(jī)硬件的限制,活塞應(yīng)力場主要以靜態(tài)分析為主,而靜態(tài)分析無法考察活塞在動(dòng)態(tài)載荷作用下的應(yīng)力變化,這就直接影響了活塞疲勞壽命預(yù)測的準(zhǔn)確性[2]。因此,本研究對(duì)某柴油機(jī)活塞進(jìn)行了瞬態(tài)分析,計(jì)算了活塞在標(biāo)定工況下(2 000 r/min)的應(yīng)力時(shí)間歷程,并求得了在此工況下活塞的疲勞壽命。

    1 計(jì)算模型

    研究對(duì)象為V型12缸增壓柴油機(jī),應(yīng)用GT-Power建立了其在標(biāo)定工況下的工作過程計(jì)算模型,計(jì)算得到柴油機(jī)的各項(xiàng)參數(shù)。

    在運(yùn)行過程中,活塞不僅受到周期性變化的高溫燃?xì)庾饔?,而且還受到活塞銷的反作用力、交變的側(cè)壓力和往復(fù)慣性力等,這些作用力的大小和作用位置難以確定,在進(jìn)行有限元分析時(shí)通常是基于假設(shè)條件對(duì)這些作用力進(jìn)行簡化,因此單獨(dú)對(duì)活塞進(jìn)行有限元分析會(huì)造成一定的誤差,其分析過程也較為復(fù)雜[3]。因此,本研究建立了包括曲柄連桿機(jī)構(gòu)和氣缸套在內(nèi)的裝配體模型,將活塞復(fù)雜的邊界條件轉(zhuǎn)化為部件內(nèi)部之間的相互作用,重點(diǎn)對(duì)活塞進(jìn)行溫度場和應(yīng)力場的計(jì)算。

    用Creo建立各部件的的三維模型,完成裝配后將三維模型導(dǎo)入Ansys Workbench中,用四面體網(wǎng)格對(duì)活塞進(jìn)行網(wǎng)格劃分,設(shè)定最大網(wǎng)格尺寸為6 mm,最小網(wǎng)格尺寸為2 mm,并對(duì)活塞頂面、活塞頂?shù)酌?、活塞銷座、裙部等處進(jìn)行加密處理,其余部件采用自動(dòng)網(wǎng)格劃分,并對(duì)部件與部件之間的連接處進(jìn)行細(xì)化。裝配體和活塞網(wǎng)格模型見圖1。網(wǎng)格模型的單元數(shù)為172 442個(gè),節(jié)點(diǎn)數(shù)為297 783個(gè)?;钊牟牧蠟殄戜X合金,其具體屬性見表1。可以看出,材料對(duì)于溫度的變化比較敏感,當(dāng)溫度達(dá)到300 ℃時(shí),材料的拉伸極限強(qiáng)度已經(jīng)不足常溫下的一半,這對(duì)活塞的疲勞壽命影響很大。

    圖1 裝配體和活塞網(wǎng)格模型

    溫度/℃導(dǎo)熱系數(shù)/W·(m·K)-1熱膨脹系數(shù)/K-1彈性模量/GPa拉伸極限強(qiáng)度/MPa20146.321.8×10-672390100150.522.6×10-668380200158.823.9×10-662325300167.224.3×10-650165400171.424.9×10-64275

    2 載荷邊界條件的確定

    2.1熱載荷邊界條件

    活塞的溫度分布采用第三類邊界條件[4]計(jì)算,即

    Φ=hA(tw-tf)。

    (1)

    式中:Φ為對(duì)流傳熱量;h為介質(zhì)的對(duì)流傳熱系數(shù);tw為活塞與周圍流場接觸面的溫度;tf為活塞周圍流場溫度。

    活塞內(nèi)部的導(dǎo)熱為非穩(wěn)態(tài)導(dǎo)熱,其導(dǎo)熱方程為

    (2)

    式中:t為溫度;τ為時(shí)間;k為活塞材料的導(dǎo)熱系數(shù);ρ為活塞材料密度;c為活塞材料比熱容;x,y,z為笛卡爾坐標(biāo)分量。

    活塞各部位的傳熱邊界條件由試驗(yàn)和經(jīng)驗(yàn)公式相結(jié)合的方式確定。

    2.1.1活塞頂面

    (3)

    (4)

    式中:Tg為燃?xì)馑矔r(shí)溫度;αg為燃?xì)馑矔r(shí)對(duì)流傳熱系數(shù);φ為曲軸轉(zhuǎn)角;φ0為終了時(shí)刻所對(duì)應(yīng)的曲軸轉(zhuǎn)角。

    圖2 瞬時(shí)溫度和瞬時(shí)對(duì)流傳熱系數(shù)

    2.1.2活塞側(cè)面

    活塞側(cè)面的換熱情況非常復(fù)雜。燃?xì)獾臒崃客ㄟ^活塞頂面?zhèn)髦粱鹆Π丁⒒钊h(huán)和裙部,再通過油膜、氣缸壁傳至冷卻水,這個(gè)過程包含了部件與不同介質(zhì)之間的對(duì)流傳熱,也包含了部件之間的導(dǎo)熱,分析難度較大,故此處將其近似處理成多層平板對(duì)流傳熱[6],總的傳熱系數(shù)由各部件的導(dǎo)熱系數(shù)和流體的對(duì)流傳熱系數(shù)共同決定。其中火力岸處傳熱系數(shù)k1和環(huán)區(qū)的傳熱系數(shù)k2的計(jì)算公式分別為

    (5)

    (6)

    式中:a1為火力岸與缸套的間隙;a2為環(huán)區(qū)上沿的間隙;b為缸套厚度;c為環(huán)中心間距;λ1為燃?xì)獾膶?dǎo)熱系數(shù);λ2為缸套的導(dǎo)熱系數(shù);λ3為活塞環(huán)的導(dǎo)熱系數(shù);hω為缸套與冷卻水之間的對(duì)流傳熱系數(shù)。

    2.1.3活塞內(nèi)腔

    本研究的柴油機(jī)活塞采用機(jī)油飛濺潤滑的方式對(duì)內(nèi)腔進(jìn)行冷卻,其傳熱系數(shù)hoil的計(jì)算公式為

    (7)

    式中:t1為缸內(nèi)燃?xì)鉁囟?;t1為活塞頂面溫度;toil為活塞頂下內(nèi)側(cè)壁面的溫度;δ為活塞頂厚度。

    采用硬度塞法對(duì)活塞表面的溫度進(jìn)行測量,測點(diǎn)布置見圖3。將測量值和仿真值進(jìn)行比較,并對(duì)邊界條件進(jìn)行修正,最后確定穩(wěn)態(tài)下活塞各部位的傳熱系數(shù)和對(duì)應(yīng)的環(huán)境溫度(見表2),各測點(diǎn)的測量值和仿真值對(duì)比見表3。

    圖3 活塞測點(diǎn)布置

    區(qū)域溫度/℃?zhèn)鳠嵯禂?shù)/W·(m2·K)-1活塞頂820615.7火力岸90350第一道環(huán)槽901920第二道環(huán)槽902000第三道環(huán)槽901900裙部90390內(nèi)腔100650

    表3 活塞測點(diǎn)溫度對(duì)比結(jié)果

    從對(duì)比結(jié)果看,除了測點(diǎn)D與測點(diǎn)E誤差略大外,其余各點(diǎn)均在合理誤差范圍內(nèi)。

    2.2機(jī)械載荷邊界條件

    裝配體受到的作用力為缸內(nèi)燃?xì)鈮毫颓S的扭矩,外部邊界條件比較簡單,但是其內(nèi)部則包含了多個(gè)連接關(guān)系,由于研究的重點(diǎn)為活塞的應(yīng)力分布,并不關(guān)注其他部件連接處的應(yīng)力情況,所以對(duì)于遠(yuǎn)離活塞的連接關(guān)系可以用運(yùn)動(dòng)副來代替接觸,這樣可以在保證精度的前提下節(jié)約計(jì)算資源[7]。同時(shí),在柴油機(jī)實(shí)際運(yùn)行過程中,部件之間的連接處會(huì)產(chǎn)生摩擦力,而摩擦力的出現(xiàn)會(huì)導(dǎo)致仿真計(jì)算收斂困難,因此,在設(shè)定邊界條件時(shí)忽略了摩擦的影響。各部件之間連接關(guān)系見表4。

    表4 裝配體各部件之間的連接關(guān)系

    3 瞬態(tài)溫度場的計(jì)算與分析

    在進(jìn)行瞬態(tài)溫度場分析時(shí),首先計(jì)算活塞的穩(wěn)態(tài)溫度場(計(jì)算結(jié)果見圖4),從而得到瞬態(tài)分析的初始條件,然后采用自動(dòng)時(shí)間步長的方法,設(shè)置初始步長為100個(gè)載荷步,最大步長為200個(gè)載荷步,將一個(gè)周期內(nèi)的瞬時(shí)對(duì)流傳熱系數(shù)和瞬時(shí)溫度施加在活塞頂面上,其他部位的邊界條件在穩(wěn)態(tài)分析的基礎(chǔ)上保持不變,由此計(jì)算得到活塞的瞬態(tài)溫度場。

    圖4 活塞穩(wěn)態(tài)溫度場

    計(jì)算完畢后,分別選取了活塞頂不同深度不同部位的溫度,以考察活塞在一個(gè)周期內(nèi)溫度的波動(dòng)情況,結(jié)果見圖5。可以看出,活塞頂面的溫度波動(dòng)最為劇烈,而隨著深度的增加,溫度的波動(dòng)逐漸變緩,當(dāng)深度超過2 mm后,溫度基本保持不變。

    圖5 活塞頂不同深度溫度時(shí)間歷程

    4 瞬態(tài)應(yīng)力場的計(jì)算與分析

    4.1機(jī)械載荷作用下的瞬態(tài)應(yīng)力場

    缸內(nèi)燃?xì)鈮毫ψ兓€由GT-Power計(jì)算得到,計(jì)算結(jié)果見圖6。將裝配體模型導(dǎo)入瞬態(tài)動(dòng)力學(xué)模塊中,在活塞頂面以載荷步的形式施加燃?xì)鈮毫?載荷步的選取與前文保持一致),同時(shí)在曲軸兩端施加2 000 r/min的恒定轉(zhuǎn)速,并約束缸套上表面沿x,y方向和缸套下定位帶表面沿z方向的移動(dòng)和轉(zhuǎn)動(dòng)。

    圖6 缸內(nèi)燃?xì)鈮毫ψ兓€

    用第四強(qiáng)度理論計(jì)算活塞受到的等效應(yīng)力,即Von-Mises應(yīng)力,計(jì)算結(jié)果見圖7與圖8。圖7示出了最高燃燒壓力時(shí)刻(367.2°)活塞機(jī)械應(yīng)力分布云圖,圖8示出了一個(gè)周期內(nèi)活塞最大機(jī)械應(yīng)力時(shí)間歷程。可以看出:

    1) 在最高燃燒壓力時(shí)刻,活塞銷座處的應(yīng)力高達(dá)171.3 MPa,是活塞結(jié)構(gòu)的危險(xiǎn)點(diǎn);而活塞頂面和活塞裙部由于分別受到燃?xì)獾臎_擊和側(cè)推力的作用,導(dǎo)致應(yīng)力值都在較高的水平;

    2) 與圖6缸壓曲線比較發(fā)現(xiàn),活塞最大應(yīng)力的變化與燃?xì)鈮毫Φ淖兓瘞缀跬?,燃?xì)獾牟▌?dòng)直接影響活塞應(yīng)力的大小,特別是在最高燃燒壓力附近,應(yīng)力變化十分顯著。

    圖7 最高燃燒壓力時(shí)刻活塞機(jī)械應(yīng)力分布云圖

    圖8 活塞最大機(jī)械應(yīng)力時(shí)間歷程

    4.2熱載荷作用下的瞬態(tài)應(yīng)力場

    將前文分析得到的瞬態(tài)溫度結(jié)果以體載荷的形式施加到活塞上,并施加相同的恒定轉(zhuǎn)速和位移約束條件,由此求得了瞬態(tài)熱應(yīng)力場。計(jì)算結(jié)果見圖9與圖10。圖9示出了最高燃燒壓力時(shí)刻活塞應(yīng)力分布云圖,圖10示出了活塞最大應(yīng)力和活塞頂面中心的應(yīng)力時(shí)間歷程。結(jié)合圖9與圖10可以看出:

    1) 活塞最大應(yīng)力出現(xiàn)在活塞頂?shù)南碌酌?,且?yīng)力值在一個(gè)周期內(nèi)基本保持不變;

    2) 活塞頂面的應(yīng)力相對(duì)較小,但波動(dòng)幅度達(dá)到9.242 MPa,這說明活塞頂面應(yīng)力受其溫度場的變化影響較大。

    圖9 最高燃燒壓力時(shí)刻活塞熱應(yīng)力分布云圖

    圖10 活塞不同部位熱應(yīng)力時(shí)間歷程

    4.3熱-機(jī)耦合作用下的瞬態(tài)應(yīng)力場

    將燃?xì)鈮毫退矐B(tài)溫度計(jì)算結(jié)果同時(shí)施加于活塞上,并施加相同的約束條件,由此計(jì)算得到活塞熱-機(jī)耦合應(yīng)力場,計(jì)算結(jié)果見圖11與圖12。圖11示出了活塞在最高燃燒壓力時(shí)刻和最小壓力時(shí)刻的熱-機(jī)耦合應(yīng)力分布云圖,圖12示出了活塞最大應(yīng)力時(shí)間歷程。結(jié)合兩圖可以看出:

    1) 在熱-機(jī)耦合作用下,不同時(shí)刻活塞的應(yīng)力分布有較大的區(qū)別,其最大應(yīng)力較單獨(dú)施加機(jī)械載荷時(shí)有所降低,熱載荷的作用在一定程度上減小了等效應(yīng)力值,而活塞總體的應(yīng)力值較單獨(dú)機(jī)械載荷作用時(shí)有所增加,環(huán)槽、活塞內(nèi)壁、活塞頂?shù)酌娴忍幍膽?yīng)力值都相應(yīng)增大;

    2) 通過與圖8、圖10的最大應(yīng)力曲線比較發(fā)現(xiàn):在熱-機(jī)耦合作用下,活塞的熱應(yīng)力在大部分時(shí)間內(nèi)占主導(dǎo)作用,但在最高燃燒壓力附近,機(jī)械載荷迅速占據(jù)主導(dǎo),這是因?yàn)樵跓釕T性的作用下,活塞總體溫度場分布基本不變,使得活塞的熱應(yīng)力基本保持不變,而燃?xì)鈮毫Φ淖兓苯佑绊懟钊麘?yīng)力的變化,不存在滯后性,所以當(dāng)燃?xì)鈮毫Ψ€(wěn)定在一個(gè)較低的水平時(shí),活塞的應(yīng)力分布與單獨(dú)作用熱載荷時(shí)的分布相似,而當(dāng)燃?xì)鈮毫眲∽兓瘯r(shí),應(yīng)力迅速響應(yīng),從而發(fā)生大幅度變化。

    圖11 活塞不同時(shí)刻熱-機(jī)耦合應(yīng)力分布云圖

    圖12 活塞最大熱-機(jī)耦合應(yīng)力時(shí)間歷程

    5 疲勞壽命預(yù)測分析

    5.1基本方法

    疲勞是指材料在循環(huán)應(yīng)力和應(yīng)變作用下,在一處或幾處逐漸產(chǎn)生局部永久性累積損傷,經(jīng)一定循環(huán)次數(shù)后產(chǎn)生裂紋或突然發(fā)生完全斷裂的過程[8]。大量試驗(yàn)研究表明,應(yīng)力幅是疲勞破壞的主要因素。

    對(duì)于常規(guī)疲勞設(shè)計(jì),一般可以分為以下5個(gè)步驟:1)根據(jù)零件可能出現(xiàn)的疲勞破壞形式確定疲勞分析方法;2)對(duì)零件進(jìn)行有限元分析,得到零件的應(yīng)力分布;3)獲取材料的疲勞壽命曲線;4)獲取零件的載荷時(shí)間歷程,并處理成載荷譜;5)選擇合適的疲勞累積損傷準(zhǔn)則,結(jié)合材料的疲勞壽命曲線進(jìn)行疲勞計(jì)算。

    5.2仿真計(jì)算

    根據(jù)熱-機(jī)耦合仿真結(jié)果,活塞受到的最大應(yīng)力為163.05 MPa,小于活塞材料的屈服極限,結(jié)合柴油機(jī)的實(shí)際運(yùn)行情況,可以判斷出活塞熱-機(jī)耦合作用下的疲勞失效形式為高周疲勞,故采用名義應(yīng)力法計(jì)算。

    5.2.1活塞材料的S-N曲線

    活塞材料的拉伸極限強(qiáng)度隨溫度變化,不同溫度下的S-N曲線見圖13,N為循環(huán)次數(shù)。

    圖13 活塞材料S-N曲線

    5.2.2載荷譜

    忽略柴油機(jī)在工作過程中的工況變化以及外載荷的隨機(jī)波動(dòng),因此活塞在每個(gè)循環(huán)下的應(yīng)力變化都是保持一致的,故可直接將一個(gè)循環(huán)的熱-機(jī)耦合應(yīng)力時(shí)間歷程作為載荷譜。

    5.2.3疲勞強(qiáng)度修正

    材料的S-N曲線實(shí)際上是以標(biāo)準(zhǔn)光滑試樣為對(duì)象在對(duì)稱循環(huán)載荷的作用下測得的,而實(shí)際情況下零件的形狀、尺寸、表面情況、受到的載荷情況等都會(huì)影響零件的疲勞強(qiáng)度,需要對(duì)這些因素進(jìn)行討論:

    1) 尺寸效應(yīng)的影響

    一般而言,零件的尺寸越大,疲勞強(qiáng)度越低,這種現(xiàn)象稱作尺寸效應(yīng),可用尺寸系數(shù)ε來表征,它代表幾何尺寸相似的零件與標(biāo)準(zhǔn)試樣的疲勞極限的比值。尺寸系數(shù)可通過查找尺寸系數(shù)曲線確認(rèn),活塞的直徑為150 mm,通過查圖可知,尺寸系數(shù)近似取0.75。

    2) 表面粗糙度的影響

    零件表面加工粗糙,相當(dāng)于存在許多微缺口,零件受到載荷時(shí)便會(huì)產(chǎn)生應(yīng)力集中,影響疲勞強(qiáng)度。為了考慮表面粗糙度對(duì)疲勞強(qiáng)度的影響,在計(jì)算過程中引入了表面加工系數(shù)β1,表面加工系數(shù)與材料拉伸極限強(qiáng)度和表面加工方式之間存在一定的關(guān)系,β1代表某種加工表面的標(biāo)準(zhǔn)試樣和標(biāo)準(zhǔn)光滑試樣的疲勞極限的比值,由于活塞頂面在加工過程中進(jìn)行了拋光處理,較為光滑,因此可近似取β1=1。

    3) 平均應(yīng)力的影響

    活塞受到的是非對(duì)稱循環(huán)載荷,平均應(yīng)力不為零,因此必須進(jìn)行平均應(yīng)力修正[9]。其修正后的等效應(yīng)力幅σA可以由下式[10]確定:

    σA=σa+ψσm。

    (8)

    式中:σa為應(yīng)力幅;ψ為平均應(yīng)力影響系數(shù),鋁合金ψ取0.335;σm為平均應(yīng)力。

    5.2.4疲勞累積損傷準(zhǔn)則

    線性累積損傷準(zhǔn)則(即Miner準(zhǔn)則)認(rèn)為,當(dāng)構(gòu)件吸收的能量達(dá)到極限值時(shí)即產(chǎn)生疲勞破壞,此時(shí)其總損傷值為1,構(gòu)件吸收的能量與不同應(yīng)力水平所對(duì)應(yīng)的循環(huán)數(shù)之間呈正比。在實(shí)際運(yùn)用中,通常對(duì)Miner準(zhǔn)則進(jìn)行修正[11],結(jié)合零部件實(shí)際情況,確認(rèn)總損傷值a,經(jīng)過修正后的Miner準(zhǔn)則表達(dá)式如下:

    (9)

    式中:D為損傷值;l為應(yīng)力水平;ni為各應(yīng)力水平下的循環(huán)次數(shù);Ni為各應(yīng)力水平下的疲勞壽命;a取0.7。

    5.3結(jié)果分析

    在Ncode designlife中構(gòu)建疲勞分析框圖,將活塞熱-機(jī)耦合瞬態(tài)分析結(jié)果導(dǎo)入輸入模塊中,調(diào)整分析引擎中的載荷譜來源為時(shí)間步輸入,選擇分析方法為MultiTemperatureCurve,并添加材料映射,最后選擇線性累積損傷準(zhǔn)則進(jìn)行疲勞計(jì)算。

    圖14和圖15分別示出了活塞的疲勞壽命云圖和活塞關(guān)鍵點(diǎn)的應(yīng)力時(shí)間歷程。從圖中可以看出:

    1) 活塞最小循環(huán)次數(shù)為8.823×107次,折合柴油機(jī)運(yùn)行時(shí)間為1 470.5 h,大于該柴油機(jī)的大修間隔期(1 000 h),活塞可以滿足使用要求;

    2) 活塞疲勞危險(xiǎn)位置為活塞銷座內(nèi)側(cè)上部,觀察該處的應(yīng)力時(shí)間歷程可以發(fā)現(xiàn),該處不僅出現(xiàn)了應(yīng)力集中,是活塞的結(jié)構(gòu)危險(xiǎn)點(diǎn),而且由于一個(gè)周期內(nèi)熱應(yīng)力和機(jī)械應(yīng)力的循環(huán)交替,其應(yīng)力幅也很大,導(dǎo)致該處極易出現(xiàn)疲勞失效;同時(shí)活塞頂面處的疲勞壽命也較短,這是因?yàn)樵撎帨囟容^高,導(dǎo)致疲勞強(qiáng)度嚴(yán)重下降,而且由于表層溫度的周期性變化和燃?xì)鈮毫Φ闹苯幼饔?,該處的?yīng)力波動(dòng)也較為顯著;

    3) 活塞頂面的應(yīng)力水平小于活塞頂?shù)酌?,而其疲勞壽命卻遠(yuǎn)低于活塞頂?shù)酌娴钠趬勖?,這是因?yàn)榛钊斆嬷苯映惺芨邷馗邏旱娜細(xì)庾饔?,?yīng)力變化較為劇烈,應(yīng)力幅較大(14.71 MPa),對(duì)疲勞壽命影響大;而活塞底面雖然應(yīng)力值較大,但應(yīng)力幅較小(6.685 MPa),相對(duì)而言受到的載荷波動(dòng)較為平緩,而且受到的熱負(fù)荷影響也相對(duì)較小,因此疲勞壽命較高。

    圖14 活塞疲勞壽命云圖

    圖15 活塞不同部位熱-機(jī)耦合應(yīng)力時(shí)間歷程

    6 結(jié)束語

    建立了由曲柄連桿機(jī)構(gòu)和缸套等部件組成的裝配體模型,對(duì)活塞進(jìn)行了瞬態(tài)分析,分別研究了活塞在熱載荷作用、機(jī)械載荷作用和熱-機(jī)耦合作用下的應(yīng)力場分布情況,求得了一個(gè)循環(huán)中活塞的應(yīng)力時(shí)間歷程。

    在瞬態(tài)分析的基礎(chǔ)上,對(duì)應(yīng)力時(shí)間歷程進(jìn)行了平均應(yīng)力修正,并考慮了活塞的尺寸效應(yīng)和表面粗糙度的影響,然后采用名義應(yīng)力法計(jì)算了活塞的疲勞壽命,分析得出活塞危險(xiǎn)位置為銷座內(nèi)側(cè)上部,該處的循環(huán)次數(shù)為8.823×107次,能滿足柴油機(jī)的使用要求,驗(yàn)證了活塞結(jié)構(gòu)的合理性。

    從應(yīng)力時(shí)間歷程入手,分析了活塞不同位置的壽命情況,進(jìn)一步驗(yàn)證了應(yīng)力幅是影響疲勞壽命的主要因素,同時(shí)溫度對(duì)活塞疲勞壽命影響也較大。

    [1] 張衛(wèi)正,劉金祥,原彥鵬,等.內(nèi)燃機(jī)失效分析與評(píng)估[M].北京:北京航空航天大學(xué)出版社,2011:150-153.

    [2] 田城.YC6MK375N-30型CNG發(fā)動(dòng)機(jī)活塞熱-機(jī)耦合分析與疲勞壽命預(yù)測[D].武漢:武漢理工大學(xué),2013:52-53.

    [3] 王文竹,程勉宏,劉剛,等.發(fā)動(dòng)機(jī)活塞組的耦合有限元分析[J].機(jī)械設(shè)計(jì)與制造,2009(4):186-188.

    [4] 姚仲鵬,王新國.車輛冷卻傳熱[M].北京:北京理工大學(xué)出版社,2001:17-25.

    [5] 楊世銘,陶文銓.傳熱學(xué)[M].北京:高等教育出版社,1998.

    [6] 樓狄明,張志穎,王禮麗.機(jī)車柴油機(jī)組合活塞的換熱邊界條件及熱負(fù)荷[J].同濟(jì)大學(xué)學(xué)報(bào),2005,33(5):664-667.

    [7] 宋少云,尹芳.裝配體仿真中用運(yùn)動(dòng)副代替接觸的有效性分析[J].機(jī)械設(shè)計(jì)與制造,2013(6):206-208.

    [8] 趙少汴.抗疲勞設(shè)計(jì)手冊(cè)[M].北京:機(jī)械工業(yè)出版社,2015.

    [9] 吳波,魏志明,王紅麗,等.柴油機(jī)連桿疲勞強(qiáng)度評(píng)價(jià)研究[J].車用發(fā)動(dòng)機(jī),2014(2):6-10.

    [10] 聞邦椿.機(jī)械設(shè)計(jì)手冊(cè)[M].北京:機(jī)械工業(yè)出版社,2014:54.

    [11] 胡定云,陳澤忠,溫世杰,等.某柴油機(jī)氣缸蓋疲勞的可靠性預(yù)測[J].車用發(fā)動(dòng)機(jī),2008(增刊):38-44.

    PredictionofFatigueLifeforDieselEnginePistonBasedonTransientAnalysis

    HE Panpan1,LIU Jianmin2,WANG Pukai1,LIU Yanbin1,KANG Qi1

    (1.Department of Mechanical Engineering,Academy of Armored Force Engineering,Beijing 100072,China;2.Department of Training,Academy of Armored Force Engineering,Beijing 100072,China)

    Taking the piston of turbocharged diesel engine as the research object, the FEA model of crank and connecting rod mechanism and cylinder assembly was built and the stress distribution of piston was calculated under thermal load,mechanical load and thermal and mechanical coupling load. Then the thermal and mechanical coupled stress field was used as fatigue load and the fatigue life of piston was computed according to nominal stress approach. The results show that the short life region was the upper part inside piston pin boss. In this area, the minimum cycle is 8.823×107or around 1 470.5 h.Accordingly, the structure of piston is reasonable and can meet the application requirements of diesel engine.

    diesel engine;piston;transient analysis;thermal-mechanical couple;fatigue life

    2017-04-12;

    2017-07-11

    國家“973”計(jì)劃項(xiàng)目(坦克裝甲車輛柴油機(jī)高原性能重大基礎(chǔ)研究,項(xiàng)目編號(hào)201697301)

    何盼攀(1993—),男,碩士,主要研究方向?yàn)檠b甲車輛動(dòng)力系統(tǒng)總體技術(shù);413901266@qq.com。

    10.3969/j.issn.1001-2222.2017.05.011

    TK423.33

    B

    1001-2222(2017)05-0057-07

    [編輯: 姜曉博]

    猜你喜歡
    傳熱系數(shù)瞬態(tài)活塞
    一種活塞擠壓式室內(nèi)通風(fēng)設(shè)備的研發(fā)與應(yīng)用
    探析寒冷地區(qū)75%建筑節(jié)能框架下圍護(hù)結(jié)構(gòu)熱工性能的重組
    高壓感應(yīng)電動(dòng)機(jī)斷電重啟時(shí)的瞬態(tài)仿真
    新型鋁合金節(jié)能窗傳熱系數(shù)和簡化計(jì)算
    十億像素瞬態(tài)成像系統(tǒng)實(shí)時(shí)圖像拼接
    基于瞬態(tài)流場計(jì)算的滑動(dòng)軸承靜平衡位置求解
    聚乳酸吹膜過程中傳熱系數(shù)的研究
    中國塑料(2015年2期)2015-10-14 05:34:24
    DC/DC變換器中的瞬態(tài)特性分析
    KS Kolbenschmidt公司的新型鋼活塞
    低噪聲活塞——降低活塞銷的拍擊噪聲
    51国产日韩欧美| 亚洲性久久影院| 国产白丝娇喘喷水9色精品| 国产精品久久久久久亚洲av鲁大| 九色成人免费人妻av| 欧美在线一区亚洲| 免费黄网站久久成人精品| 亚洲中文日韩欧美视频| 一区二区三区高清视频在线| 亚洲一级一片aⅴ在线观看| 女人十人毛片免费观看3o分钟| 免费在线观看成人毛片| 91在线观看av| 午夜福利高清视频| 国产国拍精品亚洲av在线观看| 日韩欧美精品免费久久| 一级黄色大片毛片| 在线观看av片永久免费下载| 国产av不卡久久| 国产成人a∨麻豆精品| 老司机午夜福利在线观看视频| 看片在线看免费视频| 亚洲欧美日韩无卡精品| 国产伦精品一区二区三区视频9| 蜜桃亚洲精品一区二区三区| 国产成人a区在线观看| 别揉我奶头 嗯啊视频| 天堂av国产一区二区熟女人妻| 赤兔流量卡办理| 久久久欧美国产精品| 又爽又黄a免费视频| 校园春色视频在线观看| 非洲黑人性xxxx精品又粗又长| 性插视频无遮挡在线免费观看| 国产高清有码在线观看视频| 麻豆国产97在线/欧美| 哪里可以看免费的av片| 日韩欧美免费精品| а√天堂www在线а√下载| 91久久精品国产一区二区成人| 嫩草影院入口| 亚洲国产精品成人久久小说 | 丰满乱子伦码专区| 男女边吃奶边做爰视频| 少妇猛男粗大的猛烈进出视频 | 我的女老师完整版在线观看| 亚洲内射少妇av| 国产黄色小视频在线观看| 性插视频无遮挡在线免费观看| 国产片特级美女逼逼视频| 亚洲一区高清亚洲精品| 观看免费一级毛片| 小蜜桃在线观看免费完整版高清| 大香蕉久久网| 亚洲国产高清在线一区二区三| 国产精品综合久久久久久久免费| 国产一区二区三区av在线 | 国产成年人精品一区二区| 国产精品精品国产色婷婷| 一本一本综合久久| 欧洲精品卡2卡3卡4卡5卡区| 日日摸夜夜添夜夜添小说| 成年版毛片免费区| 亚洲精品亚洲一区二区| avwww免费| 国产男人的电影天堂91| 偷拍熟女少妇极品色| 久久午夜福利片| 久久人人爽人人片av| 1024手机看黄色片| 欧美成人精品欧美一级黄| 老师上课跳d突然被开到最大视频| 看黄色毛片网站| 成年版毛片免费区| 又粗又爽又猛毛片免费看| 菩萨蛮人人尽说江南好唐韦庄 | 最近最新中文字幕大全电影3| 波多野结衣巨乳人妻| 亚洲成人av在线免费| 婷婷精品国产亚洲av在线| 99热网站在线观看| a级毛片a级免费在线| 男女视频在线观看网站免费| 亚洲人与动物交配视频| 日本三级黄在线观看| 欧美中文日本在线观看视频| 亚洲人成网站高清观看| 91狼人影院| av.在线天堂| 一区福利在线观看| 国产一区二区在线av高清观看| 99九九线精品视频在线观看视频| 亚洲婷婷狠狠爱综合网| 热99re8久久精品国产| 日本在线视频免费播放| 在线国产一区二区在线| 亚洲人成网站在线播放欧美日韩| 日本五十路高清| 色吧在线观看| av在线天堂中文字幕| а√天堂www在线а√下载| 亚洲国产精品成人久久小说 | 久久久a久久爽久久v久久| 人妻久久中文字幕网| 国产精品野战在线观看| 在线观看免费视频日本深夜| 国产真实伦视频高清在线观看| 亚洲精品色激情综合| 日本a在线网址| 亚洲精品国产av成人精品 | 男女下面进入的视频免费午夜| 久久精品国产亚洲av香蕉五月| a级毛片免费高清观看在线播放| 亚洲七黄色美女视频| 99在线视频只有这里精品首页| 亚洲欧美精品综合久久99| 99九九线精品视频在线观看视频| 色综合站精品国产| 最后的刺客免费高清国语| 韩国av在线不卡| 乱码一卡2卡4卡精品| 日本熟妇午夜| 中文字幕av成人在线电影| 91麻豆精品激情在线观看国产| 激情 狠狠 欧美| 简卡轻食公司| 亚洲图色成人| 两性午夜刺激爽爽歪歪视频在线观看| 日本-黄色视频高清免费观看| 精品久久久久久久久亚洲| 十八禁网站免费在线| 天天躁日日操中文字幕| 欧美xxxx性猛交bbbb| 老师上课跳d突然被开到最大视频| 91在线观看av| 日本黄色视频三级网站网址| 大型黄色视频在线免费观看| 少妇被粗大猛烈的视频| 色哟哟·www| 美女内射精品一级片tv| 99在线视频只有这里精品首页| 日本三级黄在线观看| 国产精品电影一区二区三区| 欧美最新免费一区二区三区| 亚洲乱码一区二区免费版| 国产高清激情床上av| 18禁在线播放成人免费| 最近在线观看免费完整版| 免费看av在线观看网站| 国产色婷婷99| 黄色日韩在线| avwww免费| 欧美一区二区亚洲| 欧美成人一区二区免费高清观看| 久久99热6这里只有精品| 波多野结衣巨乳人妻| 亚洲人成网站在线播放欧美日韩| 久久精品国产自在天天线| 久久久精品94久久精品| 国产成人影院久久av| 校园春色视频在线观看| 97超碰精品成人国产| 精品熟女少妇av免费看| 亚洲精品成人久久久久久| 日韩强制内射视频| 99久久精品国产国产毛片| 欧美丝袜亚洲另类| 欧美高清性xxxxhd video| 91在线观看av| 99久久久亚洲精品蜜臀av| 精华霜和精华液先用哪个| 淫秽高清视频在线观看| 一级毛片我不卡| 特大巨黑吊av在线直播| 久久人妻av系列| 黄色一级大片看看| 日本与韩国留学比较| 精品午夜福利在线看| 亚洲av一区综合| 18禁裸乳无遮挡免费网站照片| 秋霞在线观看毛片| 人妻丰满熟妇av一区二区三区| 免费av不卡在线播放| 18禁黄网站禁片免费观看直播| 噜噜噜噜噜久久久久久91| 1024手机看黄色片| 亚洲天堂国产精品一区在线| 伦精品一区二区三区| 日韩av不卡免费在线播放| 国产精品人妻久久久影院| 国模一区二区三区四区视频| 韩国av在线不卡| 免费一级毛片在线播放高清视频| 国产亚洲精品久久久久久毛片| 久久精品久久久久久噜噜老黄 | 人人妻人人澡人人爽人人夜夜 | 亚洲精品日韩av片在线观看| 亚洲国产欧美人成| 欧美色视频一区免费| av在线老鸭窝| 精品不卡国产一区二区三区| 少妇被粗大猛烈的视频| 俄罗斯特黄特色一大片| 最近2019中文字幕mv第一页| 天堂影院成人在线观看| 大香蕉久久网| 午夜福利18| АⅤ资源中文在线天堂| 少妇猛男粗大的猛烈进出视频 | 欧美xxxx性猛交bbbb| 国产v大片淫在线免费观看| 亚洲人成网站在线观看播放| av天堂中文字幕网| 国产精品福利在线免费观看| 色哟哟·www| 精品欧美国产一区二区三| 别揉我奶头 嗯啊视频| 久久久久免费精品人妻一区二区| 身体一侧抽搐| 亚洲av美国av| 日本免费a在线| 最近的中文字幕免费完整| 又黄又爽又刺激的免费视频.| 亚洲av二区三区四区| 午夜福利18| 有码 亚洲区| 国产欧美日韩一区二区精品| 晚上一个人看的免费电影| 免费av观看视频| 看免费成人av毛片| 欧美日韩在线观看h| 激情 狠狠 欧美| 在线免费十八禁| 国产精品伦人一区二区| 看片在线看免费视频| 久久久久久九九精品二区国产| www.色视频.com| 亚洲av熟女| 非洲黑人性xxxx精品又粗又长| 在线国产一区二区在线| 大型黄色视频在线免费观看| 久久精品综合一区二区三区| 人妻少妇偷人精品九色| 亚洲三级黄色毛片| ponron亚洲| 久久久精品94久久精品| 免费在线观看成人毛片| 国产精品不卡视频一区二区| 国产一区二区三区av在线 | 免费人成视频x8x8入口观看| 色av中文字幕| 国产探花极品一区二区| 成人三级黄色视频| 中文字幕熟女人妻在线| 欧美丝袜亚洲另类| 欧美+日韩+精品| 美女免费视频网站| 成人午夜高清在线视频| 欧美国产日韩亚洲一区| 欧美区成人在线视频| 国产av一区在线观看免费| 成年免费大片在线观看| 国产伦在线观看视频一区| 欧美不卡视频在线免费观看| 日韩av不卡免费在线播放| 我要看日韩黄色一级片| 一进一出抽搐动态| 真实男女啪啪啪动态图| 国产精品一区二区性色av| 国产一区二区激情短视频| 久久人人爽人人片av| 国产午夜精品论理片| 黄片wwwwww| 欧美xxxx黑人xx丫x性爽| 国产色爽女视频免费观看| 亚洲中文字幕一区二区三区有码在线看| 99热网站在线观看| 亚洲欧美清纯卡通| 大又大粗又爽又黄少妇毛片口| 又粗又爽又猛毛片免费看| 国产精品精品国产色婷婷| 中出人妻视频一区二区| 久久久久国产网址| 日韩制服骚丝袜av| 99在线视频只有这里精品首页| 人人妻人人看人人澡| 久久精品人妻少妇| 男女边吃奶边做爰视频| 国产精品一二三区在线看| 精品少妇黑人巨大在线播放 | 人妻丰满熟妇av一区二区三区| 12—13女人毛片做爰片一| 综合色av麻豆| 色视频www国产| 天天一区二区日本电影三级| 男女做爰动态图高潮gif福利片| 春色校园在线视频观看| 日韩大尺度精品在线看网址| 国产乱人视频| 国产精品三级大全| 干丝袜人妻中文字幕| 伦理电影大哥的女人| 国产成人91sexporn| 成人av在线播放网站| 97超碰精品成人国产| 亚洲美女黄片视频| 日韩成人av中文字幕在线观看 | 午夜福利在线在线| 天美传媒精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 全区人妻精品视频| 97在线视频观看| 久久人人爽人人片av| 国产精品国产高清国产av| 精华霜和精华液先用哪个| 一边摸一边抽搐一进一小说| 精品久久久噜噜| 国产三级中文精品| 国产高清有码在线观看视频| 老师上课跳d突然被开到最大视频| 悠悠久久av| 国产精品女同一区二区软件| 一级av片app| 欧美区成人在线视频| 精品久久久久久久久av| 少妇熟女欧美另类| 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 久99久视频精品免费| 内射极品少妇av片p| 丰满人妻一区二区三区视频av| 日韩 亚洲 欧美在线| 国产毛片a区久久久久| 在线免费观看的www视频| 最后的刺客免费高清国语| 午夜视频国产福利| 午夜福利18| 亚洲最大成人手机在线| 99久久精品一区二区三区| 少妇被粗大猛烈的视频| 欧美色视频一区免费| 日韩欧美在线乱码| 国内少妇人妻偷人精品xxx网站| 国国产精品蜜臀av免费| 国产91av在线免费观看| 精品一区二区三区视频在线观看免费| 久久久久久大精品| 成人亚洲精品av一区二区| 国产一级毛片七仙女欲春2| 亚洲性夜色夜夜综合| 内地一区二区视频在线| 久久精品91蜜桃| 黄色视频,在线免费观看| 国产精品永久免费网站| 日日啪夜夜撸| 欧美xxxx黑人xx丫x性爽| 免费av观看视频| avwww免费| 别揉我奶头~嗯~啊~动态视频| 日韩国内少妇激情av| 国产男人的电影天堂91| 老司机午夜福利在线观看视频| 色综合色国产| 最后的刺客免费高清国语| av专区在线播放| 91av网一区二区| 在现免费观看毛片| 欧美激情国产日韩精品一区| 久久99热这里只有精品18| 国产欧美日韩一区二区精品| 久久久久国产网址| 亚洲熟妇中文字幕五十中出| 国产精品一二三区在线看| 久久久久久久久久成人| 丰满人妻一区二区三区视频av| 夜夜夜夜夜久久久久| 99视频精品全部免费 在线| 18禁在线播放成人免费| 日本一本二区三区精品| 国产av在哪里看| 成人漫画全彩无遮挡| 最近最新中文字幕大全电影3| 久久久久久九九精品二区国产| 卡戴珊不雅视频在线播放| 日日摸夜夜添夜夜添小说| 黄片wwwwww| 一进一出抽搐动态| 国产精品亚洲一级av第二区| 亚洲最大成人av| 又爽又黄a免费视频| 亚洲真实伦在线观看| 亚洲精品久久国产高清桃花| a级一级毛片免费在线观看| 国产伦精品一区二区三区视频9| 亚洲国产色片| 国产白丝娇喘喷水9色精品| 久久久精品94久久精品| 亚洲精品乱码久久久v下载方式| 午夜福利高清视频| av卡一久久| 国产亚洲精品综合一区在线观看| 国产精品久久久久久精品电影| 少妇人妻精品综合一区二区 | 亚洲国产精品国产精品| 国产乱人视频| 天天躁日日操中文字幕| 赤兔流量卡办理| 国产精品国产三级国产av玫瑰| 干丝袜人妻中文字幕| 悠悠久久av| 日韩av在线大香蕉| 成人一区二区视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 在线a可以看的网站| 深爱激情五月婷婷| 大又大粗又爽又黄少妇毛片口| 国产精品av视频在线免费观看| 成人性生交大片免费视频hd| 亚洲av二区三区四区| 日韩成人av中文字幕在线观看 | 国国产精品蜜臀av免费| 国产乱人偷精品视频| 美女大奶头视频| 最近2019中文字幕mv第一页| 亚洲人成网站在线播放欧美日韩| 最近中文字幕高清免费大全6| 精品久久久久久久人妻蜜臀av| 午夜a级毛片| 秋霞在线观看毛片| 欧美日韩一区二区视频在线观看视频在线 | 我的女老师完整版在线观看| 性欧美人与动物交配| 国产免费一级a男人的天堂| 六月丁香七月| 伦精品一区二区三区| 成年女人永久免费观看视频| 日韩欧美国产在线观看| 国产一区二区三区av在线 | 可以在线观看的亚洲视频| 欧美绝顶高潮抽搐喷水| 男女边吃奶边做爰视频| 中文字幕久久专区| 热99在线观看视频| 国产精品一二三区在线看| 亚洲成人久久性| 亚洲欧美日韩卡通动漫| 午夜视频国产福利| 91麻豆精品激情在线观看国产| 蜜桃亚洲精品一区二区三区| 97碰自拍视频| 禁无遮挡网站| 天天一区二区日本电影三级| 变态另类丝袜制服| 国产v大片淫在线免费观看| 久久久久久九九精品二区国产| 夜夜爽天天搞| 中文字幕av在线有码专区| 日韩制服骚丝袜av| 性欧美人与动物交配| 国产乱人视频| 国产成人影院久久av| 伦理电影大哥的女人| 性色avwww在线观看| 国产真实乱freesex| 变态另类成人亚洲欧美熟女| 久久精品影院6| 老女人水多毛片| 久久久成人免费电影| 中国美女看黄片| 伊人久久精品亚洲午夜| 91久久精品国产一区二区三区| 午夜精品一区二区三区免费看| 99久久中文字幕三级久久日本| 亚洲久久久久久中文字幕| 国产 一区精品| 热99re8久久精品国产| 三级国产精品欧美在线观看| videossex国产| 免费无遮挡裸体视频| 校园人妻丝袜中文字幕| 尾随美女入室| 久久久精品欧美日韩精品| 成人一区二区视频在线观看| 亚洲自偷自拍三级| 国产色爽女视频免费观看| 九九久久精品国产亚洲av麻豆| 成人亚洲精品av一区二区| 亚洲欧美日韩无卡精品| 国产黄色视频一区二区在线观看 | 给我免费播放毛片高清在线观看| 黄色日韩在线| 欧美在线一区亚洲| 伊人久久精品亚洲午夜| 在线观看美女被高潮喷水网站| 亚洲性久久影院| 综合色av麻豆| 国产精品人妻久久久久久| 熟女电影av网| 国产精品爽爽va在线观看网站| 免费观看精品视频网站| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| 欧美bdsm另类| 老司机午夜福利在线观看视频| 久久人人精品亚洲av| 非洲黑人性xxxx精品又粗又长| 免费在线观看影片大全网站| 1024手机看黄色片| 不卡一级毛片| 成人av在线播放网站| 免费在线观看影片大全网站| 国产精品嫩草影院av在线观看| 亚洲专区国产一区二区| 亚洲成人中文字幕在线播放| 久久人妻av系列| 婷婷色综合大香蕉| 一进一出抽搐动态| 国产成人精品久久久久久| 给我免费播放毛片高清在线观看| 国产高清不卡午夜福利| 亚洲婷婷狠狠爱综合网| 伦理电影大哥的女人| 成人午夜高清在线视频| 欧美激情国产日韩精品一区| a级毛片a级免费在线| 人人妻人人看人人澡| 免费电影在线观看免费观看| 亚洲精品乱码久久久v下载方式| www日本黄色视频网| 精品无人区乱码1区二区| 日本免费一区二区三区高清不卡| 麻豆国产av国片精品| 久久久久久九九精品二区国产| av女优亚洲男人天堂| 永久网站在线| 午夜精品在线福利| 久久精品91蜜桃| 一级毛片我不卡| 免费观看在线日韩| 亚洲va在线va天堂va国产| .国产精品久久| 午夜福利在线观看吧| 一级黄色大片毛片| 欧美xxxx黑人xx丫x性爽| 黄色日韩在线| 国产精品久久久久久av不卡| 日韩成人伦理影院| 国产精品精品国产色婷婷| 久久久欧美国产精品| 久久久精品欧美日韩精品| 亚洲精华国产精华液的使用体验 | 五月伊人婷婷丁香| 国产v大片淫在线免费观看| 色噜噜av男人的天堂激情| 一级av片app| 搡老岳熟女国产| 九九爱精品视频在线观看| 热99在线观看视频| 精品人妻一区二区三区麻豆 | 男女视频在线观看网站免费| 黄色一级大片看看| 99久久无色码亚洲精品果冻| 三级经典国产精品| 久久久久国内视频| 日本黄色视频三级网站网址| 国产精品嫩草影院av在线观看| 麻豆久久精品国产亚洲av| 国产美女午夜福利| 欧美+亚洲+日韩+国产| 久久久久性生活片| 蜜臀久久99精品久久宅男| 国产精品久久久久久久久免| 午夜a级毛片| 直男gayav资源| 国产精品一及| 中文字幕久久专区| 麻豆国产97在线/欧美| 日韩在线高清观看一区二区三区| 亚洲最大成人av| 日本熟妇午夜| 18禁在线无遮挡免费观看视频 | 欧美激情在线99| 在线a可以看的网站| 午夜免费激情av| 悠悠久久av| 一级黄片播放器| 欧美又色又爽又黄视频| 日本黄色片子视频| avwww免费| 直男gayav资源| 免费高清视频大片| 国产激情偷乱视频一区二区| 91久久精品国产一区二区三区| 欧美人与善性xxx| 一区二区三区四区激情视频 | 熟女电影av网| 好男人在线观看高清免费视频| 男人舔奶头视频| 蜜桃亚洲精品一区二区三区| 99热网站在线观看| 国产探花在线观看一区二区| 我要看日韩黄色一级片| 日韩,欧美,国产一区二区三区 | 又粗又爽又猛毛片免费看| 精品无人区乱码1区二区| av国产免费在线观看| 欧美高清性xxxxhd video| 精品乱码久久久久久99久播| 日韩一区二区视频免费看| 中国国产av一级| 国内精品一区二区在线观看| 亚洲va在线va天堂va国产| 国产女主播在线喷水免费视频网站 | 国产精品女同一区二区软件| 成人高潮视频无遮挡免费网站|