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

    高強(qiáng)化柴油機(jī)鑄鐵缸蓋承載特性研究

    2017-03-08 05:22:01景國(guó)璽張樹勇付文清許春光文洋李鵬衛(wèi)軍朝
    車用發(fā)動(dòng)機(jī) 2017年1期
    關(guān)鍵詞:缸蓋熱應(yīng)力火力

    景國(guó)璽, 張樹勇, 付文清, 許春光, 文洋, 李鵬, 衛(wèi)軍朝

    (1. 中國(guó)北方發(fā)動(dòng)機(jī)研究所(天津), 天津 300400; 2. 中國(guó)兵器科學(xué)研究院, 北京 100089)

    高強(qiáng)化柴油機(jī)鑄鐵缸蓋承載特性研究

    景國(guó)璽1, 張樹勇1, 付文清2, 許春光1, 文洋1, 李鵬1, 衛(wèi)軍朝1

    (1. 中國(guó)北方發(fā)動(dòng)機(jī)研究所(天津), 天津 300400; 2. 中國(guó)兵器科學(xué)研究院, 北京 100089)

    研究了某柴油機(jī)鑄鐵缸蓋結(jié)構(gòu)溫度和應(yīng)力隨時(shí)間的變化行為,明確熱-機(jī)械載荷作用下缸蓋的承載規(guī)律,為缸蓋壽命預(yù)測(cè)模型建立提供依據(jù)。首先建立了缸蓋有限元分析模型,利用實(shí)測(cè)溫度和應(yīng)力數(shù)據(jù)對(duì)模型進(jìn)行了標(biāo)定,進(jìn)而基于該模型計(jì)算了標(biāo)定工況及怠速與標(biāo)定工況交替變化條件下缸蓋溫度和應(yīng)力的變化規(guī)律。仿真結(jié)果表明:缸蓋最高溫度和最大熱應(yīng)力出現(xiàn)在鼻梁區(qū)域;在工作循環(huán)內(nèi),缸蓋火力面溫度波動(dòng)幅值在30 ℃以內(nèi),由此引起的波動(dòng)熱應(yīng)力相對(duì)定常熱應(yīng)力較小,但相對(duì)高頻氣體應(yīng)力較大,故在火力面高周疲勞校核過程中必須考慮其帶來的影響;在怠速和標(biāo)定轉(zhuǎn)速交替變化工況下,鼻梁區(qū)載荷變化最明顯,溫度與應(yīng)力呈現(xiàn)反相狀態(tài),應(yīng)力幅值較高,該區(qū)域易發(fā)生低周疲勞損傷。

    柴油機(jī); 氣缸蓋; 仿真; 承載特性; 溫度; 應(yīng)力

    缸蓋是發(fā)動(dòng)機(jī)中結(jié)構(gòu)最復(fù)雜、承受機(jī)械載荷和熱載荷最嚴(yán)重的部件之一,其承載特性是影響缸蓋疲勞強(qiáng)度可靠性的重要因素[1-3]。缸蓋不僅承受螺栓預(yù)緊力的強(qiáng)約束作用,而且承受交變氣體力和熱應(yīng)力的耦合作用,同時(shí)由于發(fā)動(dòng)機(jī)服役工況范圍廣且多變,缸蓋承載規(guī)律極其復(fù)雜。研究表明,穩(wěn)態(tài)工況下螺栓預(yù)緊力、高頻溫度載荷與氣體力作用易導(dǎo)致缸蓋發(fā)生高周疲勞失效,失效部位主要出現(xiàn)在水腔隔板、氣道壁等區(qū)域。Mario Metzger的研究表明,工況大幅波動(dòng)(如停車—起動(dòng)—停車工況)易產(chǎn)生低頻波動(dòng)熱應(yīng)力,往往導(dǎo)致火力面鼻梁區(qū)發(fā)生低周疲勞失效[3]。目前,隨著發(fā)動(dòng)機(jī)功率密度的不斷提升,缸蓋承受的熱機(jī)械負(fù)荷增加,疲勞強(qiáng)度是當(dāng)前設(shè)計(jì)關(guān)注的焦點(diǎn)[4-5],因此,本研究重點(diǎn)探討了缸蓋溫度和應(yīng)力隨時(shí)間和載荷的變化規(guī)律,為缸蓋疲勞壽命預(yù)測(cè)模型建立提供依據(jù),為材料研制提供約束邊界。

    1 研究思路

    缸蓋必須同時(shí)滿足高周疲勞和低周疲勞的設(shè)計(jì)要求。高周疲勞強(qiáng)度按照無限壽命設(shè)計(jì),低周疲勞壽命與發(fā)動(dòng)機(jī)服役工況、服役壽命及用途相關(guān)。目前,主要基于發(fā)動(dòng)機(jī)臺(tái)架耐久性試驗(yàn)來考核發(fā)動(dòng)機(jī)結(jié)構(gòu)件疲勞特性,圖1示出發(fā)動(dòng)機(jī)典型臺(tái)架耐久規(guī)范的簡(jiǎn)化或等效示意,試驗(yàn)工況主要由標(biāo)定轉(zhuǎn)速穩(wěn)態(tài)工況及怠速與標(biāo)定交變工況兩部分組成,分別考核結(jié)構(gòu)件高周疲勞(HCF)和低周疲勞(LCF)性能。

    在圖1試驗(yàn)工況下,缸蓋主要承受以下作用力。1)預(yù)裝配應(yīng)力σf:主要由螺栓預(yù)緊力和氣門導(dǎo)管過盈引起;2)氣體應(yīng)力σg:主要由缸內(nèi)高頻波動(dòng)的氣體壓力引起;3)定常熱應(yīng)力σst:發(fā)動(dòng)機(jī)未工作前,結(jié)構(gòu)溫度較低且均勻,在高周疲勞考核穩(wěn)定工況(如標(biāo)定轉(zhuǎn)速工況)下,結(jié)構(gòu)溫度升高,由于熱慣性效應(yīng)溫度基本保持不變,由此形成的熱應(yīng)力即為定常熱應(yīng)力,通常該應(yīng)力已考慮預(yù)裝配應(yīng)力σf;4)高頻波動(dòng)熱應(yīng)力σht:在穩(wěn)態(tài)工況一個(gè)工作循環(huán)內(nèi),燃燒室氣體溫度、壓力和換熱邊界隨曲軸轉(zhuǎn)角高頻波動(dòng),引起缸蓋火力面表面溫度瞬態(tài)波動(dòng),并由此產(chǎn)生了高頻波動(dòng)的熱應(yīng)力;5)低頻波動(dòng)熱應(yīng)力σlt:在低周熱機(jī)疲勞考核工況下,缸蓋溫度隨著工況變化而發(fā)生大幅度波動(dòng),缸蓋處于急劇加熱和冷卻狀態(tài),由此產(chǎn)生了較大的低頻熱應(yīng)力。

    本研究主要開展以下工作:1)研究標(biāo)定轉(zhuǎn)速點(diǎn)穩(wěn)態(tài)工況下預(yù)裝配應(yīng)力σf、氣體應(yīng)力σg、定常熱應(yīng)力σst和高頻波動(dòng)熱應(yīng)力σht的分布特性和對(duì)缸蓋綜合應(yīng)力的影響規(guī)律;2)研究變工況條件下低頻波動(dòng)熱應(yīng)力σlt、溫度隨時(shí)間的變化關(guān)系;3)基于溫度和應(yīng)力分析結(jié)果,并結(jié)合發(fā)動(dòng)機(jī)臺(tái)架耐久規(guī)范進(jìn)一步構(gòu)建缸蓋疲勞壽命預(yù)測(cè)模型。

    圖1 發(fā)動(dòng)機(jī)臺(tái)架疲勞耐久試驗(yàn)考核工況示意

    主要研究思路見圖2。首先建立缸蓋有限元分析模型,利用實(shí)測(cè)溫度和應(yīng)力數(shù)據(jù)對(duì)模型進(jìn)行標(biāo)定,然后基于該模型計(jì)算標(biāo)定轉(zhuǎn)速工況穩(wěn)態(tài)條件下缸蓋溫度、應(yīng)力應(yīng)變行為,在此基礎(chǔ)上進(jìn)一步研究怠速工況與標(biāo)定工況交替變化條件下缸蓋載荷的變化特性。有限元網(wǎng)格模型見圖3,有限元模型包括缸蓋、進(jìn)排氣門、氣門導(dǎo)管、缸蓋螺栓、氣缸墊和氣缸體等部件,其中對(duì)缸蓋火力面等關(guān)鍵區(qū)域進(jìn)行了網(wǎng)格細(xì)化,模型共包含792 653個(gè)節(jié)點(diǎn),432 675個(gè)單元。

    圖2 總體研究思路

    圖3 有限元網(wǎng)格模型

    2 缸蓋溫度變化特性研究

    2.1 換熱邊界條件確定

    2.1.1 火力面熱邊界條件

    氣缸蓋火力面與高溫燃?xì)庵苯咏佑|,換熱過程以穩(wěn)態(tài)對(duì)流換熱為主。本研究基于AVLBOOST軟件采用Woshini模型計(jì)算得到了一個(gè)工作循環(huán)中的瞬時(shí)傳熱系數(shù)和燃?xì)饨橘|(zhì)溫度曲線,并基于式(2)得到了時(shí)間均化的當(dāng)量傳熱系數(shù)和燃?xì)鉁囟?,分別為1 703W/(m2·K)和840 ℃。

    ,

    (1)

    (2)

    式中:hgm和Tgm分別為時(shí)間均化的當(dāng)量傳熱系數(shù)和燃?xì)鉁囟取?/p>

    結(jié)合4氣門燃燒室空間結(jié)構(gòu),將火力面劃分為ININ,EXEX和EXIN3個(gè)區(qū)域,分別定義換熱邊界。各區(qū)傳熱系數(shù)沿著半徑方向的分布規(guī)律見圖4,圖中橫坐標(biāo)為實(shí)際位置相對(duì)缸徑的比值,縱坐標(biāo)為局部傳熱系數(shù)與當(dāng)量平均傳熱系數(shù)的比值。其中,ININ區(qū)域傳熱系數(shù)最小,EXEX區(qū)域傳熱系數(shù)最大,火力面熱邊界條件滿足式(3)要求。最后通過軟件編程實(shí)現(xiàn)了對(duì)火力面網(wǎng)格熱邊界條件的映射。

    圖4 火力面熱邊界條件分布規(guī)律

    (3)

    式中:r為氣缸半徑;b(r)為半徑r處的當(dāng)量傳熱系數(shù)。

    2.1.2 冷卻水腔表面熱邊界

    由于該發(fā)動(dòng)機(jī)采用并聯(lián)冷卻水腔設(shè)計(jì),各氣缸和單體式氣缸蓋采用單獨(dú)的進(jìn)出水口,因此僅對(duì)單缸冷卻水腔進(jìn)行分析。采用K-e模型對(duì)冷卻水進(jìn)行三維流動(dòng)數(shù)值模擬,入口邊界采用實(shí)測(cè)質(zhì)量流量(0.831 3kg/s),出口壓力為100kPa,壁面采用固定壁面邊界。通過計(jì)算獲得了缸體和缸蓋冷卻水腔表面?zhèn)鳠嵯禂?shù),并對(duì)有限元網(wǎng)格進(jìn)行了換熱邊界映射,其中缸蓋冷卻水腔映射前后傳熱系數(shù)見圖5。

    圖5 水腔表面?zhèn)鳠嵯禂?shù)

    2.1.3 其他區(qū)域熱邊界條件

    進(jìn)氣門與氣門座圈之間、排氣門與氣門座圈之間及氣門桿與氣門導(dǎo)管之間建立了接觸模型,接觸面間隙傳熱系數(shù)分別選取為10 000W/(m2·K),5 500W/(m2·K)和400W/(m2·K)[4]。其他區(qū)域換熱邊界采用經(jīng)驗(yàn)值。

    2.2 溫度場(chǎng)試驗(yàn)與模擬分析

    采用熱電偶方法進(jìn)行了缸蓋溫度場(chǎng)測(cè)試,熱電偶通過火力面鉆孔引出,測(cè)點(diǎn)布置見圖6。由于火力面輔助引線鉆孔直徑較大,阻礙火力面熱量向水側(cè)傳遞,對(duì)測(cè)溫結(jié)果有很大影響,因此,先基于試驗(yàn)測(cè)試用缸蓋結(jié)構(gòu)按照試驗(yàn)條件邊界進(jìn)行溫度場(chǎng)分析并與實(shí)測(cè)結(jié)果對(duì)比,確定火力面?zhèn)鳠嵯禂?shù)分布規(guī)律,進(jìn)而用于實(shí)際缸蓋溫度場(chǎng)預(yù)測(cè)。發(fā)動(dòng)機(jī)臺(tái)架試驗(yàn)工況條件下缸蓋溫度場(chǎng)測(cè)試和預(yù)測(cè)結(jié)果對(duì)比見圖7,由圖可知,實(shí)測(cè)最高溫度點(diǎn)位于兩個(gè)排氣門鼻梁區(qū),基于試驗(yàn)用缸蓋的多數(shù)測(cè)點(diǎn)實(shí)測(cè)值和計(jì)算值偏差在10 ℃以內(nèi),在相同熱邊界條件下實(shí)際缸蓋溫度偏低,峰值點(diǎn)溫度降低約40 ℃。

    圖6 熱電偶測(cè)點(diǎn)布置

    圖7 溫度場(chǎng)模擬結(jié)果與仿真結(jié)果對(duì)比

    在不考慮缸內(nèi)燃?xì)鉁囟群蛡鳠嵯禂?shù)隨時(shí)間變化對(duì)缸蓋整體溫度場(chǎng)影響的情況下,標(biāo)定工況點(diǎn)缸蓋溫度場(chǎng)模擬結(jié)果見圖8,高于300 ℃的區(qū)域主要位于與燃?xì)饨佑|的缸蓋底板區(qū)域,最高溫度為490 ℃,出現(xiàn)在兩排氣門鼻梁區(qū)域靠近氣門座孔邊緣處,缸蓋熱負(fù)荷極高。

    圖8 缸蓋溫度場(chǎng)分布

    2.3 工作循環(huán)內(nèi)溫度波動(dòng)分析

    溫度循環(huán)變化的燃?xì)饧捌渑c缸蓋的傳熱作用引起缸蓋火力面溫度的波動(dòng),研究表明[5]該溫度波動(dòng)向缸蓋內(nèi)部傳播,到達(dá)1~2mm時(shí)波動(dòng)幅度趨于零, 形成瞬變溫度層。三維瞬態(tài)溫度場(chǎng)可由下列微分方程式表達(dá):

    (4)

    初始條件

    (5)

    邊界條件

    (6)

    式中:T為溫度;τ為時(shí)間;k為導(dǎo)熱系數(shù);ρ為密度;c為比熱容;TW為壁面溫度;Tf為燃?xì)猱?dāng)量溫度;a為對(duì)流傳熱系數(shù)。

    矩陣常微分方程可描述為

    (7)

    本研究將前面穩(wěn)態(tài)溫度場(chǎng)計(jì)算結(jié)果作為瞬態(tài)熱分析的近似初始條件,進(jìn)而按照瞬態(tài)邊界計(jì)算了缸蓋在一個(gè)工作循環(huán)內(nèi)的溫度波動(dòng)。計(jì)算得到的缸蓋溫度波動(dòng)幅度見圖9,溫度波動(dòng)出現(xiàn)在靠近火力面約2mm的薄層內(nèi),標(biāo)定工況下溫度最大波動(dòng)幅度為26 ℃。計(jì)算得到的兩個(gè)排氣門之間特征區(qū)域溫度隨深度的變化關(guān)系見圖10,由圖可知,隨著深度增加,溫度波動(dòng)幅度逐漸減小,達(dá)到1mm后基本無波動(dòng)。

    圖9 缸蓋溫度場(chǎng)波動(dòng)幅

    圖10 火力面考察點(diǎn)溫度與深度的變化情況

    3 穩(wěn)態(tài)工況下載荷特性分析

    3.1 熱-機(jī)械耦合應(yīng)力分析

    為研究熱-機(jī)械載荷對(duì)缸蓋應(yīng)力的影響,利用有限元軟件計(jì)算了以下幾種工況:1)預(yù)裝配載荷(F),主要考慮氣門導(dǎo)管過盈裝配,加載缸蓋螺栓預(yù)緊力(130kN)和噴油器壓塊處螺栓預(yù)緊力(6 250N);2)預(yù)裝配載荷+熱載荷(F+T);3)預(yù)裝配載荷+最高燃燒壓力載荷(F+P);4)預(yù)裝配載荷+熱載荷+最高燃燒壓力載荷(F+T+P)。

    為驗(yàn)證仿真分析模型,在預(yù)裝配載荷條件下對(duì)缸蓋頂部和靠近螺栓側(cè)壁等區(qū)域進(jìn)行了應(yīng)力測(cè)試。與仿真值對(duì)比見圖11,由圖可見,多數(shù)測(cè)點(diǎn)應(yīng)力偏差較小,有限元仿真模型合理可信。

    圖11 預(yù)裝配載荷下應(yīng)力仿真與實(shí)測(cè)對(duì)比

    通過有限元應(yīng)力計(jì)算得到了預(yù)裝配載荷、熱載荷和氣體力載荷單獨(dú)作用及耦合作用下缸蓋應(yīng)力場(chǎng),結(jié)合溫度、應(yīng)力及疲勞計(jì)算結(jié)果,從缸蓋頂板、水腔、進(jìn)排氣道壁、底板冷卻水鉆孔和底板分別選取若干個(gè)特征點(diǎn)或危險(xiǎn)點(diǎn)作為應(yīng)力考察點(diǎn)(見圖12)。圖13示出各考察點(diǎn)應(yīng)力對(duì)比情況,圖14示出預(yù)裝配載荷、熱載荷、氣體力載荷對(duì)3種耦合作用下Mises應(yīng)力的貢獻(xiàn)系數(shù)。由圖13和圖14可知,預(yù)裝配載荷F對(duì)頂板、水腔和進(jìn)排氣道區(qū)域多數(shù)考察點(diǎn)的應(yīng)力貢獻(xiàn)較大,并占主導(dǎo)地位;熱載荷T作用下,由于底板溫度較高,底板向燃燒室方向膨脹,對(duì)頂部考察點(diǎn)耦合應(yīng)力有釋放緩解作用,而對(duì)底板冷卻鉆孔產(chǎn)生了較高的拉應(yīng)力,對(duì)鼻梁區(qū)產(chǎn)生了較高的壓應(yīng)力,熱載荷對(duì)底板考察點(diǎn)耦合應(yīng)力占主導(dǎo)作用。

    圖12 缸蓋應(yīng)力分析考察點(diǎn)

    圖13 不同載荷下考察點(diǎn)應(yīng)力

    圖14 不同載荷對(duì)耦合應(yīng)力的影響

    對(duì)于缸蓋結(jié)構(gòu)高周疲勞來說,預(yù)裝配載荷和熱載荷主要影響平均應(yīng)力,而氣體力載荷主要影響應(yīng)力幅。圖15示出考慮材料彈塑性影響修正后各考察點(diǎn)的應(yīng)力疲勞極限圖,由圖可見,載荷呈明顯非對(duì)稱分布,其中底板冷卻水鉆孔處部分考察點(diǎn)負(fù)載較高。

    圖15 各考察點(diǎn)應(yīng)力疲勞極限圖

    3.2 溫度波動(dòng)對(duì)熱機(jī)耦合應(yīng)力的影響

    在預(yù)裝配載荷和穩(wěn)態(tài)定常熱載荷的基礎(chǔ)上,考慮氣體力、高頻波動(dòng)熱載荷對(duì)缸蓋耦合應(yīng)力的影響,通過計(jì)算得到了其隨曲軸轉(zhuǎn)角的變化關(guān)系。圖16示出工作循環(huán)內(nèi)瞬態(tài)溫度波動(dòng)引起的高頻熱應(yīng)力的最大波動(dòng)幅值,由圖可知,應(yīng)力波動(dòng)區(qū)域與溫度瞬態(tài)變化直接相關(guān),高頻波動(dòng)熱應(yīng)力僅存在于火力面薄層區(qū)域內(nèi)[6]。取兩個(gè)排氣門之間鼻梁區(qū)D3考察點(diǎn)為研究對(duì)象,沿徑向正應(yīng)力的變化情況見圖17。由圖可見,氣體力作用導(dǎo)致壓應(yīng)力幅值進(jìn)一步增加,高頻波動(dòng)熱應(yīng)力相位與氣體力應(yīng)力明顯不同且數(shù)值較大,綜合應(yīng)力波動(dòng)幅為19.4MPa,是只考慮氣體力作用情況下應(yīng)力的1.9倍,可見,高頻波動(dòng)熱應(yīng)力對(duì)火力面區(qū)域結(jié)構(gòu)表面高周疲勞影響較大。

    圖16 熱應(yīng)力波動(dòng)幅

    圖17 考察點(diǎn)D3應(yīng)力波動(dòng)曲線

    4 變工況載荷特性分析

    按照發(fā)動(dòng)機(jī)熱沖擊規(guī)范,發(fā)動(dòng)機(jī)在圖1所示怠速和標(biāo)定轉(zhuǎn)速工況交替運(yùn)行,各工況運(yùn)行時(shí)間為3min,缸蓋溫度和應(yīng)力發(fā)生大幅度變化。在前文穩(wěn)態(tài)溫度場(chǎng)計(jì)算的基礎(chǔ)上,進(jìn)一步考慮材料溫度非線性特性,計(jì)算了變工況下缸蓋溫度和應(yīng)力隨時(shí)間的變化歷程(見圖18)。由圖18可知,在加熱過程中缸蓋底板特征點(diǎn)溫度急速增加,至60s后趨于穩(wěn)定,在加熱過程前60s內(nèi)平均加熱速率約為7.1 ℃/s,在冷卻過程前60s內(nèi)平均冷卻速率約為6.6 ℃/s;鼻梁區(qū)觀察點(diǎn)D2處應(yīng)力與溫度變化基本呈現(xiàn)反相狀態(tài),即最高溫度對(duì)應(yīng)最大壓應(yīng)力,而遠(yuǎn)離火力面的其他考察點(diǎn)溫度較低,應(yīng)力也較小。由此可見,除火力面鼻梁之外的其他區(qū)域溫度水平很低,最大應(yīng)力遠(yuǎn)小于材料屈服強(qiáng)度,而火力面鼻梁區(qū)結(jié)構(gòu)負(fù)載最高,為低周疲勞重點(diǎn)關(guān)注部位。

    圖18 變工況下溫度和應(yīng)力隨時(shí)間的變化曲線

    5 疲勞損傷與壽命預(yù)測(cè)模型建立

    在圖1所示試驗(yàn)工況下,缸蓋部件產(chǎn)生的總損傷可以表示為

    (8)

    式中:DTotal為材料發(fā)生的總損傷;DHCF為高周疲勞試驗(yàn)考核階段產(chǎn)生的總損傷;DLCF為低周疲勞試驗(yàn)考核階段所產(chǎn)生的總損傷。

    基于前文分析可知,式(8)中高周疲勞試驗(yàn)考核階段產(chǎn)生的總損傷DHCF主要由高頻波動(dòng)的氣體力和溫度載荷引起,可以描述為

    (9)

    式中:DF,HCF為高周高溫機(jī)械疲勞損傷;DC,HCF為蠕變損傷;kt為工作循環(huán)內(nèi)溫度高頻波動(dòng)熱應(yīng)力影響系數(shù),與缸蓋空間位置有關(guān),遠(yuǎn)離火力面溫度瞬態(tài)波動(dòng)層區(qū)域該系數(shù)為1,火力面區(qū)域可通過計(jì)算確定。

    在低周疲勞試驗(yàn)考核階段,缸蓋主要承受高低周疲勞載荷復(fù)合作用,材料產(chǎn)生的總損傷可以認(rèn)為是由于溫度大幅波動(dòng)而產(chǎn)生的低周熱機(jī)疲勞(TMF)損傷、高頻氣體力及熱載荷波動(dòng)引起高周疲勞損傷的疊加。其中,根據(jù)Sehitoglu理論[7],考慮疲勞與蠕變、疲勞與氧化之間的交互作用,熱機(jī)疲勞損傷DTMF主要包括疲勞損傷DF,TMF、氧化損傷DO,TMF和蠕變損傷DC,TMF三部分。因此,該階段產(chǎn)生的總損傷和低周熱機(jī)疲勞損傷可分別描述如下:

    (10)

    DTMF=DF,TMF+DO,TMF+DC,TMF。

    (11)

    由式(8)至式(11)可得,試驗(yàn)工況下缸蓋部件產(chǎn)生總損傷為

    (12)

    依照線性損傷法則,在該服役工況下缸蓋構(gòu)件的疲勞壽命方程[7]描述如下:

    (13)

    6 結(jié)論

    a) 在穩(wěn)態(tài)轉(zhuǎn)速工況下,在缸蓋火力面約1mm的薄層內(nèi)溫度存在明顯的瞬態(tài)波動(dòng),而其他區(qū)域溫度基本不變;工作循環(huán)內(nèi)瞬態(tài)溫度波動(dòng)引起的高頻熱應(yīng)力與溫度瞬態(tài)變化直接相關(guān),高頻波動(dòng)熱應(yīng)力僅存在于火力面薄層區(qū),氣體力作用導(dǎo)致該區(qū)域壓應(yīng)力進(jìn)一步增加,故高周疲勞強(qiáng)度評(píng)估時(shí)必須考慮高頻熱應(yīng)力的影響;

    b) 在怠速和標(biāo)定轉(zhuǎn)速交替變工況運(yùn)行下,缸蓋底板區(qū)域溫度發(fā)生大幅度變化,在加熱過程中缸蓋底板特征點(diǎn)溫度急速增加至60s后趨于穩(wěn)定,在加熱過程前60s內(nèi)平均加熱速率約為7.1 ℃/s,在冷卻過程前60s內(nèi)平均冷卻速率約為6.6 ℃/s;鼻梁區(qū)應(yīng)力與溫度變化基本呈現(xiàn)反相狀態(tài),而遠(yuǎn)離火力面的其他考察點(diǎn)溫度和應(yīng)力水平較低,鼻梁區(qū)為低周疲勞重點(diǎn)關(guān)注部位;

    c) 基于對(duì)缸蓋溫度和應(yīng)力變化特征的研究,提出了基于臺(tái)架服役工況的缸蓋疲勞損傷壽命預(yù)測(cè)模型,為后續(xù)壽命預(yù)測(cè)工作奠定了基礎(chǔ)。

    [1]ThomasHamm,MartinRebbert,Hermann-JosefEcker,etal.CylinderHeadDesignforHighPeakFiringPressures[C].SAEPaper2008-01-1196.

    [2]HazratMdAli,MasjukiHassan,KalamMd,etal.DeterminationofTimeVariant1D-3DTemperatureandHeatTransferDistri-butionInsidetheCoolingJacketofaSIEngineCoolingSystemafterKey-Off[C].SAEPaper2012-01-0642.

    [3]MarioMetzger,MarcLeidenfrost,EwaldWerner,etal.LifetimePredictionofEN-GJV450CastIronCylinderHeads[C].SAEPaper2014-01-9047.

    [4]HalfpennyA,AndersonR,LinX.IsothermalandThermo-MechanicalFatigueofAutomotiveComponents[C].SAEPaper2015-01-0548.

    [5]MetzgerM,SeifertT.Amechanism-basedmodelforLCF/HCFandTMF/HCFlifeprediction:Multiaxialformulation,finite-elementimplementationandapplicationtocastiron[J].TechnischeMechanik,2012,32:435-445.

    [6] 原彥鵬, 張衛(wèi)正, 程曉果,等.高強(qiáng)化內(nèi)燃機(jī)活塞瞬態(tài)溫度場(chǎng)分布規(guī)律研究[J].內(nèi)燃機(jī)工程,2005,26(4):35-38.

    [7]SehitogluH.,BoismierDA.Thermo-mechanicalfatigueofMar-M247:Part2Lifeprediction[J].TRANS.ASMEJ.Engng.MaterTechnol.,1990,112:80-89.

    [編輯: 李建新]

    Loading Characteristics of Cast Iron Cylinder Head for High-strengthened Diesel Engine

    JING Guoxi1, ZHANG Shuyong1, FU Wenqing2, XU Chunguang1, WEN Yang1, LI Peng1, WEI Junchao1

    (1. China North Engine Research Institute(Tianjin), Tianjin 300400, China;2. Ordnance Science and Research Academy of China, Beijing 100089, China)

    The changing behavior of temperature and stress for a cast iron cylinder head was investigated. Then the loading characteristics under thermal and mechanical loads were determined, which can provide the reference for building the life prediction model of cylinder head. The head finite element model was first built and calibrated with temperature and stress data and was further used for predicting the temperature and stress varying characteristics under steady working condition and altering condition between idle and rated speed. The simulation results show that the maximum temperature and thermal stress locate near the exhaust valve bridge region on the firing deck. The amplitude of temperature fluctuation in the firing deck region is less than 30℃, but it will still lead to thermal stress. The produced thermal stress is negligible compared with the steady thermal stress, but important to have to be considered in high cycle fatigue evaluation compared with high frequency gas stress. The load varies obviously in firing deck region under the altering condition between idle and rated speed. The temperature and stress are in out of phase, the stress amplitude is very high, and therefore the firing deck region generates the low cycle fatigue easily.

    diesel engine; cylinder head; simulation; loading characteristic; temperature; stress

    2016-01-26;

    2016-11-02

    景國(guó)璽(1984—),男,副研究員,博士,研究方向?yàn)閮?nèi)燃機(jī)現(xiàn)代設(shè)計(jì)方法、疲勞強(qiáng)度設(shè)計(jì)等;okjgx@163.com。

    10.3969/j.issn.1001-2222.2017.01.008

    TK413.2

    B

    1001-2222(2017)01-0042-07

    猜你喜歡
    缸蓋熱應(yīng)力火力
    WNS型鍋爐煙管管端熱應(yīng)力裂紋原因分析
    火力全開
    火力全開! 廣州上半年20條村改造,投入超800億!
    模擬缸蓋工藝在當(dāng)代內(nèi)燃機(jī)制造業(yè)中的應(yīng)用
    高性能連體缸蓋關(guān)鍵鑄造工藝研究及應(yīng)用
    《火力與指揮控制》投稿須知
    采用單元基光滑點(diǎn)插值法的高溫管道熱應(yīng)力分析
    缸蓋氣門導(dǎo)管座圈的加工及改進(jìn)
    鑄造缺陷對(duì)A356 鋁合金缸蓋疲勞性能的影響
    上海金屬(2015年5期)2015-11-29 01:13:43
    基于流熱固耦合的核電蒸汽發(fā)生器傳熱管熱應(yīng)力數(shù)值模擬
    综合色丁香网| videos熟女内射| 欧美老熟妇乱子伦牲交| 男女啪啪激烈高潮av片| 国产亚洲最大av| 精品少妇内射三级| av电影中文网址| 搡老乐熟女国产| freevideosex欧美| 赤兔流量卡办理| 久久av网站| 看免费成人av毛片| 久久鲁丝午夜福利片| 久久毛片免费看一区二区三区| 国产淫语在线视频| 国产av一区二区精品久久| 欧美老熟妇乱子伦牲交| 久久国产精品大桥未久av| 国产成人精品久久久久久| 一级片免费观看大全| 国产精品无大码| 免费观看a级毛片全部| 亚洲精品,欧美精品| 视频在线观看一区二区三区| 亚洲精品视频女| 男的添女的下面高潮视频| 国产免费福利视频在线观看| 亚洲欧美一区二区三区黑人 | 国产毛片在线视频| 黄色一级大片看看| 亚洲成人av在线免费| 欧美激情高清一区二区三区 | 91午夜精品亚洲一区二区三区| 国产精品 国内视频| 亚洲av成人精品一二三区| 国产精品欧美亚洲77777| 老司机影院成人| 99久久人妻综合| 大香蕉久久网| 久久精品久久久久久噜噜老黄| 观看美女的网站| 韩国高清视频一区二区三区| 不卡视频在线观看欧美| 国产熟女午夜一区二区三区| 久久97久久精品| 香蕉精品网在线| 亚洲久久久国产精品| 欧美 日韩 精品 国产| 亚洲国产精品国产精品| 肉色欧美久久久久久久蜜桃| 一级片'在线观看视频| 国产精品熟女久久久久浪| 午夜福利视频在线观看免费| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕色久视频| 水蜜桃什么品种好| 久久精品国产亚洲av涩爱| 日本午夜av视频| 日本爱情动作片www.在线观看| 亚洲精品中文字幕在线视频| 人人妻人人澡人人爽人人夜夜| 又大又黄又爽视频免费| 亚洲精品国产av成人精品| 久久久久久人人人人人| 国产有黄有色有爽视频| 午夜久久久在线观看| 亚洲欧美中文字幕日韩二区| 黄片小视频在线播放| 欧美日韩成人在线一区二区| 亚洲av电影在线观看一区二区三区| 麻豆乱淫一区二区| 大香蕉久久网| 大码成人一级视频| 亚洲欧美日韩另类电影网站| 看免费成人av毛片| 国产精品亚洲av一区麻豆 | 国产亚洲精品第一综合不卡| 香蕉丝袜av| 女人精品久久久久毛片| 水蜜桃什么品种好| 免费看不卡的av| 国产黄频视频在线观看| 观看美女的网站| 18禁动态无遮挡网站| 女人久久www免费人成看片| 亚洲欧洲国产日韩| 99久久精品国产国产毛片| 亚洲欧美色中文字幕在线| 欧美xxⅹ黑人| 国产亚洲午夜精品一区二区久久| 性色av一级| 亚洲成av片中文字幕在线观看 | 下体分泌物呈黄色| 亚洲,欧美,日韩| 只有这里有精品99| 亚洲婷婷狠狠爱综合网| 新久久久久国产一级毛片| 国产av码专区亚洲av| 国产男人的电影天堂91| 精品少妇一区二区三区视频日本电影 | 国产爽快片一区二区三区| 国产一区二区激情短视频 | 国产午夜精品一二区理论片| 美女主播在线视频| 两个人看的免费小视频| 在线亚洲精品国产二区图片欧美| 欧美xxⅹ黑人| 黑人猛操日本美女一级片| 日韩中文字幕欧美一区二区 | 国产午夜精品一二区理论片| 国产极品粉嫩免费观看在线| 日韩免费高清中文字幕av| 一级爰片在线观看| 国产精品 国内视频| 成人18禁高潮啪啪吃奶动态图| 婷婷色av中文字幕| 精品国产一区二区三区久久久樱花| 日韩大片免费观看网站| 国产 精品1| 亚洲视频免费观看视频| 国产精品秋霞免费鲁丝片| 春色校园在线视频观看| 亚洲国产看品久久| 两性夫妻黄色片| 观看美女的网站| 亚洲美女搞黄在线观看| 男的添女的下面高潮视频| 青春草亚洲视频在线观看| 黄色 视频免费看| 免费大片黄手机在线观看| 国产麻豆69| 在现免费观看毛片| 精品少妇黑人巨大在线播放| 日韩中文字幕视频在线看片| 国产亚洲一区二区精品| 人体艺术视频欧美日本| 精品一区在线观看国产| 黑人巨大精品欧美一区二区蜜桃| 国产精品 国内视频| 在线观看美女被高潮喷水网站| 熟妇人妻不卡中文字幕| 母亲3免费完整高清在线观看 | 日韩一区二区视频免费看| 2022亚洲国产成人精品| 午夜福利视频精品| 亚洲国产毛片av蜜桃av| 午夜av观看不卡| 精品卡一卡二卡四卡免费| 国产又爽黄色视频| 久久99精品国语久久久| 91精品国产国语对白视频| 一个人免费看片子| 黄片无遮挡物在线观看| 亚洲四区av| 亚洲精华国产精华液的使用体验| 国产成人91sexporn| 老司机影院成人| 尾随美女入室| 中国三级夫妇交换| 最近的中文字幕免费完整| 在线精品无人区一区二区三| 一级片'在线观看视频| 在线观看www视频免费| 国产一区二区激情短视频 | 亚洲,欧美,日韩| 极品人妻少妇av视频| a级毛片黄视频| 人人妻人人爽人人添夜夜欢视频| 日韩不卡一区二区三区视频在线| 一边摸一边做爽爽视频免费| 国产有黄有色有爽视频| 亚洲欧洲精品一区二区精品久久久 | a 毛片基地| 精品少妇黑人巨大在线播放| 亚洲精品国产一区二区精华液| 韩国高清视频一区二区三区| 亚洲av日韩在线播放| 国产精品女同一区二区软件| 丝袜喷水一区| 在线观看免费日韩欧美大片| 热re99久久国产66热| 建设人人有责人人尽责人人享有的| 亚洲国产日韩一区二区| 男女边吃奶边做爰视频| 亚洲图色成人| 日韩一本色道免费dvd| 国产1区2区3区精品| 巨乳人妻的诱惑在线观看| 汤姆久久久久久久影院中文字幕| 国产免费一区二区三区四区乱码| 人人妻人人爽人人添夜夜欢视频| 国产成人精品婷婷| 亚洲av在线观看美女高潮| 黄频高清免费视频| 91在线精品国自产拍蜜月| 午夜福利视频在线观看免费| 少妇被粗大猛烈的视频| 麻豆乱淫一区二区| 成人毛片60女人毛片免费| 在线观看美女被高潮喷水网站| 老女人水多毛片| 超色免费av| 亚洲国产最新在线播放| 国产不卡av网站在线观看| 日韩中字成人| 国产日韩欧美视频二区| 满18在线观看网站| 大片免费播放器 马上看| 成年动漫av网址| 最黄视频免费看| 五月天丁香电影| 亚洲久久久国产精品| 日韩中文字幕视频在线看片| 亚洲图色成人| 成人黄色视频免费在线看| 亚洲综合精品二区| 黄色配什么色好看| 男女国产视频网站| 精品亚洲成a人片在线观看| 18在线观看网站| 王馨瑶露胸无遮挡在线观看| 亚洲国产日韩一区二区| 亚洲av电影在线观看一区二区三区| 高清在线视频一区二区三区| 国产精品国产三级专区第一集| 美国免费a级毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 在线免费观看不下载黄p国产| 十八禁网站网址无遮挡| 大香蕉久久成人网| 亚洲一码二码三码区别大吗| 亚洲成人av在线免费| 一级毛片 在线播放| 午夜91福利影院| 丁香六月天网| 免费高清在线观看日韩| 久久久久久伊人网av| 久久午夜福利片| 亚洲精品一区蜜桃| 精品一品国产午夜福利视频| 91午夜精品亚洲一区二区三区| 制服人妻中文乱码| 蜜桃国产av成人99| 亚洲国产日韩一区二区| 精品一区二区三区四区五区乱码 | h视频一区二区三区| 成年美女黄网站色视频大全免费| 男女高潮啪啪啪动态图| 国产无遮挡羞羞视频在线观看| 日韩,欧美,国产一区二区三区| 国产精品亚洲av一区麻豆 | 国产精品久久久久久久久免| 超色免费av| 欧美国产精品va在线观看不卡| 亚洲国产日韩一区二区| 久久久久久久久免费视频了| 午夜免费观看性视频| 亚洲av.av天堂| 热re99久久精品国产66热6| 亚洲一码二码三码区别大吗| 99久久精品国产国产毛片| 捣出白浆h1v1| 欧美日韩视频高清一区二区三区二| 国产xxxxx性猛交| 亚洲欧美成人综合另类久久久| 青春草国产在线视频| 久久99一区二区三区| 亚洲精品av麻豆狂野| 蜜桃在线观看..| 欧美精品国产亚洲| 波多野结衣av一区二区av| 欧美日韩视频精品一区| 黄网站色视频无遮挡免费观看| 国产成人精品久久久久久| 精品少妇内射三级| 伦理电影大哥的女人| 国产人伦9x9x在线观看 | 久久午夜综合久久蜜桃| 午夜福利一区二区在线看| 一本色道久久久久久精品综合| 亚洲精品国产av蜜桃| freevideosex欧美| 少妇 在线观看| 中文字幕人妻丝袜制服| 一区二区三区乱码不卡18| 不卡视频在线观看欧美| 观看美女的网站| 少妇被粗大的猛进出69影院| 午夜久久久在线观看| 老女人水多毛片| 国产极品粉嫩免费观看在线| 亚洲国产成人一精品久久久| 国产熟女午夜一区二区三区| 丝袜在线中文字幕| 波多野结衣一区麻豆| 亚洲精品成人av观看孕妇| 成年美女黄网站色视频大全免费| 中文字幕精品免费在线观看视频| 精品国产乱码久久久久久男人| 卡戴珊不雅视频在线播放| 国产成人一区二区在线| 中国国产av一级| 国产精品香港三级国产av潘金莲 | 精品国产乱码久久久久久小说| 黄网站色视频无遮挡免费观看| av网站免费在线观看视频| 国产精品久久久久成人av| 国产白丝娇喘喷水9色精品| 亚洲情色 制服丝袜| 赤兔流量卡办理| 久久久国产一区二区| 新久久久久国产一级毛片| 国产欧美日韩综合在线一区二区| 三上悠亚av全集在线观看| 国产精品成人在线| 亚洲成人一二三区av| freevideosex欧美| 国产乱来视频区| 成人毛片a级毛片在线播放| 国产精品久久久久成人av| 国产高清不卡午夜福利| 性色avwww在线观看| 国产一区二区激情短视频 | 免费观看av网站的网址| 国语对白做爰xxxⅹ性视频网站| 女人久久www免费人成看片| 精品少妇久久久久久888优播| 久久亚洲国产成人精品v| www.熟女人妻精品国产| 久久久精品免费免费高清| 午夜激情久久久久久久| 精品一区二区三卡| 69精品国产乱码久久久| 天天躁日日躁夜夜躁夜夜| 国产 精品1| 亚洲精品中文字幕在线视频| 久久鲁丝午夜福利片| 又大又黄又爽视频免费| 丰满少妇做爰视频| 久久精品国产亚洲av天美| 国产女主播在线喷水免费视频网站| 黄色怎么调成土黄色| 伊人久久大香线蕉亚洲五| 美女脱内裤让男人舔精品视频| 中文字幕色久视频| 秋霞在线观看毛片| 国产精品一区二区在线观看99| 夫妻午夜视频| 三上悠亚av全集在线观看| 成人毛片a级毛片在线播放| 国产精品久久久久久av不卡| 精品视频人人做人人爽| 美女主播在线视频| 中文乱码字字幕精品一区二区三区| 亚洲久久久国产精品| 精品午夜福利在线看| 欧美亚洲日本最大视频资源| 天堂中文最新版在线下载| 黄片播放在线免费| 亚洲欧美色中文字幕在线| 国产精品无大码| 老司机亚洲免费影院| 国产欧美亚洲国产| 欧美中文综合在线视频| 日日啪夜夜爽| 久久精品亚洲av国产电影网| 久久国内精品自在自线图片| 老汉色av国产亚洲站长工具| 啦啦啦在线免费观看视频4| 欧美+日韩+精品| 国产日韩一区二区三区精品不卡| 亚洲成av片中文字幕在线观看 | 久热久热在线精品观看| 亚洲美女黄色视频免费看| 黄网站色视频无遮挡免费观看| 免费高清在线观看视频在线观看| 精品国产超薄肉色丝袜足j| 91aial.com中文字幕在线观看| 啦啦啦啦在线视频资源| 99香蕉大伊视频| 狠狠精品人妻久久久久久综合| 久久久久精品久久久久真实原创| 亚洲国产精品国产精品| 各种免费的搞黄视频| 亚洲国产毛片av蜜桃av| 午夜久久久在线观看| 精品人妻在线不人妻| 一级毛片电影观看| 深夜精品福利| 亚洲 欧美一区二区三区| 国产成人精品福利久久| 一级毛片我不卡| 最近中文字幕2019免费版| 婷婷色综合www| 国产成人精品久久久久久| 中文乱码字字幕精品一区二区三区| 国产女主播在线喷水免费视频网站| a 毛片基地| 亚洲精品自拍成人| 中文字幕av电影在线播放| 欧美日韩亚洲高清精品| 亚洲四区av| 丁香六月天网| 老鸭窝网址在线观看| 亚洲国产最新在线播放| 久久鲁丝午夜福利片| 亚洲国产av新网站| 欧美亚洲日本最大视频资源| 2022亚洲国产成人精品| 最近中文字幕高清免费大全6| 亚洲国产精品国产精品| 超色免费av| 99久国产av精品国产电影| 久久久久国产网址| av一本久久久久| 在线观看一区二区三区激情| 99香蕉大伊视频| 少妇熟女欧美另类| 免费观看av网站的网址| 午夜日韩欧美国产| 欧美激情 高清一区二区三区| 午夜激情久久久久久久| 97在线人人人人妻| 日韩伦理黄色片| 久久精品久久久久久久性| 日韩av不卡免费在线播放| 2021少妇久久久久久久久久久| 色94色欧美一区二区| 观看美女的网站| 国产午夜精品一二区理论片| 在线观看三级黄色| 三上悠亚av全集在线观看| 国产精品无大码| 免费黄色在线免费观看| 成人毛片60女人毛片免费| 亚洲精品国产av成人精品| 女性生殖器流出的白浆| 青青草视频在线视频观看| 80岁老熟妇乱子伦牲交| av又黄又爽大尺度在线免费看| 久久久久久人人人人人| 成人手机av| 少妇 在线观看| 美女福利国产在线| 日本黄色日本黄色录像| 王馨瑶露胸无遮挡在线观看| 边亲边吃奶的免费视频| 麻豆精品久久久久久蜜桃| 五月伊人婷婷丁香| 波多野结衣av一区二区av| 国产男女内射视频| 国产成人一区二区在线| 王馨瑶露胸无遮挡在线观看| 交换朋友夫妻互换小说| 亚洲精品中文字幕在线视频| 99热网站在线观看| 人妻 亚洲 视频| 黄色怎么调成土黄色| 欧美最新免费一区二区三区| 欧美日韩av久久| 免费高清在线观看视频在线观看| 免费黄网站久久成人精品| 久久久久视频综合| kizo精华| 精品少妇内射三级| 边亲边吃奶的免费视频| 欧美亚洲日本最大视频资源| 色吧在线观看| 18+在线观看网站| 久久99热这里只频精品6学生| 亚洲经典国产精华液单| 亚洲一区中文字幕在线| 国产精品无大码| 久久精品国产鲁丝片午夜精品| 又粗又硬又长又爽又黄的视频| 亚洲第一青青草原| 亚洲婷婷狠狠爱综合网| 老鸭窝网址在线观看| 欧美亚洲日本最大视频资源| 国产精品亚洲av一区麻豆 | 久久久久久久久久人人人人人人| 午夜激情av网站| 国产 一区精品| 日日啪夜夜爽| www.熟女人妻精品国产| av有码第一页| 国产精品久久久久久精品电影小说| 街头女战士在线观看网站| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品国产一区二区精华液| 99re6热这里在线精品视频| 亚洲婷婷狠狠爱综合网| 久久久亚洲精品成人影院| 老汉色∧v一级毛片| 国产成人免费观看mmmm| 宅男免费午夜| 天天躁日日躁夜夜躁夜夜| a级毛片在线看网站| 久久久久久久大尺度免费视频| 99久久中文字幕三级久久日本| 午夜免费男女啪啪视频观看| 啦啦啦视频在线资源免费观看| 久久精品国产亚洲av高清一级| 激情视频va一区二区三区| 亚洲美女视频黄频| 成年女人在线观看亚洲视频| 啦啦啦中文免费视频观看日本| 午夜福利在线观看免费完整高清在| 日韩一区二区三区影片| av片东京热男人的天堂| 宅男免费午夜| av电影中文网址| 汤姆久久久久久久影院中文字幕| 伊人久久国产一区二区| 国产视频首页在线观看| 久久精品人人爽人人爽视色| 九色亚洲精品在线播放| 久久精品久久久久久噜噜老黄| 国产精品香港三级国产av潘金莲 | 国产高清国产精品国产三级| 久久久国产精品麻豆| 97人妻天天添夜夜摸| 一区福利在线观看| 人人澡人人妻人| 欧美激情高清一区二区三区 | 亚洲在久久综合| 高清欧美精品videossex| 国产黄频视频在线观看| tube8黄色片| 黄片播放在线免费| 天堂中文最新版在线下载| 久久久久久免费高清国产稀缺| 十分钟在线观看高清视频www| 熟女av电影| 精品久久久精品久久久| 国产精品成人在线| 一级爰片在线观看| 久久99热这里只频精品6学生| 大话2 男鬼变身卡| 久久久久久久大尺度免费视频| 日韩av免费高清视频| 精品人妻熟女毛片av久久网站| 在现免费观看毛片| 少妇被粗大猛烈的视频| 人妻 亚洲 视频| 丝袜美足系列| 亚洲精品自拍成人| 极品少妇高潮喷水抽搐| 9色porny在线观看| av.在线天堂| 国产成人欧美| 啦啦啦在线观看免费高清www| 成年人免费黄色播放视频| 国产精品久久久久久精品电影小说| 99热国产这里只有精品6| 久热这里只有精品99| 毛片一级片免费看久久久久| 婷婷色麻豆天堂久久| 亚洲精品一区蜜桃| 国产一区二区 视频在线| 国产有黄有色有爽视频| 国产 精品1| 少妇人妻 视频| 久久狼人影院| 亚洲成人手机| 亚洲av中文av极速乱| 久久久久久久国产电影| av一本久久久久| 人体艺术视频欧美日本| 婷婷色av中文字幕| 国产男女内射视频| 国产精品一区二区在线观看99| 欧美成人午夜精品| 日韩伦理黄色片| 国产成人精品一,二区| 啦啦啦在线免费观看视频4| 久久免费观看电影| 99久久人妻综合| 美女xxoo啪啪120秒动态图| 日本-黄色视频高清免费观看| 国产精品久久久久久久久免| 国产精品三级大全| 精品国产一区二区三区久久久樱花| 在线观看www视频免费| 免费不卡的大黄色大毛片视频在线观看| 黄片小视频在线播放| 日本欧美国产在线视频| 久久久久精品久久久久真实原创| kizo精华| 在现免费观看毛片| 亚洲美女视频黄频| 人妻一区二区av| 国产精品久久久久久av不卡| 在线观看免费日韩欧美大片| 亚洲三级黄色毛片| 好男人视频免费观看在线| 两性夫妻黄色片| 亚洲精品久久成人aⅴ小说| 亚洲成人手机| 亚洲美女黄色视频免费看| 精品一品国产午夜福利视频| 午夜91福利影院| 成人黄色视频免费在线看| 久热这里只有精品99| 亚洲精品国产av成人精品| 日本av免费视频播放| 国产精品久久久av美女十八| 国产成人av激情在线播放| 亚洲一码二码三码区别大吗| 久久精品国产鲁丝片午夜精品| 黄色怎么调成土黄色| 91午夜精品亚洲一区二区三区| 亚洲国产毛片av蜜桃av| 美女福利国产在线| 街头女战士在线观看网站|